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Abstract of the Dissertation 

QCD Resummation Techniques 

by 

Tibor Kucs 
Doctor of Philosophy 
in 
Physics 

State University of New York 
at Stony Brook 
2004 

Dissertation Director: Prof. George Sterman 

The primary aim of high-energy QCD phenomenology is the de- 
termination of cross sections for particle coUisions. One of the 
fundamental properties of QCD, the asymptotic freedom, suggests 
that the coupling constant in this high-energy regime is small. Pro- 
vided that the coefficients of the perturbative expansion are small 
enough, the perturbation theory should give reliable results. How- 
ever, in many quantities of interest the smallness of the expansion 
coefficients is violated due to large logarithmic enhancements. In 
this case the perturbation series cannot be truncated at fixed order 
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and, instead, it must be resummed. The development of such re- 
summation algorithms is the main subject of the research presented 
in this thesis. 

In the first part, we propose a resummation technique applicable 
to the Regge limit, which is defined for elastic scattering as the 
region of large energies and small momentum transfer. We develop 
a new systematic procedure for this limit in perturbative QCD to 
arbitrary logarithmic order. The formalism relies on the IR struc- 
ture and the gauge symmetry of the theory. We identify leading 
regions in loop momentum space responsible for the singular struc- 
ture of the amplitudes and perform power counting to determine 
the strength of these divergences. Using a factorization procedure 
introduced by Sen, we derive a sum of convolutions in transverse 
momentum space over soft and jet functions, which approximate 
the amplitude up to power-suppressed corrections. A set of evolu- 
tion equations generalizing the BFKL equation and controlling the 
high energy behavior of the amplitudes to arbitrary logarithmic 
accuracy is derived. The general method is illustrated in the case 
of leading and next-to-leading logarithmic gluon reggeization and 
BFKL equation. We confirm the standard results at LL accuracy. 
At NLL order, we find an agreement with the reggeization conjec- 
ture up to two loops. However, starting at three loop order, we 
identify contributions violating the Regge ansatz. In addition, we 
calculate the evolution kernel determining the high-energy behav- 
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ior of the non-reggeized term in the scattering ampHtude. 

In the second part, we focus our attention to another intriguing 
problem of high-energy QCD, the resummation associated with soft 
radiation in dijet events which is comphcated by the presence of 
non-global logarithms. We introduce a set of correlations between 
energy flow and event shapes that are sensitive to the flow of color 
at short distances in jet events. These correlations are formulated 
for a general set of event shapes, which includes jet broadening and 
thrust as special cases. We illustrate the method for e+e~ dijet 
events, and calculate the correlation at leading logarithm in the 
energy flow and at next-to-leading-logarithm in the event shape. 
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Chapter 1 
Introduction 



The birth of modern theory of strong interactions dates back to the early 
sixties when Gell-Mann and Zweig, [1] proposed the quark model to describe 
the properties of hadrons. This model suggests that all strongly interacting 
objects are composed of more elementary constituents named quarks. The 
first experimental support of this idea came in the late sixties at SLAG in 
a Deeply Inelastic Scattering. After this, the top priority became the hunt 
for a theory describing the dynamics of quarks. The distinctive feature of 
quarks is that they carry internal quantum number called color. Each quark 
can be in three color states. In the early seventies, Pritzsch, Gell-Mann and 
Leutwyler, Ref . [2] , proposed a gauge theory describing the interaction between 
these quarks, called Quantum Ghromodynamcis (QGD). The mediator of this 
interaction, analogous to a photon in QED, is a gluon. At the classical level, 
QGD is invariant under the color SU{3) local gauge transformations. The 
quarks transform in the fundamental representation and the gluons in the 
adjoint representation of this group. The process of quantization introduces 
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a gauge-fixing and tlie ghost terms into the Lagrangian. It satisfies a global 
residual gauge symmetry, called BRS symmetry, Ref. [3]. This symmetry 
proves to be very powerful in deriving identities between Green functions which 
are exact, i.e. they hold to all orders in perturbation theory. 

The use of QCD as a perturbation theory is justified due to the property 
of asymptotic freedom discovered by 't Hooft, Gross, Wilczek and Pohtzer [4]. 
This property says that with increasing energies the coupling characterizing 
the interaction between the quarks and gluons decreases. 

At the Ultra- Violet (UV) spectrum of internal momenta, the 4-D field the- 
ory is plagued by infinities. They appear since the current theory is only an 
effective theory valid in a certain energy regime. At sufficiently high energies it 
must, presumably, be embedded into some more fundamental theory. Never- 
theless the theory is internally consistent, since we can remove these infinities 
by the process of renormalization, [5] . 

On the opposite end of momentum spectra, we encounter another type of 
divergences due to the Infra-Red (IR) region of soft momenta and coUinear 
momenta. These divergences occur only for Green functions with external 
particles on-shell, when at least one of the particles is massless. Their origin 
is due to the degeneracy of states occurring in the soft and coUinear limit, 
since we cannot possibly distinguish soft emissions and coUinear splittings from 
situations when these emissions and splittings are absent. This observation 
suggests, that although IR divergences appear on a diagram by diagram basis, 
they cancel in properly averaged quantities. Namely, one needs to sum over 
all indistinguishable states. Already in 1937, Bloch and Nordsieck, Ref. [6], 
showed that this is, indeed, the case in QED when the summation over final 
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states is performed. 

In QCD the situation is more complicated due to the self-couphng of gluons. 
In this case the KLN theorem, Ref. [7], which extends the summation over 
final state degeneracies to initial states as well, comes to our rescue. It is these 
quantities, which are free of IR divergences in QCD. From the experimental 
point of view this solution to the IR problem has a caveat, however. We can 
hardly expect to prepare our initial states in collision experiments to satisfy the 
conditions of the KLN theorem. Instead, we take one step further and invoke 
factorization theorems, Ref. [8]. It is exactly these theorems that enable us to 
turn perturbative QCD to a predictive calculational tool. 

Factorization theorems claim that it is possible to separate short and long 
distance physics in physical quantities. A cross section for hadronic colli- 
sion can be written as a convolution over functions describing long distance 
dynamics, which can be either Parton Distribution Functions (PDFs) or Frag- 
mentation Functions (FFs), and functions describing short distance physics, 
which are partonic cross sections. The former have a physical interpretation of 
probabilities to find partons inside hadrons in case of PDFs, or to find hadrons 
inside partons in the case of FFs. These quantities cannot be completely de- 
termined from perturbative QCD and need to be fitted to experimental data. 
However, perturbative QCD enables us to find the evolution of the density 
functions with energy scale, Ref. [9]. These evolution equations resum large 
logarithmic corrections in this energy scale. They are the consequence of the 
factorized form for the cross section and the renormalization group equation 
stating that the physical quantities cannot depend on the scale at which we 
make the separation of short and long distance dynamics. Actually, under 
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very general assumptions, we can claim that whenever there is a factorization, 
there is a corresponding resummation. We will encounter concrete examples 
of this statement in later chapters. The main feature of the distribution func- 
tions is that they are process independent. This suggests that after fitting long 
distance functions in one process, we can use them in certain other processes 
involving the same type of hadron. 

The second ingredient of the factorized cross section is the partonic cross 
section, which quantifies the interaction of the underlying partons in the hard 
process. This quantity is IR safe and calculable within perturbation theory 
provided the coefficients of the perturbative expansion are small. This is, how- 
ever, not always the case. The coefficients are usually enhanced at kinematic 
edges of phase space. The reason for this is easily understood. The cancellation 
of IR divergences happens between virtual and real corrections. The integra- 
tion in the virtual corrections spans all the energy scales: soft, hard and UV. 
The last one is removed by renormalization, so only the region between the soft 
and the hard scale remains. The phase space for real corrections depends on 
the kinematics considered. If we are completely inclusive then the integration 
region is the same as in the case of virtual corrections and there is a perfect 
cancellation between the two. However, once we impose kinematic constraints 
on the final state, then the cancellation between the two is incomplete and we 
can be left with large logarithms spoiling the perturbative expansion. In this 
case we need to resum the perturbation series. 

The development of resummation procedures for various processes is the 
main topic of this work. It consists of two parts. In the first part we propose 
a systematic resummation technique applicable in the Regge limit, Ref. [10] 
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- [11]. This limit concerns almost forward elastic scattering of particles. In 
the second part of this thesis we pursue resummation of soft radiation accom- 
panying final state jets in lepton collisions, Ref. [12]. 
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Part I 
Regge Resummation 
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Chapter 2 
The Method 



2.1 Introduction 

The study of semi-hard processes within the framework of gauge quantum 
field theories has a long history. For reviews see Refs. [13]- [15]. The defin- 
ing feature of such processes is that they involve two or more hard scales, 
compared to Aqcd, which are strongly ordered relative to each other. The 
perturbative expansions of scattering amplitudes for these processes must be 
resummed since they contain logarithmic enhancements due to large ratios of 
the scales involved. One of the most important examples is elastic 2 — > 2 
particle scattering in the Regge limit, s ^ |i| (with s and t the usual Mandel- 
stam variables). It is this process that we investigate in the next two chapters. 
We extend the techniques developed in Refs. [16] and [17] and devise a new 
systematic method for evaluation of QCD scattering amplitudes in the Regge 
limit to arbitrary logarithmic accuracy, Ref. [10]. 

The problem of the Regge limit in quantum field theory was first tackled in 
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the case of fermion exchange amphtude within QED in Ref. [18]. Here it was 
found that the positive signature amphtude takes a reggeized form up to the 
two loop level in Leading Logarithmic (LL) approximation. In Ref. [19] the 
calculations were extended to higher loops, and the imaginary part of the Next- 
to-Leading Logarithms (NLL) was also obtained. The analysis in Refs. [18] 
and [19] was performed in Feynman gauge. It was realized in Ref. [20] that a 
suitable choice of gauge can simphfy the class of diagrams contributing at LL. 
The common feature of all this work was the use of fixed order calculations. 
To verify that the pattern of low order calculations survives at higher orders, a 
method to demonstrate the Regge behavior of amplitudes to all orders is nec- 
essary. This analysis was provided by A. Sen in Ref. [16], in massive QED. Sen 
developed a systematic way to control the high energy behavior of fermion and 
photon exchange amplitudes to arbitrary logarithmic accuracy. The formalism 
relies heavily on the IR structure and gauge invariance of QED and provides 
a proof of the reggeization of a fermion at NLL to all orders in perturbation 
theory. 

The resummation of color singlet exchange amplitudes in non-abelian gauge 
theories in LL was achieved in the pioneering work of Ref. [21], where the 
reggeization of a gluon in LL was also demonstrated. The evolution equations 
resumming LL in the case of three gluon exchange was derived in Ref. [22] . In 
Ref. [17], n-gluon exchange amplitudes in QCD at LL level were studied and a 
set of evolution equations governing the high energy behavior of these ampli- 
tudes was obtained at LL. A different approach was undertaken in Ref. [23]. 
Here n ^ m amplitudes were studied in SU(2) Higgs model with spontaneous 
symmetry breaking. Starting with the tree level amplitudes, an iterative proce- 
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dure was developed, which generates a minimal set of terms in the perturbative 
expansion that have to be taken into account in order to satisfy the unitarity 
requirement of the theory. See also Ref. [24]. The extension of the BFKL 
formalism to NLL spanned over a decade. For a review see Ref. [25]. The 
building blocks of NLL BFKL are the emissions of two gluons or two quarks 
along the ladder, Ref. [26] , one loop corrections to the emission of a gluon along 
the ladder, Ref. [27], and the two loop gluon trajectory, Refs. [28], [29], [30] 
and [31]. The particular results were put together in Ref. [32]. In Ref. [33], 
the trajectory for the fermion at NLL was evaluated by taking the Regge limit 
of the explicit two loop partonic amplitudes, Ref. [34]. 

Besides the NLO perturbative corrections to the BFKL kernel a variety 
of approaches have been developed for unitarization corrections, Refs. [35- 
37], which extend the BFKL formalism by incorporating selected higher-order 
corrections. The procedure proposed in this work, Refs. [10] and [11], in a 
way, places these approaches in an even more general context. In principle, 
it makes it possible to find the scattering amplitudes to arbitrary logarithmic 
accuracy and to determine the evolution kernels to arbitrary fixed order in the 
coupling constant. The formalism contains all color structures and, of course, 
the construction of the amplitude to any given level requires the computation 
of the kernels and the solution of the relevant equations. 

The first part of the thesis is organized as follows. In Chapter |21 we develop 
the general algorithm. In Sec. 12.21 we discuss the kinematics of the partonic 
process under study and the gauge used. In Sec. 12.31 we identify the leading 
regions in internal momentum space, which produce logarithmic enhancements 
in the perturbation series. After identifying these regions, we perform power 
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counting to verify that the singularity structure of individual diagrams is at 
worst logarithmic. The leading regions lead to a factorized form for the ampli- 
tude (First Factorized Form). It consists of soft and jet functions, convoluted 
over soft loop momenta, which can still produce logarithms of s/|)f:|. In Sec. 
12.41 we study the properties of the jet functions appearing in the factorization 
formula for the amplitude. We show how the soft gluons can be factored from 
the jet functions. In Sec. 12. 51 we demonstrate how to express systematically the 
amplitude as a convolution in transverse momenta. In this form all the large 
logarithms are organized in jet functions and the soft transverse momenta in- 
tegrals do not introduce any logarithms of s/\t\ (Second Factorized Form). We 
derive evolution equations that enable us to control the high energy behavior 
of the scattering amplitudes. 

In Chapter El we illustrate the general method valid to all logarithmic 
accuracy in the case of LL and NLL in the amplitude and we examine the 
evolution equations at the same level. In Sec. we resum the amphtude at 
LL and we find the LL gluon Regge trajectory. In Sec. 13.21 we analyze the 
amplitude at NLL order. We confirm the LL BFKL evolution equation. In 
Sec. 13. 3[ we address the problem of NLL evolution equations and the gluon 
reggeization at this accuracy. We confirm the Regge hypothesis at two loop 
level. However, we identify contributions violating the Regge ansatz starting 
at three loop order. 

Some technical details are discussed in appendices lA.ll - IB. 21 The first 
appendix treats power counting for regions of integration space where internal 
loop momenta become much larger than the momentum transfer. In Appendix 
IA.2I we illustrate the origin of special vertices encountered due to the resum- 
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mation. In Appendix IA.3I we show a systematic expansion for the amphtude 
leading to the first factorized form. In Appendix IA.4I we list the Feynman 
rules used throughout the text. In Appendix IA.5I we demonstrate the origin 
of extra soft momenta configurations (Glauber region) which need to be con- 
sidered in the analysis of amplitudes in the Regge limit. In Appendix lB.il we 
study some symmetry properties of jet functions and finally in Appendix IB.2[ 
we give details on the derivation of the color octet NLL evolution equations. 

2.2 Kinematics and Gauge 

We analyze the amplitude for the elastic scattering of massless quarks 

q{pA, ta, Xa) + q'iPB, rs, Xb) q{pA - q, ri, Xi) + q'{pB + q, ^2, A2), (2.1) 

within the framework of perturbative QCD in the kinematic region s ^ —t 
{Regge limit), where s = {pa +PbY ^"^^ t = q^ are the usual Mandelstam vari- 
ables. We stress, however, that the results obtained below apply to arbitrary 
elastic two-to-two partonic process. We pick process 1)2.11) for concreteness 
only. The arguments in Eq. 1)2.11) label the momenta, pj, the colors, rj, and 
polarizations, Aj, for i = A,B,1,2, respectively, of the initial and final state 
quarks. We choose to work in the center-of-mass (cm.) where the momenta 
of the incoming quarks and the momentum transfer have the following com- 
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ponents ^ 



Pa 

Pb 
Q 



1,0-, 0, 



0+,Wx,0 



(0+,0-,gx). 



(2.2) 



Strictly speaking = ±|t|/v2s, so the components vanish in the Regge 
hmit only. 

In the color basis 



1 



2Nr 



(2.3) 



with Nf. the number of colors, we can view the amplitude for process ()2.H) as 
a two dimensional vector in color space 



A 



^8 



(2.4) 



where A-i and Ag, are defined by the expansion 



A 



l)rArB,ri r2 



+ As {b. 



•r-A rB, ri T-2 • 



(2.5) 



Since the amplitude is dimensionless and all particles are massless, its compo- 



"'^We use light-cone coordinates, v = [v^ ,v — ± v^)/^/2. 
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nents can depend, in general, on the following variables 



Ai = Aii —,—,as{fi^),e 



) 



for i = 1, 8, 



(2.6) 



where /z is a scale introduced by regularization. We use dimensional regu- 
larization in order to regulate both infrared (IR) and ultraviolet (UV) di- 
vergences with D = A — 2e the number of dimensions. Choosing the scale 
fi"^ = s, the strong coupling, ^^(/i), is small. However, in general, an in- 
dividual Feynman diagram contributing to the process ()2.H1 at r-loop order 
can give a contribution as singular as (s/t) al~^Hn'^^(—s/t). In Sec. 12.5.31 we 
will confirm that there is a cancellation of all terms proportional to the i-th 
logarithmic power for z = r + 1, . . . , 2r at order a^'"*"^ in the perturbative ex- 
pansion of the amplitude. Hence at r loops the amplitude is enhanced by a 
factor (s/t) al~^^ W{—s/t), at most. In order to get reliable results in pertur- 
bation theory we must, nevertheless, resum these large contributions. In the 
k-th non-leading logarithmic approximation one needs to resum all the terms 
proportional to (s/t) a^"*"^ luf^^ {~s/t), j = 0, . . . , k at r-loop level. 

We perform our analysis in the Coulomb gauge, where the propagator of a 
gluon with momentum k has the form 



—i 6, 



'ab 



Nap{k,k) 
k"^ + ie 



= —i 6, 



'ab 



k"^ + ie 



1 



( 



9al3 — 



kfj + k^ kp kct kfj 




in terms of the vector 



k = k — {k ■ 1]) 1], 



(2.8) 
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with 

an auxiliary four-vector defined in the partonic cm. frame. The numerator of 
the gluon propagator satisfies the following identities 

k ■ k 

k''Na/3{k,k) = 0. (2.10) 

The first equality in Eq. ()2.1()|1 is the statement that the nonphysical degrees of 
freedom do not propagate in this gauge. For use below, we list the components 
of the gluon propagator: 



N+~{k) = N~+{k) 
N++{k) = N-'{k) 
N^\k) = N'^{k) 



k+k- 



k-k' 

^ {k- -k+)k' 

2k -k 

k'y 



rC ' rC 



We note that these are symmetric functions under the transformation k^ — >■ 
—k^, except for the components A^^ * = A^* ^, which are antisjTumetric under 
this transformation. It was demonstrated in Ref. [38] that QCD is renormal- 
izable in Coulomb gauge, by considering a class of gauges which interpolates 
between the covariant (Landau) and the physical (Coulomb) gauge. 
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2.3 Leading Regions, Power Counting 

In order to resum the Regge logarithms, we need to identify the regions of 
integration in the loop momentum space that give rise to singularities in the 
hmit i/s — > 0. We follow the method developed in Refs. [39,40], which begins 
with the identification of the relevant regions in momentum space. 

2.3.1 Singular contributions and reduced diagrams 

The singular contributions of a Feynman integral come from the points in loop 
momentum space where the integrand becomes singular due to the vanishing 
of propagator denominators. However, in order to give a true singularity the 
integration variables must be trapped at such a singular point. Otherwise we 
can deform the integration contour away from the dangerous region. These 
singular points are called pinch singular points. They can be identified with 
the following regions of integration in momentum space, 

1. soft momenta, with scaling behavior k'^ ~ c ^/s for all components 

2. momenta coUinear to the momenta of the external particles, with scaling 
behavior 

~ a/s, k~ ~ A v^, ~ A^/^ ^/s for the particles moving in the + 
direction and 

k'^ ~ Xy/s, k~ ~ ^/s, ~ A^/^ ^/s for the particles moving in the — 
direction, 

3. so-called Glauber or Coulomb momenta, Ref. [41], with scaling behavior 
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a) b) c) 

Figure 2.1: The reduced diagrams a) and c) contributing to the amphtude. 
Diagram b) represents a decomposition of diagram a) for the purpose of power 
counting. 

~ 0"^ ^/s, \k_i_\ ~ (T ^/s, where A < a"^ < cr, and where the scahng 
factors A, a satisfy the strong ordering A ^ a ^ 1 (The origin of this 
region is illustrated in Appendix IA.5I ). 

4. hard momenta, having the scaling behavior /c^ ~ y/s for all components. 

The extra gauge denominators l/{k ■ k) originating from the numerators of 
the gluon propagator, Eq. ()2.7|) . do not alter the classification of the pinch 
singular points mentioned above. Actually, only the subsets 1 and 3 in the 
above classification can be produced due to the extra gauge denominators. 

With every pinch singular point, we may associate a reduced diagram, 
which is obtained from the original diagram by contracting all hard lines (sub- 
set 4) at the particular singular point. As shown in Refs. [39, 40, 42] the 
reduced diagram corresponding to a given pinch singular point must describe 
a real physical process, with each vertex of the reduced diagram representing 
a real space-time point. This physical interpretation suggests two types of 
reduced diagrams contributing to the process !^'2.1^ . shown in Fig. 12.11 
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The jet A{B) contains lines whose momenta represent motion in the + (— ) 
direction. The hnes included in the blob S' and the lines coming out of it are 
all soft (configurations 1 and 3 in the classification of loop momenta described 
above). These two oppositely moving (virtual) jets may interact through the 
exchange of soft lines, Fig. 12.1b . and/or they can meet at one or more space- 
time points, Fig. 12.1b . 

Having found the most general reduced diagrams giving the leading be- 
havior of the amplitude for process ()2.H) in the Regge limit, we can estimate 
the strength of the IR divergence of the integral near a given pinch singular 
point. First we restrict ourselves to cases involving subsets 1 and 2 from the 
classification of loop momenta above. To do so, we count powers in the scaling 
variables A and a. 

The scaling behavior of these loop momenta implies that every soft loop 
momentum contributes a factor cr^, every jet loop momentum gives rise to 
the power A^, every internal soft boson (fermion) line provides a contribu- 
tion (T^^ (c"^^) and every internal jet line (fermionic or bosonic) scales as A~^. 
In addition, there can be suppression factors arising from the numerators of 
the propagators associated with internal lines and from internal vertices. As 
pointed out in Ref . [39] , in physical gauges each three-point vertex connecting 
three jet lines is associated with a numerator factor that vanishes at least lin- 
early in the components of the transverse jet momenta, and therefore provides 
a suppression A^/^. 

We are now ready to estimate the power of divergence corresponding to 
the reduced diagrams describing our process. First we restrict ourselves to 
the case shown in Fig. 12.1k . As indicated schematically in Fig. 12.1b . we can 
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perform the power counting for the jets and for the soft part separately. All 
soft propagators and all soft loop momenta are included in the soft subdiagram 
S. The superficial degree of IR divergence of the reduced diagram R from Fig. 
12.1b and Fig. 12.1b can then be written as 



where the external lines and loops of S' are included in S. For uj{R) > 
the overall integral is finite, while uj{R) < corresponds to an IR divergent 
integral. When uj{R) = 0, the integral diverges logarithmically. Here we set 
A ~ 0" for power counting purposes. We come back to the effect of relaxing 
this condition in connection with a discussion of item 3, Glauber regions, in 
our list of singular momentum configurations. 

2.3.2 Power counting 

In this subsection, we consider the case when all vertices in a diagram are 
elementary only, that is, without contracted sub-diagrams carrying large loop 
momenta. In Appendix lA.ll we show that our conclusions are unchanged by 
contracted vertices. 

We perform the power counting for the soft part S first. Let f,b be the 
number of fermion, boson lines external to S' and let E = f+b. The superficial 
degree of divergence for S, found by summing powers of a, can be written 



uj{R) 



(2.12) 



uj{S) = 4(E - 2) - 26 - / + 2 + u{S'), 



(2.13) 
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where the first term is due to loop integrations linking S' to the jets, while 
the second and the third terms originate from propagators associated with 
the bosonic and fermionic lines, respectively, connecting the jets A, B and 
the soft part S' . The term +2 is introduced because we are resumming only 
leading power corrections proportional to s/t and therefore we exclude the 
overall factor s/t from the power counting. Since the lines entering S' are soft, 
we obtain the superficial degree of divergence for S' simply from dimensional 
analysis. It is given by 



Combining Eqs. ()2.13|) and ()2.14|) . the superficial degree of infrared divergence 
for the soft part 5* is then 



Before carrying out the jet power counting, we introduce some notation. 
Let Ea be the number of soft lines attached to jet A\ I is the total number 
of jet internal lines; Va is the number of a-point vertices connecting jet lines 
only; Wa has a meaning similar to Va-, with the difference that every vertex 
counted by Wa has at least one soft line attached to it. These are the vertices 
that connect the jet A to the soft part 5*. Finally, L denotes the number of 
loops internal to jet A. As noted above, we will perform the power counting 
for the case when the scaling factor for the soft momenta, a, is of the same 
order as the scaling factor for jet A momenta. When the scaling factors are 



uj{S') =4-6-3//2. 



(2.14) 



oj{S) = 6 + 3//2-2. 



(2.15) 
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different we encounter sub divergencies, which can be analyzed the same way 
as described below. We also assume that there are no internal and external 
ghost lines included in the jet function. Later we will discuss the effect of 
adding ghost lines. 

The superficial degree of divergence for jet A can now be expressed as 

uj{A) = 2L- I + V3/2. (2.16) 

The last term represents the suppression factor associated with the three point 
vertices. We denote the total number of vertices internal to jet A by 

^ = 5ZK + ^a)- (2.17) 

a 

Next we use the Euler identity relating the number of loops, internal lines and 
vertices of jet A 

L = I-v + l, (2.18) 
and the relation between the number of lines and the number of vertices 

2J + Ea + 2 = ^a(t;„ + w„). (2.19) 

a 

Using Eqs. ()2.16|) - ()2.19|) we arrive at the following form for the superficial 
degree of divergence for jet A 

u{A) = 1-{Ea + W3)/2 + 5^(« - A){v^ + w^)/2. (2.20) 

Q>5 
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Since every vertex counted by Wa connects at least one external soft line, we 
have the condition 

Ea>W3 + J2wc. (2.21) 

The equality holds when there is no vertex with two or more soft lines attached 
to it. Combining Eqs. (|2.2(Jj) - (j2.2H) we arrive at the following lower bound on 
the superficial degree of divergence for jet A: 

uj{A) > 1-E^ + J]^,;J2 + J](a-4)(t;„ + ti;„)/2. (2.22) 
The third and the last term in Eq. ()2.22j) are always positive or zero and hence 

u{A) >1-Ea. (2.23) 

A similar result holds for jet -B, and therefore the superficial degree of coUinear 
divergence for jets A and B is 

u{A) +uj{B) >2 - E, (2.24) 

with E = Ea + Eb as in Eq. ()2.13|) . Combining the results for soft and jet 
power counting, Eqs. (j2.15|) and (j2.24p . respectively in Eq. (j2.12|) . we finally 
obtain the superficial degree of IR divergence for the reduced diagram in Fig. 

uj{R) > f/2. (2.25) 

This condition says that we can have at worst logarithmic divergences, pro- 
vided no soft fermion lines are exchanged between the jets A and B. We 
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can therefore conclude that a reduced diagram from Fig. 12.1b containing ele- 
mentary vertices can give at worst logarithmic enhancements in perturbation 
theory. In order for the divergence to occur, the following set of conditions 
must be satisfied: 

1. There is an exchange of soft gluons between the jets A and B only, with 
no soft fermion lines attached to the jets. 

2. The jets A and B contain 3 and 4 point vertices only, see Eq. ()2.22j) . 

3. Soft gluons are connected to jets only through 3 point vertices, Eq. 
(I2.22|l . and at most one soft line is attached to each vertex inside the 
jets, Eq. ^rm . 

4. In the reasoning above we have assumed that there is no suppression 
factor associated with the vertices where soft and jet lines meet. In 
order for this to be true, the soft gluons must be connected to the jet 
A{B) lines via the +(— ) components of the vertices. 

Next we consider adding ghost lines to the jet functions. As we review in 
Appendix IA.41 the propagator for a ghost line with momentum k is propor- 
tional to l/{k -k). Hence every internal ghost line belonging to the jet gives 
a contribution which is power suppressed as 1/s. Since the numerator factors 
do not compensate for this suppression, we can immediately conclude that the 
jet functions cannot contain internal or external ghost lines at leading power. 

So far we have not taken into account the possibility when the soft loop 
momenta are pinched by the singularities of the jet lines. This situation allows 
different components of soft momenta to scale differently. For example, a 



22 



minus component of soft momentum can scale as the minus component of jet 
A momentum A, while the rest of the soft momentum components may scale 
as 0", where A ^ a ^ 1. The origin of these extra pinches is illustrated in 
Appendix IA.51 

Let us see what happens when we attach the ends of a gluon line with this 
extra pinch to jet A at one end and the soft subdiagram S at the other end. 
The integration volume for this soft loop momentum scales as Ao"^. The soft 
gluon denominator gives a factor cr~^. If this soft gluon is connected to the soft 
part at a 4-point vertex, there is no new denominator in the soft part. On the 
other hand, if the soft gluon is attached to the soft part via a 3-point vertex 
then the extra denominator including the numerator suppression factors scales 
as a"^. The new jet line scales as A~^ as long as the condition A^/^ > cr is 
obeyed; otherwise, we have the scaling for the extra jet line. For A^/^ > a 
the Glauber region produces logarithmic infrared divergence. When A^^^ < a, 
the overall scaling factor A/cr^ indicates power suppressed contribution. 

Let us now investigate another possibility, when the soft gluon connects jet 
A and jet B directly and its momentum is pinched by the singularities of the 
jet A and the jet B lines. Denoting the scaling factors of jet A and jet B as Xa 
and Xb, respectively, the integration volume provides the factor A^A^cr^ and 
the soft gluon denominator contributes the power a~^. The extra jet A and 
jet B denominators scale as A^^ and A^^, provided A^^ > a and A]/^ > a. For 
-^A B ^ ^ both extra jet denominators provide the scaling factor cr~^. When 

1 /2 

^AB ~ '^y power counting suggests logarithmically divergent integrals. 

We have therefore verified that when the softest component of a soft line 
satisfies the ordering o"^ < A < a, the Glauber (Coulomb) momenta produce 
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logarithmically IR divergent integrals and need to be taken into an account 
when identifying enhancements in perturbation series. The analysis demon- 
strated above for the case of one Glauber gluon can be extended to the situa- 
tion with arbitrary number of Glauber gluons. This follows from dimensional 
analysis, in a similar fashion as the treatment of purely soft loop momenta 
above. 

We conclude that the reduced diagram in Fig. 12.1b is at most logarithmi- 
cally IR divergent, modulo the factor s/\t\. The reduced diagram in Fig. 12.1b 
looses one small denominator compared to the reduced diagram in Fig. 12. lb 
and since we are working in physical gauge, this loss cannot be compensated 
by a large kinematical factor coming from the numerator. Hence the reduced 
diagram in Fig. 12.1b is power suppressed compared to the reduced diagram in 
Fig. 12.1b . and we do not need to consider it at leading power. 

Finally, let us discuss the scale of the soft momenta. In the case of soft 
exchange lines, each gluon propagator supplies a factor l/{a^s), which we 
want to keep at or below the order t in the leading power approximation. 
Thus the size of the scale is fixed to be cr ~ \/\t\/Vs- the case of soft 
lines which are attached to jet A or to jet B only, the scaling factor lies in the 
interval {^/\t\/y/s,l). In the case of Glauber momenta, we again need cr ~ 
A^^/ v^. Then the condition A^^^ > o", which is necessary for the logarithmic 
enhancement, implies that the scaling factors for + and — components of the 
Glauber (Coulomb) momenta can go down to \t\/s, the scale of the small 
components of jet momenta. Additionally, we should note that soft and jet 
sub-diagrams that do not carry the momentum transfer may approach the 
mass shell (A, a ^ 0). Such lines produce true infrared divergences, which 
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a) b) 

Figure 2.2: Jet A moving in the + direction (a) and jet B moving in the — 
direction (b). 

we assume are made finite by dimensional regularization to preserve the gauge 
properties that we will use below. The same power counting as above shows 
that these divergences are also at worst logarithmic. 



2.3.3 First factorized form 

The analysis of the previous subsection suggests the following decomposition 
of the leading reduced diagram from Fig. 12. lb . Let us denote the {n + 2)- 
point and (m + 2)-point Green functions, IPl in external soft gluon lines, 
corresponding to jet A, JIa)1\'.'.."iZ^Pa, V',h,---, K), Fig. and to jet B, 
J(B)ll['.'.lZ^PB,q,mPi^ ■ ■ ■ ^Pm), Fig. Od, respectively. The jet function jj^j 
('^(B)) ^^^'^ depends on the color of the incoming and outgoing partons r^, 
ri {tb, T2), as well as on their polarizations A^, Ai (As, A2), respectively. In 
order to avoid making the notation even more cumbersome we do not exhibit 
this dependence explicitly. In addition the dependence of jj^^^ and j''^^^ on 
the renormalization scale /i and the running coupling Q;s(/i) is understood. 
The jet functions also depend on the following parameters: the gauge fixing 
vector 7], Eq. ()2.9|1 . of the Coulomb gauge, the four momenta of the external 
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soft gluons attached to jet A (B), ki, . . . ,kn (pi, . . . ,Pm), and the Lorentz and 
color indices of the soft gluons attached to the jet A (B), /ii, . . . , oi, . . . , a„ 
(z/i, . . . , Um', bi, . . . , hm)- The momenta of the soft gluons attached to the jets 
A and B satisfy the constraints Yll=i = 1 ^j=i P:i ~ 1' 

According to the results of the power counting, the soft gluons couple to jet 
A via the minus components of their polarizations, and to jet B via the plus 
components of their polarizations. Therefore, only the following components 
survive in the leading power approximation 

JA^''''"^'"{pA,q,V,VB]ki,. . . ,kn) = ^n^^'j JiA)l\'.'..'Z^PA,(l,V',h,- ■ ■ ,k„), 

jir^^'"'^"^iPB,q,V,VA;Pl,- ■ ■,Pm) = (fl^A^ JIb)uI'.'.'.uZ^PB,Q,V-,Pi, ■ ■ ■ ,Pni), 

(2.26) 



.1=1 



where we have defined light-like momenta in the plus direction va = (1, 0, 0^) 
and in the minus direction vb = (0, 1, 0_l). We can now write the contribution 
to the reduced diagram in Fig. 12.1k . and hence to the amplitude for process 
fj2.1|) . in the form 



''n—l \ „ /m—1 



n,m-^ \i=l / \i=l / 

X V. ^A, Vb] ki,..., pi, . . . , p^) 

X /^^''■■"'-{pB,q,V,VA;Pu...,Pm), (2.27) 

where the sum over repeated color indices is understood. Corrections to Eq. 
(I2.27|l are suppressed by positive powers of t/s. The jet functions Ja,b are de- 
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fined in Eq. ()2.26p in the leading power accuracy. The internal loop momenta 
of the jets A, B and of the soft function S are integrated over. The soft function 
will, in general, include delta functions setting some of the momenta ki, . . . ,kn 
and color indices of jet function Ja to the momenta Pi, . . . ,Pm and 

to the color indices 61, . . . , &m of jet function Jb- The construction of the soft 
function S is described in Appendix IA.3I For a given Feynman diagram there 
exist many reduced diagrams of the type shown in Fig. 12.1k . and one has to be 
careful in systematically expanding this diagram into the terms that have the 
form of Eq. ()2.27|1 . This systematic method can be achieved using the "tulip- 
garden" formalism first introduced in Ref. [44] and used in a similar context 
in Ref. [16]. For convenience of the reader we summarize this procedure in 
Appendix I A. 31 

Let us now identify the potential sources of the enhancements in ln(s/|t|) of 
the amplitude given by Eq. ()2.27|) . If we integrate over the internal momenta of 
the jet functions then we can get ln((p^ ■ from Ja and ln((pB ■ v)^/^) 

from Jb- In addition, according to the results of the power-counting, Eq. 
(I2.23|l . we know that the jet function with n external soft gluons diverges as 
1/A"~^. After performing the integrals over the minus components of the ex- 
ternal soft gluon lines attached to jet A and over the plus components of the 
external soft gluons connected to jet B, these divergent factors are potentially 
converted into logarithms of ln((p^ ■ ri)^/\t\) and \n{{pB ■ vY/ltl), respectively. 
Our goal will be to separate the full amplitude into a convolution over param- 
eters that do not introduce any further logarithms of the form ln(s/|t|). This 
task will be achieved in Sec. 12.5.11 In the following section, we analyze the 
characteristics of the jet functions. 
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a) 



b) 



Figure 2.3: a) Decoupling oi a. K gluon from jet A. b) Leading contributions 
resulting from the attachment of a G gluon to jet A. 

2.4 The Jet Functions 

In this section we study the properties of the jet functions A, B given by Eq. 
()2.26|) since, as Eq. ()2.27|) suggests, they will play an essential role in later 
analysis. Since the methods for both jet functions are similar we restrict our 
analysis to jet A only; jet B can be worked out in the same way. In Sec. 
12.4.11 we examine the properties of jet A when the minus component of one 
of its external soft gluon momenta is of order \/\t\. In Sec. 12.4.21 we find the 
variation of jet A with respect to the gauge fixing vector r], and finally in Sec. 
12.4.31 we examine the dependence of jet A on the plus component of a soft 
gluon momentum attached to this jet. 

2.4.1 Decoupling of a soft gluon from a jet 

According to the results of power counting above, soft gluons attach to lines in 
jet A via the minus components of their polarization. Following the technique 
of Grammer and Yennie [45] we decompose the vertex at which the jth gluon 
is connected to jet A. We start with a trivial rewriting of J a in Eq. ()2.2(i|l 
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We now decompose the metric tensor into the form g'^'^ = K^'^{kj) + G^'^{kj) 
where for a gluon with momentum kj attached to jet A, K'^'^ and G'^'^ are 
defined by 



va ■ kj — ie 

G^'^ikj) = g^"' - K^'^ikj). (2.29) 

The K gluon carries scalar polarization. Since the jet A function has no 
internal tulip-garden subtractions (they are contained in the soft function S), 
we can use the Ward identities of the theory [46], which are readily derived 
from its underlying BRS symmetry [47], to decouple this gluon from the rest 
of the jet A after we sum over all possible insertions of the gluon. The result 
is 



Va ■ kj — ie 

n 

(-«^«r'"'"0 j^""^^"^-''--^ ■■■""(pa, q,VB,V, ki,...,ki + kj, . . . , kj, . . . k^). 



X 



(2.30) 



The notation a^- and kj indicates that the jet function does not depend 

on the color index aj and the momentum kj, because they have been factored 
out. In Eq. ()2.H()|1 . gs is the QCD coupling constant and f^'"-'"-^ are the 
structure constants of the SU{3) algebra. The pictorial representation of this 
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equation is shown in Fig. 12.3b . The arrow represents a scalar polarization and 
the double line stands for the eikonal line. The Feynman rules for the special 
vertices and the eikonal lines in Fig. 12.3b are listed in Appendix IA.4I Strictly 
speaking the right-hand side of Eq. (|2.3(J|) and Fig. 12.3b contain contributions 
involving external ghost lines. However, from the power counting arguments 
of Sec. 12.3.21 we know that when all lines inside of the jet are jet-like, the 
jet function can contain neither external nor internal ghost lines. Therefore 
Eq. ()2.30p is valid up to power suppressed corrections for this momentum 
configuration. 

The idea behind the K-G decomposition is that the contribution of the 
soft G gluon attached to the jet line in the leading power is proportional 
to v'^G^^v'^ = 0. In order to avoid this suppression, the G gluon must be 
attached to a soft line. The general reduced diagram corresponding to the G 
gluon attached to jet A is depicted in Fig. 12.3b . The lines coming out of S 
as well as the lines included in it are soft. The letter G next to the jth gluon 
in Fig. 12.3b reminds us that this gluon is a G-gluon attaching to jet J{A)^i via 
the G~^^{kj) vertex. 

The reasoning described above applies to the case when all components of 
soft momenta are of the same order. In the situation of Coulomb (Glauber) 
momenta, this picture is not valid anymore, since the large ratio k±/k~ com- 
ing from the component can compensate for the suppression due to the 
attachment of the G part to a jet A line via the transverse components of the 
vertex. 
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2.4.2 Variation of a jet function with respect to a gauge 
fixing vector 77 

In this subsection we find the variation of the jet function jj^"''' with respect 
to a gauge fixing vector rj. The motivation to do this can be easily under- 
stood. We consider the jet function with one soft gluon attached to it only, 
J^\pA,q,VB,f])- Let us define 

^a=Pa-V and (b = VVb- (2.31) 

In these terms, jet function J^^ can depend on the following kinematical 
combinations: Ja\pai(1iVbiV) = Ja\^a^ Pa ■ vb, Cb, t). Using the identity 
Pa - vb = 2^aC-B and the fact, that the dependence of Ja on the vector vb is 
introduced trivially via Eq. ()2.26p . we conclude that 

ji\pA,q,VB,il) = CbJ'a\U t). (2.32) 

Our aim is to resum the large logarithms of In(pJ^) that appear in the pertur- 
bative expansion of the jet A function. In order to do so, we shall derive an 
evolution equation for p\ dJ^^ /dp\. Since pa appears in combination with rj 
only, we can trace out the p~^ dependence of J^^ by tracing out its dependence 
on 77. This can be achieved by varying the gauge fixing vector 77. The idea 
goes back to Collins and Soper [44] and Sen [43] . We will generalize the result 
to in Sec. EE^ 

We consider a variation that corresponds to an infinitesimal Lorentz boost 
in a positive + direction with velocity 6f3. Thus, for the gauge fixing vector 
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a) b) 

Figure 2.4: The result of a variation of jet function J^"^ with respect to a 
gauge fixing vector. 

7] = (1,0,0,0) 2, Eq. (Q, the variation is: Sr] = fjS(3 = (0,0,0,1)5/3. It 
leaves invariant the norm r^^ = 1 to order 0{6(3). The precise relation between 
the variation of the jet A function with respect to p\ and 61]" is 

a 7(1) f) t(1) a t{1) f) t(1) 

dp\ 7]°' OQb 7]°' 

We have used the chain rule in the first equality and the simple relation 

(b dJ^^ /d(B = Ja^, following from Eq. ()2.32|) . in the second one. 

In order for Eq. (j2.33p to be useful, we need to know what the variation of 

jet A with respect to the gauge fixing vector 7] is. The result of this variation 

for J^"^ is shown in Fig. 12.41 It can be derived using either the formalism of 

the effective action, Ref. [48], or a diagrammatic approach first suggested in 

Ref. [44] and performed in axial gauge. We give an argument how Fig. 12.41 

arises in Appendix IA.21 Here we only note that the form of the diagrams in 

Fig. 12.41 is a direct consequence of a IPI nature of the jet functions. The 
^For the moment we use Cartesian coordinates. 
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explicit form of the boxed vertex 



-iS''{k) = -i{r]-kfi'^ + fi-kr]'^), (2.34) 

as well as of the circled vertex is given in Fig. IA.2I of Appendix IA.41 while 
their origin is demonstrated in Appendix IA.2I The dashed lines in Fig. 12.41 
represent ghosts, and these are also given in Fig. IA.2I of Appendix I A. 41 The 
four vectors rj, given in Eq. ()2.9|) . and 



appearing in Eq. ()2.34|) are defined in the partonic cm. frame, Eq. ()2.2|) . We 
list the components of Sfj,N^°'{k) 

kl-k^ 



S,{k)N^-ik) = k^[-f^±^ 
S,{k)N^\k) = (2.36) 



for later reference. 

In Fig. 12.41 we sum over all external gluons. This is indicated by the sum 
over i. In addition, we sum over all possible insertions of external soft gluons 
{ii, . . . ,in^} G {!,..., n}\{i}. This summation is denoted by the symbol vr. 
We note that at lowest order, with only a gluon i attached to the vertical 
blob in Fig. 12.4b . this vertical blob denotes the transverse tensor structure 
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depending on the momentum ki of this gluon 

t (k^g^^ - ktk^ . (2.37) 

It is labeled by a gluon line which is crossed by two vertical lines, Fig. IA.2I 
The ghost line connecting the boxed and the circled vertices in Fig. 12.4b can 
interact with jet A via the exchange of an arbitrary number of soft gluons. We 
do not show this possibility in Fig. 12.4b for brevity. 

Let us now examine what the important integration regions for a loop with 
momentum k in Fig. 12.4b are. The presence of the ghost line and of the 
nonlocal boxed vertex requires that in the leading power the loop momentum 
k must be soft. It can be neither collinear nor hard. This will enable us to 
factor the gluon with momentum k from the rest of the jet according to the 
procedure described in Sec. 12.4.11 

2.4.3 Dependence of a jet function on the plus compo- 
nent of a soft gluon's momentum attached to it 

In this subsection we want to find the leading regions of the object k+d/j'^/dk+. 
This information will be essential for the analysis pursued in the next sections. 
For a given diagram contributing to we can always label the internal loop 
momenta in such a way that the momentum kj flows along a continuous path 
connecting the vertices where the momentum kj enters and leaves the jet 
function When we apply the operation kjd/dkj on a particular graph 

corresponding to it only acts on the lines and vertices which form this 
path. The idea is illustrated in Fig. 12.5b . The gluon with momentum k at- 
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taches to jet A via the three-point vertex vi. Then the momentum k flows 
through the path containing the vertices vi,V2,V3 and the hues li,l2- The 
action of the operator k~^d/dk^ on a hne or vertex which carries jet-like mo- 
mentum gives a negligible contribution, since the + component of this lines 
momentum will be insensitive to k^. In order to get a non-negligible contribu- 
tion, the corresponding line must be soft. In Fig. 12.5b . lines li and I2 must be 
soft in order to get a non-suppressed contribution from the diagram after we 
apply the k^d/dk~^ operation on it. This, with the fact that the external soft 
gluons carry soft momenta, also implies that the lines ^3, . . . , le must be soft. 
This reasoning suggests that in general a typical contribution to kj'dJ^^^ /dkj' 
comes from the configurations shown in Fig. 12.5b . It can be represented as 

X /f {pA,q,V,VB;k[,...,k'^,). (2.38) 



The function j contains the contributions from the soft part S and from 
the gluons connecting the jet ■* and S in Fig. 12.5b . The jet function ■* 
has fewer loops than the original jet function J^""*. Now applying the operation 
kj'd/dk^ to Eq. ()2.38p . the operator kj'd/dkj' acts only to the function 
Hence we can write 



r) 

L+ ^ jin)ai ...a„ 



ki, ■ ■ ■ , kn', k^, . . . , k^i) " (p^, q, rj, Vb', ki, . . . , k^,). 

(2.39) 
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a) 



kj 
b) 



Figure 2.5: a) Momentum flow of tfie external soft gluon inside of jet A. b) 
Typical contribution to k'^dJ^^/dk'j'. 

We conclude that the contribution to k'^dJ^J^'^ /dk'^ can be expressed in terms 

(n') 

of jet functions J\ which have fewer loops than the original jet function. 



2.5 Factorization and Evolution Equations 

We are now ready to obtain evolution equations which will enable us to resum 
the large logarithms. First, in Sec. 12.5. H we will put Eq. (j2.27|) into what 
we call the second factorized form. Then, in Sec. 12.5.21 we derive the desired 
evolution equations. In Sec. I2.5.3[ we will show the cancellation of the double 
logarithms and finally in Sec. I2.5.4[ we demonstrate that the evolution equa- 
tions derived in Sec. 12. 5. 21 are sufficient to determine the high-energy behavior 
of the scattering amplitude. 
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2.5.1 Second factorized form 

The goal of this subsection is to rewrite Eq. ()2.27j) into the following form [16] 



^n~l \ /m—1 



n,m \j=l / \j=l 

X r^^ (pa, g,?7,fB; ki±, kn±; M) 

X rg"^'"-''"(pB,g,r7,t;A;Pi±,...,P^±;M), (2.40) 



where F^'' and are defined as the integrals of the jet functions J^"'' and 
over the minus and plus components, respectively, of their external soft 
momenta, with the remaining light-cone components of soft momenta set to 
zero, 



/ rM 

^^"^"••■'^"(PA,g,r?, 



"-i / I'M \ 

vb] ki±, . . . , kn±; M) = Y\yJ '^^^ J 
X JA^^'"'''"iPA,q,V,'VB;kix,. . . ,kn±,kt = 0, ...,/c+ = 0,k^,...,k-), 

X J^ir^^'"'''"{pB,q,V,VA]Pl±, ■ ■ ■ ,Pni±,PY = 0,...,p;; = 0,pt,. . . ,pt^)■ 

(2.41) 

In Eq. ()2.4()|1 . S' is a calculable function of its arguments and M is an arbitrary 
scale of the order ^/\t\. The functions F^^^ and S' depend individually on this 
scale, but the final result, of course, does not. Based on the discussion at the 
end of Sec. I2.3.3[ one can immediately recognize that all the large logarithms 
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are now contained in the functions Ta and Tb- The convolution of Ta, Tb and 
S' is over the transverse momenta of the exchanged soft gluons. Since these 
momenta are restricted to be of the order the integration over transverse 
momenta cannot introduce ln(s/|t|). This indicates that at leading logarithm 
approximation the factorized diagram with the exchange of one gluon only 
contributes. In general, when we consider a contribution to the amplitude at 
L = La + Lb + Ls' loop level, where La, Lb and Ls> is the number of loops 
in Ta, Tb and S', respectively, we can get L — Ls' logarithms of s/\t\ at most. 
Hence, the investigation of the s/t dependence of the full amplitude reduces 
to the study of the p\ and dependence of Ta and Tb, respectively. We 
formalize this statement at the end of Sec. 12.5.31 after we have proved that Ta 
(Tb) contains one logarithm of p\ (p^) per loop. 

Let us now show how we can systematically go from Eq. ()2.27|) to Eq. 
(imUl . We follow the method developed in Ref. [16]. We start from Eq. (IIITTI) 
and consider the k~ integrals over the jet function Ja for fixed kf, ki±: 




i=l 



{pA,q,V,VB;ki,...,kn), 



(2.42) 



where Ra is given by the soft function S and the jet function Jb, 




,{n,m) 

ai... an,bi... f), 



(g, T], va, f b; fci, . . . , fc„; Pi, . . . , Pm) 



X j: 



r(m) fei... b, 
B 



{pB,q,V,VA;Pl, ■ ■ ■,Pm)- 



(2.43) 
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We next use the following identity for Ra- ^ 

n-l 

RA{K,...,k-__,) = Ra{K = o,...,k-_, = o)l[e{M-\k-\) 

1=1 

n-l 

+ [RA{ki, . . . , k-, = 0, . . . , = 0) 

i=l 

- RAik^, K-1, K = o,..., A;-_i = 0) eiM - \k;\)] 

n-l 

X Yl e{M-\kj\). (2.44) 
j=i+i 

We have suppressed the dependence on the color indices and other possible 
arguments in Ra for brevity. The scale M can be arbitrary, but, as above, we 
take it to be of the order of \/\t\. The first term on the right hand side of Eq. 
()2.44|) has all k~ =0. The rest of the terms can be analyzed using the K-G 
decomposition discussed in Sec. 12.4.11 Consider the (z = 1) term, say, in the 
square bracket of Eq. ()2.44|1 inserted in Eq. ()2.42j) . Let us denote it Ai. In the 
region ^ M the integrand vanishes. On the other hand, for \k^\ ~ M 
we can use the K-G decomposition for the gluon with momentum ki. The 
contribution from the K part factorizes and the integral over the component 
has the form 

/AU- 
[i?l--(fcr,fc2- = 0,...,A;-_i = 0) 
Va ■ K,i 

- eiM - \k^\) R^-^-ik^ = 0, . . . , = 0)] 

n-l / „M n-l 

X Y.i'dsr'''" ll^kjJ^r"^''^—"ipA,q,V,VB; 

k2,...,ki + ki,...,kn)). (2.45) 

■^Recall that kn = q — (ki + . . . + fc„_i), so fc„ is not an independent momentum. 
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Eq. ()2.45|) is valid when all the lines inside the jet are jet-like. In that case 
the contributions from the ghosts are power suppressed. The contribution 
corresponding to a G gluon comes from the region of integration shown in 
Fig. 12.3b . It can be expressed in the form of Eq. (j2.42j) involving some jj^* ■* 
with fewer loops than in the original J^^\ and an i?^ with more loops than in 
the original Ra- Then we can repeat the steps described above with this new 
integral. 

Every subsequent term in the square bracket of Eq. ()2.44j) can be treated 
the same way as the first term. This allows us to express the integral in Eq. 
()2.42|1 in terms of k~ integrals over some J^" ^s, which have the same or fewer 
number of loops than the original Ja^\ 

/M ri'-l 
n K ~ (pa, q, V, vb; k[,..., k'^,) . (2.46) 

We now want to set /Cj = in order to put Eq. ()2.42|) into the form of Eq. 
()2.40|) . To that end, we employ an identity for J^" ^ (we again suppress the 
dependence on the color indices for brevity) 

, . . . , /c^, , k^j^^ . . . , 

(2.47) 



Ja iPA,q,v,VB;k[,...,k'^,) 



Ja'^ (PA,q,V,VB;k'^^ = 0, = 0,k'i 

n'-l + 

+5: 



i=l 



dl 



d j{n') I J I 

J A [PA,q,V,VB] 



k' + k' + ]+ k' + - a k' + - (]\ 
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Substituting the first term of Eq. ()2.47|) into Eq. ()2.46|) . we recognize the def- 
inition for Ta, Eq. ()2.4H) . We have shown in Sec. l2.4.Hl that the contributions 
from the terms proportional to dJ^ /dl'l in Eq. ()2.47|) can be expressed as 
soft-loop integrals of some J^" \ again with fewer loops than in J^^ \ When 
we substitute this into Eq. ()2.46|) we may express the resulting contribution 
in terms of integrals which have the form of Eq. ()2.42j) . We can now repeat all 
the steps mentioned so far, with this new integral. By this iterative procedure 
we can transfer the k~ integrals in Eq. ()2.42|) to J^""* and also set kf = inside 
J^^. In a similar manner, we can analyze the integrals in Eq. ()2.27j) . and 
express them in terms of defined in Eq. 1)2.411) . This algorithm, indeed, 
leads from the first factorized form of the considered amplitude, Eq. ()2.27|) . 
to the second factorized form, Eq. ()2.4()j) . 

2.5.2 Evolution equation 

We have now collected all the ingredients necessary to derive the evolution 
equations for quantities defined in Eq. ()2.41|1 . Consider F^''. We aim to find 
an expression for p\dT^^'' / dp^. As discussed in Sec. 12.4.21 this will enable 
us to resum the large logarithms of In(pjj) and eventually the logarithms of 
ln(s/|t|). According to Eq. ()2.41|) . in order to find p\dr^^^ /dp\, we need to 
study PaOJ^J"^ /dp\. Using the identities Pa ■ vb = '^C.aCb, Pa ■ h = 2^yi6, 
where = k^r]^ and ^a, Cb are defined in Eq. (|2.3ip . we conclude that 

= Q J(«) (e^, {^J-i, t, {q^ ■ k,^}--,\ {k,^ ■ k,^}lJ^,) . (2.48) 
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From this structure, using the chain rule, we derive the following relation 
satisfied by J^"^ which generalizes Eq. ()2.3Hj) to J4"'' with arbitrary number 
of external gluons. 



Pa- 



dJ' 



in) 



(n) 



n-l 



drj' 



dJ 



in) 



1=1 



dk- 



+ 



dJ 



in) 



(2.49) 



Now, we integrate both sides of Eq. ()2.49j) over 11^=1 "^^7) 

kj' = 0. Then, using the definition for Eq. ()2.41|) . the left hand side is 

nothing else but p^dT^^^ / dp\. The first term on the right hand side of Eq. 

is simply -fj'^dT^^'^ /dr]''. Noting that CbO/a^/OCb = n./^\ the last 
term gives simply n F^'' . For the middle term, we use integration by parts 



n-l / i-M \ n-\ rj j-(n) n-l , „M x n-l 

n(/.,*7)E*.-|^ = n(/„«7)E 

7=1 V^-Af / Cfci \J-M J 



dk; 



(n) 



J = l ' t=v - J-- 

n~l /n-l \ 

E J ^ iU^^Jj M [j^A\K = +M,...) + f^\k; = -M,...) 

-(n-l)Fj\ 



(2.50) 



Combining the partial results, Eqs. (j2.49|) and (j2.5(J|) . we obtain the following 
evolution equation 



Pa 



5fJ) 



dpj 



n-l pM /n-l 



i=l 



M 



M 



4HK = +M,...)+jriK = -M,... 



in)(U- 



in) 



, r(n) _ -g dVA 



(2.51) 
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The jet function J^' in the first term of Eq. ()2.51|1 is evaluated at {kf = 0}"^^ 
and the k~s are integrated over for j = 1, . . . , n — 1 and j ^ i- The first term 
in Eq. ()2.5H) can be analyzed using the K-G decomposition for gluon i since 
the k~ is evaluated at the scale M ~ \/\t\. The outcome of the last term 
in Eq. ()2.5ip has been determined in Sec. 12.4.21 Fig. 12.41 ^. As a result we 
have all the tools necessary to determine the asymptotic behavior of the high 
energy amplitude for process ()2.1|) . To demonstrate this, we will rewrite Eq. 
(I2.51|) into the form where on the right hand side there will be a sum of terms 
involving \ convoluted with functions which do not depend on p\. Let us 
proceed term by term. 

Again, the K-G decomposition applies to the first term in Eq. 1)2.511) be- 
cause the external momenta are fixed with k~ = ±M. Using the factorization 
of a K gluon given in Eq. ()2.30|) it is clear that the contributions from the K 
gluons cancel for J^^s evaluated at k~ = +M and k^ = — M. Hence only the 
G gluon contribution survives in this term. Its most general form is shown in 
Fig. 12.3b . Before writing it down let us introduce the following notation. For 
a set of indices {1,2, . . . ,n}\{i} consider all the possible subsets of this set, 
with 1, 2, . . . , (n — 1) number of elements. Let us denote a given subset by vr, 
its complementary subset vf, the number of elements in this subset as and 

in its complementary as = (n — 1) — n^. With this notation, we can write 

^Strictly speaking we have analyzed rj'^ d j'^^^ / d t]" , but because of the relationship be- 
tween J^"-* and r^-* given by Eq. (|2.41|l . once we know t)"' d j')^^ / d rj" we also know 
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the ith contribution to the first term in Eq. ()2.51|) in the form 



J' 



(n) ai ... an 



A 



(K = +M, . . .) + j^^' {k- = -M, ...)=$: y n ^2.)- 



'^aiai.^...ai^^ bi...fejv 



(k- = +M, k-, . . . , kl^;kt = 0, A;+ = 0, . . . , 
ki±, ±, • • • , ki^^ ±; /i, • • • , ^Af; Q', ^7) 



0; 



X J 



{nfi+N)aj ...aj 6i...b]v 

1 ""TT 



hj- ..... /c- , k- — .... , k- — , 



kii±,---,kin^±;lu---jN;PA,q,v) + (^i ^ -^)- 



(2.52) 



In Eq. (j2.52p . the summation over repeated indices is understood. We sum 
over all possible subsets vr. In other words, we sum over all possible attach- 
ments of external gluons to jet function Ja and to the soft function S. The 
elements of a given set tt are denoted 11,12, ... , in^,. The elements of a comple- 
mentary set 7f are labeled 11,12, .. . .,in^. The number of gluons connecting S 
and Ji"^+^^ is N. 

Following the procedure described in Sec. 12.5.11 with Ra in Eq. ()2.42j) 
replaced by S in Eq. ()2.52|1 . we can express the contribution from a G gluon 
in the first term of Eq. ()2.51|) in a form 

n-1 I'M /n-l \ 

j(n)a,...an^j^- ^ _ _ ^ ^(n) a, ... a„ = -M, . . .) 

i=l " - \ / 
« m 

m j=l 

X T^r^''-''"{pA,q,vJi±,...Jm±;M). (2.53) 
The function /C^"^'™) does not contain any dependence on pA. It can contain 
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delta functions setting some of the color indices bt, as well as transverse mo- 
menta li± of equal to color indices Oj and transverse momenta ki± of 

p(n) 
J- A ■ 

Next we turn our attention to the last term appearing in Eq. (j2.5ip . The 
contribution to this term has been depicted graphically in Fig. 12 .41 Consider 
the term in Fig. 12.4b . It can be written in a form 

^ X] ^'a^ a,, ... a,„^ b ) ? • • • ? ! K = 0, /c^^ = 0, . . . , = 0; 

ki±,ki-^±, . . . , ki^^ _|_ ; / = fcj + /cj^ + . . . + ki^^ 

X ^ ""^ {k-.^, . . . , k-.^_ ; fc^^ = 0, . . . , = 0; /Ci^ _L, . . . , _l; 

/ = /c, + A;^^ + . . . + fci^^ ; Pa, q, v)- (2-54) 

In Eq. fl2.54|l . we have used the same notation as in Eq. ()2.52|1 . Momentum / 
connects S' with J^''^^\ Following the same procedure as in Sec. 12.5.11 with 
Ra appearing in Eq. ()2.42j) replaced by S' introduced in Eq. ()2.54|1 . we can 
express this contribution in a form given by Eq. ()2.53|) with a different kernel 
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The contribution from Fig. 12.4b can be written 




0; 




- / _ _ _|_ _|_ 

\ n' ' n ' ' ins- ' 



(2.55) 



The flow of momenta k and I is exhibited in Fig. 12. 4b . The momentum k flows 
through the boxed vertex and the ghost line shown in Fig. 12.4b which forces 
this momentum to be soft, so that lines k and / are part of the function S". 
Since the line with momentum k is soft, then all gluons attaching to J^"'"^^'' 
in Eq. ()2.55|) are soft and we can again apply the procedure described in Sec. 
12.5. II to bring the contribution in Fig. 12.4b into the form given by Eq. ()2.5Hj) 
with a different kernel, of course. 

In summary, we have demonstrated that all the terms on the right hand 
side of Eq. 1)2.511) can be put into the form given by Eq. ()2.53j) . This indicates 
that Eq. ()2.5H) . indeed, describes the evolution of F^'' in InpJ^ since it can be 
written as 




m j=l 



. , -n(™) bi...bm / r I \ 



(2.56) 
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The kernels /C^"'"^) do not depend on p\. As indicated above, they can contain 
delta functions setting some of the color indices hi, as well as transverse mo- 
menta li± of r^'* equal to color indices and transverse momenta of T^'*. 
The systematic use of this evolution equation enables us to resum large loga- 
rithms In(pJ^) at arbitrary level of logarithmic accuracy. Analogous equation 
is satisfied by r^. It resums logarithms of ln(p^). 

2.5.3 Counting the number of logarithms 

Having derived the evolution equations for F^'', Eqs. ()2.51|) and ()2.56p . it does 
not take too much effort to show that at r-loop order the amplitude contains at 
most r powers of ln(s/|t|). We follow the method of Ref. [16]. We have argued 
in Sec. 12.5.11 that the power of ln(s/|t|) in the overall amplitude corresponds 
to the power of In(pJ^) in F^"*. So we have to demonstrate that at r-loop 
order F^''^'', where F^'^'' represents a contribution to F^'* at r-loop level, does 
not contain more than r logarithms of In(pJ^). We prove this statement by 
induction. First of all, the tree level contribution to F^''^'* is proportional to 
the expression 

where t°''is are the generators of the SU{3) algebra in the fundamental rep- 
resentation. The sum over {ii, . . . , in} indicates that we sum over all possi- 
ble insertions of the external soft gluons. Eq. ()2.57p is evaluated at {k^' = 
0}^^i. Expanding the denominators in Eq. ()2.57j) we obtain the expression 
—2p\{k^_^ -|- . . . + k^,) — [ki^ + . . . -|- ki.)\ + ie. We see that the poles in k~ 
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planes are not pinched and therefore the k~ integrals cannot produce In(p^) 
enhancements. 

Next we assume that the statement is true at r-loop order, and show that 
it then also holds at (r + l)-loop level. To this end we consider the evolution 
equation, Eq. ()2.5H) . and examine {p\d/dp\ — 1) F^^^''^^"'^'*. Its contribution 
is given by the first and the third term on the right hand side of Eq. fl2.51|) . 
As already mentioned, the first term in Eq. (j2.51|) can be analyzed using 
K-G decomposition. The contributions from the K terms cancel each other 
while the contribution from the G gluons are given by the kind of diagram 
shown in Fig. 12.3b . The latter, however, can be written as a sum of soft loop 

in' r') 

integrals over J)^ ' with r < r, since we loose at least one loop in the original 
j{n,r+i) g^^^ momentum integration. This is demonstrated in Eq. 

fj2.52|) . Following the procedure described in Sec. 12.5.11 we may express these 
contributions as transverse momentum integrals of some F^ \ see Eq. fl2.53|) . 
These contain at most r' < r logarithms of In(pJ^). The contribution from the 
third term in the evolution equation, Eq. ()2.51|) . is given by the diagrams 
depicted in Fig. 12.41 These are again soft loop integrals of some ''^'^ with 
r' < r, and they can be expressed as transverse momentum integrals of F^ \ 
see Eqs. fl2.54|) and ()2.55p . which have, therefore, at most r logarithms of 
In(pJ^). Since both terms on the right hand side of Eq. (|2.51|) have at most r 
logarithms of ln{p\), then also p^dT^^'^^^^ /'^Pa most r logarithms of 

In(p^) at (r + l)-loop level. This immediately shows that r^'^+^) itself cannot 
have more than (r + 1) logarithms of In(p^) at (r + l)-loop level. 

This result enables us to formally classify the types of diagrams which con- 
tribute to the amplitude at the A;-th nonleading logarithm level. As has been 
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shown in Sec. 12.5.11 we can write an arbitrary contribution to the amphtude 
for process ()2.H1 in the Regge hmit in the second factorized form given by 
Eq. ()2.40|) . Consider an r-loop contribution to the amphtude and let L^, Lb 
and Ls be the number of loops contained in Ta, Tb and S. Since Ta (Tb) 
can contain La (Lb) number of logarithms of p\ (p^) at most, the maximum 
number of logarithms, A/'maxLog, we can get is 

^maxLog = — Ls- (2.58) 

This indicates that when evaluating the amplitude at the /c-th nonleading 
approximation, we need to consider diagrams where 1,2,..., {k+1) soft gluons 
are exchanged between the jet functions J a and Jb- 

2.5.4 Solution of the evolution equations 

Having obtained the evolution equations, Eqs. ()2.51|1 and ()2.5(i|l . we discuss 
how to construct their solution. Our starting point is Eq. ()2.56p . In shorthand 
notation it reads 

p\ r?-) = Y.Y1 ® (2.59) 

op A ^,=0 n' 

at r-loop level. Indices n and ra', besides denoting the number of external 
gluons of the jet function, also label the transverse momenta and the color 
indices of these gluons. The symbol ® in Eq. ()2.59|) denotes convolution over 
the transverse momenta and the color indices. Note that Eq. (j2.59|) holds for 
V A with the overall factor p\ divided out {Va = ^a/Pa)- We have proved. 
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in Sec. 12.5.31 that T^'^ can contain at most r logarithms of In(pJ^) at r-loop 
level. Therefore the most general expansion for is 

rr^-Ec^^ln^(pl)- (2.60) 

j=0 

If we want to know r^'^-* at N'^LL accuracy (fc = is LL, = 1 is NLL, etc.), 
we need to find all c^p^^ such that r — j < k. The coefficients c^"''"'' in Eq. 
(I2.60|l depend on the transverse momenta and the color indices of the external 
gluons. Using the expansion for and Eq. ()2.60|) . in Eq. ()2.59|) 

and comparing the coefficients with the same power of In(pJ^), we obtain the 
recursive relation satisfied by the coefficients c^''^^ 

r—l n+r—r' 

jc'f''^= ^ Yl '''"^'^ ® cf_Y'^ (2.61) 

r'=j — l n'=l 

In Eq. (12.611) . we have used that, in general, l<n'<n + r — r\ 

We now show that Eq. ()2.61|) enables us to determine all the relevant 
coefficients Cj'''^^ of order by order in perturbation theory at arbitrary 
logarithmic accuracy. We start at LL, k = 0, and consider n = 1. At r-loop 
level we need to find the coefficient ci^'^\ It can be expressed in terms of lower 
loop coefficients using Eq. ()2.6ip and setting j = r and n = 1 

rc^^) = ^^(i."';i)®4"_'r'^- (2.62) 

n'=l 

In Sec. 13.11 we will prove that the one loop kernel satisfies /C-^'^'^-' = 0, Eq. 
fl3.6|) . This implies that in Eq. ()2.62j) the coefficient ci^'^^ is expressed in terms 
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of lower loop coefficient c^J.i~ and hence, we can construct the coefficients 
at arbitrary loop level once we compute Cq^'^'' , the coefficient corresponding to 
the tree level jet function r^^'^\ 

Next we construct all for n > 1 at LL accuracy. Let us assume that 
we know all cl" '^'^ for all r and for n' < n. We apply Eq. ()2.6H) for j = r 

r c^"-'-) = J2 ''^ ® ci!.'r'\ (2.63) 

n'=l 

In Sec. I3.2l we will show that the evolution kernel in Eq. ()2.(i3j) obeys /C^"'"''^) = 
e{n-n') £("''^''1), Eq. dSSHl), where e{n-n') is the step function. This implies 
that the sum over n' in Eq. ()2.63|) terminates at n' = n. Isolating this term 
in Eq. ()2.(i3|l . we can write 

n-l 

r 4^^'^) = ^("•'^^i) ® ctr'^ + Yl ''^ ® ^i-'r'^- (2-64) 

n'=l 

So after we calculate the tree level coefficient c'q^'^\ we can construct all the 
coefficients Cr"'*^^ using Eq. ()2.64|) order by order in perturbation theory, since 
according to the assumption we know Cr" '^^ for all r and for n' < n. This 
proves that we can construct the jet functions at LL, = 0, for all n to all 
loops. 

We now assume that we have constructed all the jet functions at the N^LL 
accuracy for a given A; > and we will show that we can determine all the jet 
functions at the N'^+^LL level. We start with n = 1. Using Eq. ()2.fil|) with 
n = 1, j = r — (/c + 1), isolating the term with r' = r — 1 in the sum over r' 
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and using /C(i'"''^) = ^^^z A^(^'^;^), we arrive at 

(r-fc-i)e2-i = ^^^'^^^^®ci^:ri+ E ^(^'"'-^')®c;:ri (2.65) 

r'=r~k—2 n'=l 

After we evaluate the coefficient Cg^''^^^^ (impact factor), Eq. ()2.65|) implies that 
we can calculate the coefficients c^r-k-i oi'der by order in perturbation theory, 
because, according to the induction assumption, we know all the coefficients 
since they are at most N'^LL. Once the coefficients of F^"* are determined 
at N'^"'"^LL level, we assume that we know all the coefficients of T^^ ^s for n' < n. 
We want to show that we can now construct all the coefficients for at 
N*'"'"-^LL accuracy. First we need to calculate c^''''^^\ Then we use Eq. ()2.6H) 
to express the coefficient c,[,"'^^_j^, isolating the terms with r' = r — 1 and n' = n, 
as 

n-l 
n'=l 

2 n+1 — r' 

+ E E ® ^^-r (2-66) 

r'=r—k—2 n'=l 

The terms appearing in the sum over r' in Eq. ()2.66|) are known according to 
the assumptions since for them r' — {r — k — 2) < k. We also know, according to 
the induction assumptions, the contributions to the second term of Eq. ()2.66|) . 
since they have n' < n. Therefore, we can construct c^"'^Li order by order in 
perturbation theory. This finishes our proof that we can determine the high 
energy behavior of F^^ at arbitrary logarithmic accuracy. Note that to any 
fixed accuracy only a finite number of fixed-order calculations of kernels and 
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coefficients Cq ' must be carried out. In a similar way we can construct a 
solution for 

Once we know the high energy behavior for T^^^ and T^^\ then the second 
factorized form, Eq. ()2.40|) . implies that we also know the high energy behavior 
for the overall amplitude. Because a jet function F*^"^ is always associated with 
at least n — 1 soft loop momentum integrals in the amplitude, we infer from 
Eq. ()2.58p that if we want to know this amplitude at N^LL accuracy, it is 
sufficient to know T^^'^ (rg"^) at N^+i-"LL (N^+^-^LL) level forn < if + 1 
{m < K + 1). We note, however, that to construct these functions according 
to the algorithm above, it may be necessary to go to slightly larger, although 
always finite, values of n and m. Let us describe how this comes about, starting 
with the basic recursion relations for coefficients, Eq. ()2.6H) . 

We assume that for fixed n on the left-hand side of Eq. ()2.61|) . the loga- 
rithmic accuracy k is bounded by the value necessary to determine the overall 
amplitude to Kth nonleading logarithm: k = r — j < K + 1 — n, which we 
may rewrite as n + r — {K + 1) < j < r. On the right-hand side of Eq. 
()2.61|) we encounter the coefficients of the jet functions with n' external lines, 
satisfying the inequality n' < n + r — r' < n + r — {j — 1). Combining these 
two inequalities, we immediately obtain that n' < K + 2. Then, for any given 
number of external gluons n' on the right-hand side, we encounter a level of 
logarithmic accuracy k' = r' — (j — 1) < n + r — n' — {j — 1) < K -\-2 — n'. This 
reasoning indicates that, in general, we will need all ■* (F^ ^) at N^"*'^~"'LL 
level for n' < K + 2 (m' < K + 2), when evaluating the ampli- 
tude at N-^LL accuracy. We note that for fermion exchange in QED it was 
shown in Ref. [16] that only contributions with n' < K + 1 are nonzero. As we 
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will witness in Chapter El two-loop calculations appear to indicate, Ref. [11], 
that QCD requires the full range of n' identified above, starting at NLL. 

2.6 Conclusions 

We have established a systematic method that shows that it is possible to 
resum the large logarithms appearing in the perturbation series of scattering 
amplitudes for 2^2 partonic processes to arbitrary logarithmic accuracy in 
the Regge limit. Up to corrections suppressed by powers of \t\/ s, the amplitude 
can be expressed as a sum of convolutions in transverse momentum space over 
soft and jet functions, Eq. ()2.4()|1 . All the large logarithms are organized in 
the jet functions, Eq. ()2.4H1 . They are resummed using Eqs. ()2.5H1 and/or 
()2.56|) . The evolution kernel K, in Eq. ()2.56|) is a calculable function of its 
arguments order by order in perturbation theory. This is the central result of 
our analysis. 

The derivation of the evolution equations and the procedure for finding 
the kernels were given above in Coulomb gauge. Clearly, it will be useful and 
interesting to reformulate our arguments in covariant gauges. In addition, the 
connection of our formalism to the effective action approach to small-x and 
the Regge limit, Refs. [35,36] should provide further insight. 
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Chapter 3 

The Applications 



In the previous chapter, we have developed the general formalism for obtain- 
ing the high-energy behavior of the scattering amplitude for process ()2.1|) at 
arbitrary logarithmic accuracy. In this chapter, we apply these techniques to 
study this amplitude at LL and NLL level. 

a) b) 

Figure 3.1: Diagrams contributing to the amplitude at NLL approximation: 
factorized one gluon exchange diagram (a) and non-factorized two gluon ex- 
change diagram (b). 
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3.1 Amplitude at LL 

According to Eq. ()2.58|) . the amplitude at LL comes solely from the factorized 
diagram shown in Fig. 13.1b .. but without any gluon self-energy corrections. 
The jet A, containing lines moving in the plus direction, and jet B, consisting 
of lines moving in the minus direction, interact via the exchange of a single soft 
gluon. This gluon couples to jet A via the — component of its polarization and 
to jet B via the + component of its polarization. Since v\Nap{q,ifjv^^ = 1, 
we can write at LL 



Asbs = -jJA^''{PA,q,v) j'k^''{PB,q,v), (3-1' 



where fog is the color basis vector corresponding to the octet exchange, defined 
in Eq. ()2.3p . Using s = 2p\p]^, the logarithmic derivative of the amplitude 
can be expressed as 

dAs _ laj^ 1 (i)adJ^ 



dins t d\np\ t d\np^ 

In Sec. 12.4.21 Eq. ()2.33p . we have derived an evolution equation resumming 
ln(]9^) in We note that J^^ = T^^\ and that (j2.33j) is a special case of the 
evolution equation ()2.51|) . The diagrammatic representation of the first term 
on the far right hand side of Eq. ()2.33p . which follows from Fig. 12.41 in the 
case when we have one external soft gluon attached to a jet function, is given 
by the diagrams in Fig. 13.21 Diagram in Fig. 13.2b corresponds to Fig. 12.4b 
and the diagrams in Figs. 13.2b and 13.2b correspond to Fig. 12.4b for n = 1. 
The diagrams in Figs. 13.2b and 13.2b are in the factorized form, while the one 
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a) b) c) 

Figure 3.2: Diagrammatic representation of the evolution equation for jet J^^ 
at LL. 

in Fig. 13.2b is not. 

As discussed in Sec. I2.4.2| power counting shows that the loop momentum 
k in Fig. 13.2b must be soft. This implies that we can make the following 
approximations. First, since at LL all internal lines of the jet A are coUinear 
to the + direction, we can neglect the dependence of J\ , i.e. we may 
set = inside J^^\ Also, we can pick the plus components of the vertices 

(2) 

where the soft gluons attach to the jet J\ . A short calculation, which uses 
the Feynman rules for special lines and vertices listed in Appendix I A. 41 gives 
the contribution to Fig. 13.2b in a form 

Fig.E2K = 

(3.3) 

where we have defined = gsfi"^- Using Eqs. 1)2.111) and ()2.36p for the 
components of the gluon propagator and the boxed vertex, respectively, it 
is easy to see that in the Coulomb (Glauber) region, fc" ^ A;^ ~ the 
integrand in Eq. ()3.3p becomes an antisymmetric function of k^ and that 
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a) b) c) 

Figure 3.3: Diagrams determining the contributions to the gluon trajectory 
at the order a^. 

therefore the integration over vanishes in this region. 

In the soft region, where all the components of soft momenta are of the 
same size we can use the K-G decomposition for the soft gluon with 

(2) 

momentum k attached to J\ . At LL, however, there cannot be any soft 

('2) _____ I— 

internal lines in J)^ in Eq. ()3.3|) . since, as discussed in Sec. 12.5.31 only 
integrals over coUinear momenta can produce powers of lnpj[. Therefore, at 
LL, only the K gluon contributes, because the G gluon must be attached to 
a soft line. The K gluon can be decoupled from the rest of the jet using 
the Ward identities, Eq. ()2.3()|) . Their application in Eq. ()3.3|1 gives 

X v%N^,{q-k)v''^/^^\pA,q,v)- (3-4) 

We have used the identity facb fdcb = N^dad = CA^ad in Eq. Eq. 
now gives a factorized form for Fig. 13.2b . Since the contributions in Figs. 13.2b 
and 13.2b are already in the factorized form, we can immediately infer that the 
gluon reggeizes at LL. Combining the terms from Fig. 13.21 in Eq. ()2.33j) . we 
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obtain the evolution equation at leading logarithm 



d 

p\-Q-fJA^''{pA,q,v) = a{t) JA^''{pA,q,v)- (3-5) 

Using the notation for evolution kernels introduced in Sec. I2.5.4[ Eq. ()3.5p 
implies that 

/C(^'2'^) = 0. (3.6) 

InEq. (HSl) 

+ + + (3.7) 

is the gluon trajectory up to the order ag, and a^a\t), a^\t) and a^c\t) are 
its contributions given in Figs. 13.3b - 13.3b . respectively, 



X Vf,/3uik,-q,q - k)vs, 



/3 

d^k 1 



(3.8) 



In Eq. (13. 8|) . V^py{k, —q, q — k) stands for the momentum part of the three- 
point gluon vertex. After contracting the tensor structures in Eq. ()3.8p . using 
the explicit form for V^/ji., va, vb, 5*^ (Eq. ()2.34j) ) and for the components of 
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the gluon propagator, Eq. (|2.1ip . we obtain for a^abci^)^ 



" ^ ' ^2 J (27r)^ {ko + ks) A;2 {k - qf {k ■ If [k - q) ■ {k - q) 



X 



1 r d^k 



{k - qf + 2P kl (k - q)^ ■ q^ + 2kl kl k^ ■ (k - q)^ + 2kl kl {k - q)\] 



^ ^«;c4 / |^_^_iL__ . (3,9) 



Next, we perform the and k^ integrals in Eq. ()3.9p . For ai^^()f:), these 
integrals are UV/IR finite. However in the case of a^^],(t), the integral is 
linearly UV divergent. In order to regularize this energy integral, we invoke 
split dimensional regularization introduced in Ref. [49]. The idea is to regu- 
larize separately the energy and the spatial momentum integrals, i.e. to write 
d'^kE d^^kid^^k for Euclidean loop momenta kE- The dimensions Di and 
D2 are given by Di = l—2ei and D2 = 3—2e2, with ej 0+ for j = 1, 2. Since 
the energy integral for a^c^ (t) is scaleless, it vanishes in this split dimensional 
regularization. The energy integrals in a^ali^) straightforward. 

All the k^ integrals can be expressed as derivatives with respect to k'^ 
and/or (A; — q)^ of a single integral 



/(„, 6) . r . ' = 1./'' ^ - 



^/kl + a2 {kl + IP) bVW^ 



(3.10) 
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The result of these integrations over is 



D-2 



A;2 1 



X { I{\k±\Ak±- q±\) 



[(k-q)l-kl]^ [(A:-g)i-A:i]V' 
= 0. (3.11) 

Combining the results of Eq. ()3.11|) and Eq. ()3.7p . we obtain the standard 
expression for the gluon trajectory at LL 

a{t) = 1 + C^a.^.^^ J ■ (3.12) 

We can now simply solve the evolution equation ()3.2j) . to derive the factorized 
(reggeized) form for the amplitude in the color octet 

A8(s,t,a,) = s"Wi8(t,a,). (3.13) 

The amplitude factorizes into the universal factor s"^*-*, which is common for 
all processes involving two partons in the initial and final state and dominated 
by the gluon exchange, and the part Ag, the so-called impact factor, which is 
specific to the process under consideration. 
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3.2 Amplitude at NLL 

At NLL level the contribution to the amplitude comes from both the one gluon 
exchange diagram, Fig. 13.1b . and from the two gluon exchange diagram, Fig. 
13.1b . At this level, both singlet and octet color exchange are possible in the 
latter. Including the self-energy corrections to the propagator of the exchanged 
gluon (taking into account the corresponding counter-terms), we can write the 
contribution from the diagram in Fig. 13.1b as follows, 

^ -i (p^, g, r/) (^N-^{q, v) + jv'^^ < n,.(g, r/) t;^ /^^^ {ps, q, v), 

(3.14) 

where Il^i,{q,ri) stands for the one loop gluon self-energy. We now put this 
contribution into the first factorized form, Eq. ()2.27j) . isolating the plus po- 
larization for jet A, and the minus polarization for jet B. At NLL in the 
amplitude, we need the soft function S**^^'^-*, Eq. ()2.27|) with n = m = 1, 
to accuracy 0{as)- Using the tulip-garden formalism described in Appendix 
IA.3| the contribution to the first term on the right hand side of Eq. ()3.14|) is 
given by the subtractions shown in Fig. 13.41 In accordance with the notation 
adopted in Appendix IA.3| the dashed lines indicate that we have made soft 
approximations on gluons that are cut by them. A dashed line cutting a gluon 
attached to jet A{B) means that the gluon is attached to the corresponding jet 
through minus(plus) component of its polarization. Since = in the Regge 
limit, Eq. ()2.2|) . we have N'^^{q) = g^^. This implies that the contributions 
between the diagrams in Fig. 13.4b and in Fig. I3.4l i as well as between the 
diagrams in Fig. 13.4b and in Fig. 13. 4f cancel each other. Therefore only the 
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+ 



b) 



c) 





Figure 3.4: Expansion of the one gluon exchange amphtude at NLL using the 
tuhp garden formahsm. 

zeroth-order soft function diagram in Fig. 13.4b survives in the factorized form, 
Eq. (ITTTI) . 

For the two gluon exchange, Fig. 13.1b . we only need the lowest order soft 
function at NLL in the amplitude (and LL in singlet exchange). The expression 
for the two gluon exchange diagram in Fig. 13.1b takes the form, Eq. ()2.27j) . 



j{2) ab^f^+ = 0,k-, k^) S{k+, k-, k^) J^^> ""(A;- = 0, k+, k 



where S{k) is given by 



(2) a 6/, - 



(3.15) 



Q(k) = 1 ^"^(^) N- + {q-k) 
^ ' 2\ k^ + te {q-kf + ie' 



(3.16) 



We have suppressed the dependence of the functions appearing in Eq. ()3.44p 
on other arguments for brevity. At NLL accuracy we are entitled to pick the 
plus Lorentz indices for jet function J a and the minus indices for jet function 
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Jb only. We can also set /c+ = in Ja and fc" = in Jb since all the loop 
momenta inside the jets are collinear. Eq. ()3.44j) represents the first factorized 
form, Eq. for the amplitude A^'^\ 

Next, we follow the procedure described in Sec. 12.5. ll to bring the amplitude 
into the second factorized form, Eq. ()2.40|) . We employ an identity based on 
Eq. for the function S{k) defined in Eq. 

s{k+,k-) = s{k+ = o,k- = o)e{M -\k+\)e{M -\k-\) 

+ [s{k+, k~ = 0)- s{k+ = o,k- = 0) e{M - \k+\)] e{M - \k-\) 
+ [s{k+ = 0, k~) - s{k+ = o,k~ = 0) e{M - \k-\)] e{M - \k+\) 

+ [{S{k+, k-) - S{k+, k- = 0) e{M - - 

{s{k+ = 0, k-) - s{k+ = o,k- = 0) e{M - \k-\)} e{M - \k+\)]. 

(3.17) 

The contribution from the first term in Eq. ()3.17j) gives immediately the 
second factorized form with F^^^ and defined in Eq. ()2.4H) for n = m = 2. 

We now discuss the rest of the terms in Eq. ()3.17j) . which can be analyzed 
using the K-G decomposition, since, by construction, there is no contribution 
from the Glauber region. At the current accuracy only the i^-gluon con- 
tributes. After substituting the second term of Eq. ()3.17j) into Eq. ()3.44|1 . 

(2) 

we can factor the gluon with momentum k from jet . However, it is easy 
to verify, using the definitions for K and G gluons, Eq. ()2.29j) . the Ward 
identities, Eq. ()2.30|1 . and the explicit components of the gluon propagator, 
Eq. ()2.11|) . that the k~^ integral is over an antisymmetric function. As a re- 
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suit, this contribution vanishes. In a similar fashion, the contribution from 
the third term in Eq. ()H.17|1 . after used in Eq. ()H.44|1 . vanishes, since now we 

(2) 

can factor the soft gluon with momentum k from jet J\ and the k~ integral 
is over an antisymmetric function. 

In the case of the last term in Eq. ()3.17|) . after used in Eq. ()3.44|) . we 
can factor the soft gluon with momentum k from both jets J)^ and . The 
integrals of the soft function S{k) over k~^ and k~ are then 

^2 I'M 17+ 17- 

S{k^,q;M)^CAj§^ J ^^^S{k\k-,k^,q). (3.18) 

As usually, we leave the transverse momentum integral undone. The 1/k^ 
and l/k~ in the integral above are given by the Principal Value prescription 
because there is no contribution from the Glauber region. Since the amplitude 
is independent on the choice of scale M, we can evaluate it at arbitrary scale. 
We choose to work in the limit M — 0. In this limit the contribution to 
the integral comes from the imaginary parts of the gluon propagators in Eq. 
()3.16|) . —iTT6{k'^) and —m5{{k — qY). The integration is then trivial and Eq. 
()3.18|) becomes 

S(k,, ,) . lim S(k,, M) . -C. (3.19) 

Combining the partial results of the analysis described above in Eq. ()3.44p . we 
arrive at the second factorized form for the double gluon exchange amplitude. 
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Fig. mb, 



+ 



/ 
/ 



(2^ 




(2) ab 



ik±)j7^S{k+ = 0,k- = 0,k^)T 



(3.20) 



Using Eq. ()3.14j) for A^^^ and Eq. ()3.2()|1 for A^'^'' , we obtain the amplitude for 
the process ()2.ip at NLL accuracy 



In Eq. (13.2111 ■ we have used the exphcit form for S{k^ = 0,k^ = 0,k±), 
which can be easily identified from Eq. ()3.1(i|) . We have also used the integral 
representation of the gluon trajectory given in Eq. ()3.12|) . 

In order to determine the high energy behavior of the amplitude in Eq. 
()3.21|) . we need to examine the high energy behavior of F^'* or r^"* at NLL 
and the evolution of or F^-* at LL. In this section, we restrict the discussion 

(2) 

of evolution equations to LL level, and hence we analyze the behavior of T\ 
only. We will address the study of NLL jet evolution, and gluon reggeization 
at this level in the next section. 

We use the evolution equation given by Eq. ()2.51|1 in order to determine 
the LL dependence of F^-* on \n{p\). In our special case of the two gluon 




{PB,q,v) 



(3.21) 
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a) 



b) 






(2) 

Figure 3.5: Diagrams determining the evolution of TV. 



exchange amphtude, it reads 



d 



Pa 



dp] 



1 r 



(2) ab 



M 



jf"\k- = +M,k+ = 0,A;j 



+ j'x^^'^k- = -M,k+ = 0,kj 



d 



p(2) ab 
^ A 



(3.22) 



The first term in Eq. ()3.22j) can be analyzed using the K-G decomposition. 

(2) 

The contributions from the i^'-gluon cancel between the J\ {k~ = +M) and 
J^\k^ = — M). The contributions from the G gluon, which we now discuss, 
are shown in Figs. 13.5b and 13.5b . 

Since the gluon with momentum q — k in Fig. 13.5b cannot be in the Glauber 



region, we can use K-G decomposition on it. The K part factors from J 



(3) 
A ; 



while the G part does not contribute at LL. After factoring out the gluon with 



momentum q — k and performing the approximations on the jet function J 



(2) 
A ' 
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the contribution to Fig. 13.5b for k = 



+M is 



Fig. = -tglfaecfdeb ^ J Si{k+ = 0,k~ = +M, k^, I) 

X jf ^V = 0,r,/^), (3.23) 



where we have defined 

Next we follow the established procedure. First, we write 

Si{k,i) = Si{k,i- = o)e{M-\i-\)+ [Si{k,i) - Si{k,i- = o)e{M-\i-\)] . 

(3.25) 

When we use the second term of Eq. (I3.25|l in Eq. ()3.23|1 . we can factor 

(2) 

the gluon with momentum I from J\ . Since the resulting integrand is an 
antisymmetric function under the simultaneous transformation M —M, 
—l"^, the contributions on the right hand side of Eq. ()3.22p evaluated 
for k^ = +M and k~ = —M cancel each other. Therefore we can write, using 
Eq. (jT^ in Eq. (jS2Sl), 

Fig. IS = -^g!faecU,j^ I / d/+ 

X Si {k+ = o,k~ = +M, k^j- = 0, /+, h) rj^ + . . . , 

(3.26) 

where by dots we mean the term which is canceled after we take into account 
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(2) (2) 

the contributions to both J\'{k~ = +M) and J\ {k^ = —M) on the right 
hand side of Eq. (IT^ . 

Next, we perform the integral in Eq. ()H.2fi|) . As we have aheady men- 
tioned above, since the final result does not depend on the scale M, we can 
choose arbitrary value of M. We have chosen to perform the calculation in 
the limit M ^ 0. Then the only non- vanishing contribution comes from the 
imaginary part of the propagator — k)"^ +«e], — "ivr 5{2Ml~^ + k)\). For 
this term the integration is trivial and we obtain 

M{Y\g.^) = -asfaeche, j ^^2 \l 1-) + ■ ■ ■ ^ (3-27) 

which gives an M-independent contribution to the right hand side of Eq. 

(IS2S1). 

We follow the same steps when dealing with the diagram in Fig. 13.5b . 
whose soft sub-diagram is given by 

X V,s^{q-l,l-k,k-q). (3.28) 

First we use the identity ()3.25p for 5*2. The contribution due to the second term 
in Eq. (|3.25|) vanishes, after the gluon with momentum / has been factored 
from Jj^ , due to the antisymmetry of the integrand. Hence again, as in the 
case discussed above, only the term given by S2{1~ = 0,Z+,/_L,fc) contributes. 
In the limit M — > 0, the contribution comes from the imaginary part of the 
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same denominator as in the case of Fig. 13.5b . The result is 



M(Fig.E3b)= - asfaecfdeb J 



(2vr)^-2 II (/ - q)l (k - 1)1 
X (A^i/i - k± ■ l±l]_ - k± ■ q±l]_ - k]_l± ■ q± + 2k± ■ l±l± ■ q±) 

X T^2^'\h) + ... . (3.29) 



Combining the results of Eqs. ()3.27p and ()3.29|) . we obtain the expression for 
the surface term in Eq. ()3.22|1 

M [ jf ^\k- = +M, k+ = 0, k^) + jf ^\k- = -M, k+ = 0, k±) ] 



2Q^s faecfbed 



{2n)D-^ \ll {k - 1)1 (/ - q)l {k - 1)1 II (g - 1)1 
^T^a"\Il)- (3.30) 

Next, we analyze the contributions to the term 7]°" d / d 7]°" T^"^ in the evo- 
lution equation ()3.22|1 . The contributing diagrams are shown in Figs. 13. 5h - 
13. 5t '. Note that for every diagram in Figs. 13.5b - 13. 5t '. we have also diagrams 
when a loop containing the boxed vertex is attached to the external gluon with 
momentum fc, instead of to the external gluon with momentum q — k. 

In Fig. 13.5b . we have to consider all the possible insertions of external 
gluons with momenta k and q — k. We have six possibilities. The contribution 
shown in Fig. 13.5b is proportional to (omitting the color factor) 

A1.-(V ™ r Al.- f ^""^ N~^^{l) S,{l) jl-k)^ jq-iy 
dk (Fig.linfc)oc/ dk j J^^^^J^J—f^^k =0). 

(3.31) 



M J~M 



Since the integrand is an antisymmetric function under k — > —k and 
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a) 



b) 



<zi:> 



l-k< 



□- 



d) 



q — k 




+ 




q — k 



Figure 3.6: Contributions to the diagram in Fig. 13. 5f when the gluon coming 
out of the boxed vertex is attached to the soft hne (a) and when either or both 
gluons with momenta k and / are K gluons and they are factored from the jet 
(b-f). 

— the integral in Eq. ()3.3ip vanishes. The same antisymmetry property 
holds for the remaining five diagrams and therefore, there is no contribution 
from them. 

Let us next focus on the diagram in Fig. 13. 5f . When the gluon with 

(3) 

momentum I attaches to a soft line inside of the jet J\ , the contribution 
takes the form shown in Fig. 13.6b . If it attaches to a jet line, its contribution 
can be written as 



Fig.rm 



-9sfi 



bed 



d^/ 



(27r) 



D 



X J 



(3) acd 



{k' 



0,k-,ki,l+ = o,r,/i 



(3.32) 
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with the soft function 

Ssik, l) = {q- k) _ (g-A;-/)^ • 

We use the identity for this soft function S'3, obtained from Eq. ()3.17|) by the 
replacement A;"*" , 



S3{r,k-) 





83(1- ^ 


= 0,k- = 0) 9{M 


+ 


[S3{1-.. 


,k~ = o)-s,ii- 


+ 




= o,k-)-S3ii- 


+ 




',k') - S3{r,k' 






= o,k-)-s,ir 



- \1-\)9{M- \k-\) 

= o,k~ = o)e{M- \i-\)]e{M- \k-\) 
= o,k" = o)e{M -\k~\)]e{M -\r\) 

= 0)^(M- |A;-|)} 

= o,k- = o)e{M -\k-\)}e{M 

(3.34) 



to treat the soft gluons with momenta k and / attached to jet J\ . The con- 
tribution from the first term in Eq. ()3.34|1 . when used in Eq. ()3.32|1 . vanishes 
since the integrand S^lk^ = k^ = 0,k_i_,l~ = 0,1^, l±) is an antisymmetric 
function of as can be easily checked using Eqs. ()2.1H) . ()2.36|) and ()B.35|) . 
We can apply the K-G decomposition on the gluon with momentum / when 

treating the second term in Eq. ()3.34|) used in Eq. ()3.32p . At LL only the 

(3) 

K gluon contributes. It can be factored from the jet function J\ with the 
result shown in Figs. 13.6b and 13. 6h . In a similar way we can treat the gluon 
with momentum k in the third term of Eq. ()3.34|) . After we factor this gluon 
from the jet J\ , we obtain the contributions shown in Figs. 13. (iH and 13. (ib . 
In the case of the last term in Eq. ()3.34p . we can factor out both soft gluons 
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with momenta k and / from jet . The result of this factorization is shown 
in Fig. 13 ■6f . 

Next, we note that the combination of the diagrams in Figs. l3.5l i. l!rHfe and 
13.6b is the same as the result encountered in the analysis of the LL amplitude, 
Fig. E31 We write 

/M 
dk- (Fig.EIll + Fig.E3l3 + Fig.ESb) = a^^\q± - k±)T'-]^^''\pA,q, k^). 
■M 

(3.35) 

where a^^\q — k) in Eq. ()3.35p is given by the diagrams in Fig. 13.31 with an 
external momentum q — k = (0~*",0~,g_L — k±). In the case when the gluon 
coming out of the boxed vertex attaches to an external gluon with momentum 
k, we evaluate the one loop trajectory a^^\k±) in Eq. ()3.35p . 

To complete the analysis, we have to discuss the diagrams in Figs. 13. (jb . 
and 13. 6b - I3.6f . In the region ~ we can factor the gluon with momentum 
/ from the jet function J^^\l'^ = 0,l~,l±) in the case of the diagram in Fig. 
13.6b . The resulting k~ and integral is over an antisymmetric function of 
k~ and Z^, and therefore it vanishes. So the only contribution comes from the 
Glauber region, where we can set 1^=0 outside Ja {l'^ = 0, As above, 

we perform the and k~ integrals in the limit M — > 0. The integrand does 
not develop a singularity in k~ and/or strong enough to compensate for the 
shrinkage of the integration region J^^^ dk~ when M — 0. Hence the diagram 
in Fig. 13.6b does not contribute in the limit M — ^ 0. In a similar way as for 
the diagram in Fig. 13.6b . none of the diagrams in Figs. 13.6b - I3.6t ' contribute. 
The diagrams in Figs. 13.6b - 13.6b vanish in the M — > limit, while in the case 
of the diagram in Fig. 13. 6t ' the k~ and integral is over an antisymmetric 
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function of k~ and 

At this point we have discussed all the contributions appearing on the right 
hand side of the evolution equation ()3.22|) . Combining the partial results given 
by Eqs. (|3.35|) and (|3.3(J|) in Eq. (j3.22p . we arrive at the evolution equation 

(2) 

governing the high energy behavior of T)^ 

(^Pa - ''\PA^q,k±) = 2a s faecfbed J ^2^)^-2 ^^^^ ^' 

( kl {k-q)l ql 



ll{k-l)l (/ - g)i (A: - Oi 



(3.36) 



Projecting out onto the color singlet in Eq. ()3.36|) . we immediately recover 
the celebrated BFKL equation [21]. 

3.2.1 Evolution of T^""^ at LL 

We can now generalize Eq. ()3.36p to the case of . The evolution kernel in 
this case contains, besides a piece diagonal in the number of external gluons, 
also contributions which relate jet functions with different number of external 
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gluons 

, d 



Pa 



2as'Y^fa,ebifajebj J (^2n)^-2 (2n)D-'i + " ^i-L ~ ^. 



p(n) ai...6i...fej...a„/ + , , , V. \ 

A 1^/;^, y, ...,(,,_!_,..., J 

n 

+ 5^ X r!,") g, . . . , 

i=l 

n-1 

+ E ...... (3.37) 

n'=l 

where ®_l denotes a convolution in transverse momentum space. The last term 
in Eq. ()3.37p corresponds to the configurations when one or more external 
gluons attach to a gluon or a ghost lines forming the one loop kernel derived 
for r^"*. Using the notation of Sec. I2.5.4[ we can write Eq. ()3.37|) at r-loop 
order in a form 

- 1) rj'^) = g/c("'"'^i) ® rj''^-^). (3.38) 



It corresponds to Eq. ()2.64p of Sec. 12.5.41 when written in terms of the coef- 
ficients Cr"'*^^ introduced in Eq. ()2.(j()|l . From Eq. ()3.38|1 we immediately see 
that the following property of the one loop kernel 

^(n,n',l) ^ _ £(nX;l)^ ^3 39) 
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is satisfied. We recall that this step was essential in demonstrating that the 
set of evolution equations, Eq. ()2.5H1 . forms a consistent system, refer to the 
paragraph above Eq. ()2.64|) . 

The term diagonal in the number of external gluons in Eq. (|3.37|) coin- 
cides with the evolution equation derived in Ref. [17]. Our formalism, besides 
enabling us to go systematically beyond LL accuracy, indicates that even at 
LL, in addition to the kernels found in Ref. [17], the kernel has contributions 
which relate jet functions with different number of external gluons. 



3.3 Gluon reggeization at NLL 

In this section, we study the NLL evolution equation and the gluon reggeiza- 
tion. In Sec. 13.3.11 we identify the contributions to the amplitude at NLL 
and negative signature channel. Sec. 13.3.21 gives a detailed derivation of the 
relevant evolution equation. 

3.3.1 Color octet and negative signature amplitude 

The statement of gluon reggeization means that the scattering amplitude in the 
Regge limit, dominated by the gluon exchange (color octet) in the t-channel, 
and antisymmetric under s ^ u (negative signature), takes the form 



^8 ('^5 ^) rA_rB,rir2 



rA 



T2rB 1 



(3.40) 



where -^^(t) and F^{t) are the impact factors, which depend on the properties 
of the scattered particles. The Regge trajectory a{t) determines the depen- 
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Figure 3.7: Projection of the amplitude onto the negative signature channel. 

dence on the energy ^/s. It is a universal function of t meaning that it does 
not depend on the particles involved in the scattering. In order to project onto 
the negative signature channel, besides the amplitude for qq' qq' scattering, 
Eq. ()2.1|) . we consider the amplitude for the process 

q{pA, ri, Xa) + q'iPB, rB, \b) q{PA - q, rA, Ai) + q'{pB + q, r2, A2). (3.41) 

The amplitude for scattering (j3.41|) . A{s,t), can be obtained from the ampli- 
tude for the process ()2.1|) . A{s,t), by the crossing symmetry s ^ u. This 
means that A{s,t) = A{u,t) ~ A{—s,t). Hence, if we define 

A^-\s,t)^^ (A(s,t)-i(s,t)), (3.42) 

as indicated in Fig. 13.71 the amplitude A^~^ will have a negative signature by 
construction. Since s ~ —u in the Regge limit, the logarithmic derivative of 
this amphtude is: 

dA^-\ , 1 / , dA^ ^\ , ^ , 

= 2 («g^(M)-.g^(M)j. (3-43) 

We now isolate the negative signature contribution to the amplitude for 



77 



process at NLL. The projection of the set of two gluon exchange diagrams, 
Fig. Kllb . onto the negative signature takes the form: 

X X J^^^"^(A;- = 0,A;+,fc^), (3.44) 

where S{k) has been defined in Eq. (jH.K-ij) and J\ is the jet function corre- 
sponding to the anti-quark in the process ()3.41|) moving in the plus direction. 
In Appendix IB. II we show that the jet functions J^""* and J^"'' are related by 
the symmetry: 

(3.45) 

in the Regge limit. Applying this identity to the jet function with two external 
gluons (n = 2), the helicity conserving part of the jet function (A^ -|- Ai = 
2Xa = ±1) obeys: 

j'^\k+ = 0, -k-, k^) = f^\k+ = 0, k-, k^). (3.46) 

We apply this result to analyze Eq. ()3.44|) . Using Eqs. ()3.16|) and ()3.46|) . 
we easily see that in the Glauber region fc~ ^ ~ k± the integrand in Eq. 
()3.44|1 is over an antisjnumetric function of k~ and therefore the k~ integral 
vanishes. This indicates that there is no contribution from the Glauber region 
when analyzing a soft gluon with momentum k in Fig. 13.2b and therefore this 

(2) (2) 

gluon can be factored from the jet functions J\ and J^' . As a consequence, 
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the high-energy behavior of the amphtude for the process ()2.H) . projected 
onto the negative signature channel, is determined by the jet functions with 
one external gluon only, J^^^ and J^j^\ It also means that ^ is automatically 
in the color octet channel. 

Combining the contributions from the one and the two gluon exchange 
diagrams shown in Fig. 13. ![ the negative signature amplitude A^g~^ takes the 
form: 

X jiS^''iPA,q,v), (3.47) 

where n+_(g, 77) is the one-loop gluon self-energy and a^^\t) is the one-loop 
contribution to the gluon trajectory given in Eq. ()3.12|1 . Eq. ()3.47|) indicates 
that in order to determine the high energy behavior of the scattering ampli- 
tude, we need to determine the high-energy behavior of J^^ and If we 
show that J^^ and J^^ satisfy: 



n 7(1) 

= ^it)JA\ (3.48) 

at NLL, then the solution to the evolution equation ()3.48|1 used in Eqs. ()3.43j) 
and ()3.47|1 implies the Regge ansatz, Eq. ()3.40|) . This would provide a proof 
of gluon reggeization at NLL to all orders in perturbation theory in QCD, 
Ref. [28]. We study the evolution equation ()3.48j) in the next section. 
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at NLL. 



3.3.2 NLL evolution equation for 

We now derive the evolution equation for jet function J^^"" at NLL. We restrict 
ourselves to J^^ since J^^ can be analyzed in the same way. Applying the 
general evolution equation, Eq. (j2.5ip . for n = 1 and using the results in Fig. 
\2A\ we arrive at the contributions to {d/d \np\ — 1)J^^^^ shown in Fig. 13.81 
The diagram in Fig. I3.8l i. is already in a factorized form so it does not need to 
be analyzed further. The blob in the bottom part of this diagram represents 
all possible two loop corrections. 

In Fig. 13.8b . the bottom dark blob denotes the one loop gluon-self energy. 
The gluon with momentum k can be analyzed the same way as in the case of 
LL, Sec. 13.11 with a result in a factorized form. 

In order to prove reggeization of a gluon at NLL, we need to show that the 
soft gluon with momentum k in Fig. 13.8b and the soft gluons with momenta 
k and I in Fig. 13.8b can be factored from the jet functions J^^"* and 
respectively. Let us now analyze these particular diagrams. The contribution 
to Fig. 13.8b can be written 

Fig.^E^ = -igsfabc j -^:^S^y{k,q) f^2^'ll{pA,q-,k). (3.49) 
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a) b) c) 

Figure 3.9: Contributions to Fig. 13.8b from the K gluon at NLL. 
Next, we use the following identity for the integrand in Eq. ()3.49p 



+ 
+ 



S+ + {k) jlJ)'l)c(k) ~ S^ + {k) j[}-^Xc^k^ — 0, 



S,Ak) jflTM - S^^k) J[%::{k) . (3.50) 



The first term in Eq. (j3.5(jp . after being inserted into Eq. (j3.49|) . can be 
analyzed using the K-G decomposition, since there is no contribution from 
the Glauber region due to the antisymmetry of S+ + {k~^ , k~ = 0,k±) under 
the transformation k~^ — ^ —k^. The contributions from the K gluon, after 
applying the Ward identities, can be expressed in a form shown in Fig. 13.91 
The diagram in Fig. 13.9b is the same as in the case of LL. At NLL, we also have 
to insert the gluon and ghost self-energies where appropriate. This diagram 
is in a factorized form and therefore consistent with a gluon reggeization at 
NLL. 

As for diagrams in Figs. 13.9b and 13.9b . the gluon with momentum / can 
be factored from the jet function J^^ as shown in Appendix IB. 21 Thus the 
contributions from the K gluon in Figs. 13.9b - 13.9b can be written in a factor- 
ized form and they are in an agreement with the gluon reggeization at NLL. 
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a) b) c) 




d) e) f) 

Figure 3.10: Soft contributions to Figs. 13.8b and 13.8b . 

Next, we analyze the contributions from the G gluon with momentum k cor- 
responding to the first term in Eq. ()3.50p and to Fig. 13.8b . They come from 



the diagrams shown in Figs. l3.1(Jb . l3.1(Jb and l3.1(Jl i since the G gluon must 
be attached to a soft line. In Appendix IB. 21 we show that the soft gluons 
with momenta / and k — I can be factored from the jet function J\ in Fig. 
13.10b . the soft gluons with momenta / can be factored from the jet functions 

(rx\ _______ _______ 

J\ in Figs. 13.10b and I3.10d and therefore the contribution from the G gluon 
corresponding to Fig. 13.8b is in a factorized form. 

Now, we examine the gluons with momenta k and/or q — k m the second 
and the third terms of Eq. ()3.50p . after used in Eq. ()3.49|) . They must attach 
to soft lines and therefore the typical contributions to these terms come from 
diagrams shown in Figs. 13.10b - I3.10b l. First we analyze the contributions to 
the second term in Eq. ()3.50j) . As shown in Appendix IB.2[ the contribution 
to the diagram in Fig. 13.10b is in a factorized form. However, there are 
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contributions to the diagrams in Fig. 13.10b and Fig. 13.10b . coming from 
the region when all soft gluons external to J\ are Glauber, which cannot be 
written in a simple factorized form. These contributions, Eq. ()B.20|) and Eq. 
f)B.24|) . are rather in a form of a convolution in transverse momentum space. 

As far as the last term in Eq. ()3.50|) is concerned, its contributions come 
from the diagrams in Figs. 13.10b - I3.10H . As we demonstrate in Appendix 
IB.2| the soft gluons external to jjf ^ in Figs. I3.1()k and 13.10b . as well as the 
soft gluons external to J^^"* in Figs. 13.10b and I3.10t i can be factored from the 
the jet functions. 

Finally, we need to examine the diagram in Fig. 13.8b which can be written 

as 

X '""{pA, q; = 0,k', k^, /+ = 0, 

+ J -^SUl. q) J^I^ '\PA, q; = 0, h) 

+ J ^^SUK q) Jf '\VA, q; k+ = 0, fc-, k^). (3.51) 

The contribution to the first term in Eq. ()3.51|) comes from a diagram shown 
in Fig. 13.8b where the soft gluons with momenta k, I and q — k — l attach to jet 
lines inside the jet function At NLL we can set fc+ = /+ = inside jf 
and also the soft gluons attach to the plus components of the jet's vertices. As 
shown in Appendix IB.2[ in the region when both gluons with momenta k and 
/ are Glauber, we cannot factor out these gluons and the contribution from 
this double Glauber region is not in a factorized form. 
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The contribution to the last two terms in Eq. ()3.51|) comes from the graphs 
shown in Figs. KllOb and IH.lOf . respectively. The gluon with momentum 
q — k — l attaches to soft gluons with momenta fc, Fig. 13.10b . and /, Fig. I3.10t '. 
respectively. In Appendix IB. 2 1 we demonstrate that the gluons external to J^^\ 
can be factored from the jet functions J\ . 

As a result of the previous reasoning, we arrive at a conclusion that the 
evolution equation for the jet function J^^'' at NLL takes the form: 

X vf'"'\pA,q-k^,l^). (3.52) 

The kernel /C*^^'^-* is given by all the diagrams where it was possible to factor 
the soft gluons from the jet functions. This piece is in an agreement with 
gluon reggeization at NLL according to Eq. ()3.48|) . Besides this term, there 
is a contribution in Eq. ()3.52j) corresponding to the configurations when all 
soft gluons external to the jet function J\ are Glauber. Combining the re- 
sults of Eqs. flR20|l . (Ib:24I) and (|R39|l in Appendix |B21 we can find out the 
explicit form of the kernel /C^^'^-*. The presence of this term in Eq. ()3.52j) 
suggests the breakdown of gluon reggeization at NLL level. We have to note 
the following, however. At two loop level, these non-factorizing terms vanish. 
The non-factorizing piece in every diagram comes from the double pinched 
region, corresponding to the two delta function terms originating from the 
two denominators on the fermion line present in the tree-level jet function. 
They, however, do not produce any dependence on transverse momenta and 
therefore this term is symmetric in external momenta. Hence if we inter- 
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change the two external momenta flowing out of J\ and into the three point 
gluon-gluon-gluon or ghost-gluon-ghost vertex, Figs. I3.10b . imib and 13. 8b . the 
integrand becomes an antisymmetric function under this exchange and the in- 
tegral vanishes. However, when the jet function contains loop corrections, it 
will, in general, not be symmetric under the interchange of external momenta 
only. In order to get the symmetry, we need to interchange the colors as well. 
Nevertheless, the above argument shows that the standard analysis and the 
calculations which have been performed at the two loop level, Refs. [26] - [33], 
are not in a contradiction with our results. 

3.4 Conclusions 

As an illustration of the general algorithm we have demonstrated it in an 
action at NLL in both the amplitude and the evolution equations. First, in 
Sec. 13.11 we have performed the resummation of the amplitude at LL. We 
have found the standard expression for the gluon Regge trajectory, Eq. ()3.12|1 . 
Then, in Sec. 13.21 we have identified the class of diagrams contributing to the 
amplitude at NLL level. They involved one and two gluon exchange diagrams. 
Fig. 13.11 The resummation of the two gluon exchange diagram in Sec. 13. 2[ 
needed to be performed at LL. It has lead to the celebrated BFKL equation, 
Eq. ()3.3fij) . The one gluon exchange contribution had to be resummed at 
NLL. Together with this resummation, we have addressed the question of 
gluon reggeizaiton at NLL in Sec. 13.31 We have found that the majority of the 
terms in the evolution kernel are in an agreement with this conjecture to all 
orders in perturbation theory. However, we have also identified contributions 
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violating the Regge ansatz starting at the three loop level and calculated the 
corresponding evolution kernels, Eqs. ()B.20jl . ()B.24|1 . ()B.39j) . We have not 
been able to prove that these contributions decouple from the jets. A further 
study of symmetries involving jets with three external gluons should shed more 
light on the presence or absence of the non-factorized terms, Ref. [11]. 
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Part II 

Resummation in Dijet Events 
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Chapter 4 

Event Shape / Energy Flow 
Correlations 

4.1 Introduction 

The agreement of theoretical predictions with experiment for jet cross sections 
is often impressive. This is especially so for inclusive jet cross sections at high 
Pt, using fixed-order factorized perturbation theory and parton distribution 
functions [52]. A good deal is also known about the substructure of jets, 
through the theoretical and experimental study of multiplicity distributions 
and fragmentation functions [53], and of event shapes [54-56]. Event shape 
distributions [57-59] in particular offer a bridge between the perturbative, 
short-distance and the nonperturbative, long-distance dynamics of QCD [60] . 

Energy flow [61] into angular regions between energetic jets gives informa- 
tion that is in some ways complementary to what we learn from event shapes. 
In perturbation theory, the distribution of particles in the final state reflects 
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interference between radiation from different jets [53], and there is ample evi- 
dence for perturbative antenna patterns in interjet radiation at both e+e" [62] 
and hadron coUiders [63,64]. Energy flow between jets must also encode the 
mechanisms that neutralize color in the hadronization process, and the tran- 
sition of QCD from weak to strong coupling. Knowledge of the interplay 
between energy and color flows [65,66] may help identify the underlying event 
in hadron collisions [67], to distinguish QCD bremsstrahlung from signals of 
new physics. Nevertheless, the systematic computation of energy flow into 
interjet regions has turned out to be subtle [68] for reasons that we will review 
below, and requires a careful construction of the class of jet events. It is the 
purpose of this work to provide such a construction, using event shapes as a 
tool. 

In this paper, we introduce correlations between event shapes and energy 
flow, "shape/flow correlations", that are sensitive primarily to radiation from 
the highest-energy jets. So long as the observed energy is not too small, in 
a manner to be quantified below, we may control logarithms of the ratio of 
energy flow to jet energy [66,69]. 

The energy flow observables that we discuss below are distributions asso- 
ciated with radiation into a chosen interjet angular region, fl. Within Q we 
identify a kinematic quantity Qn = sQ, at cm. energy Q, with s ^ 1. Qa 
may be the sum of energies, transverse energies or related observables for the 
particles emitted into Q. Let us denote by Q the complement of Q. We are 
interested in the distribution of Qci for events with a fixed number of jets in 
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This set of events may be represented schematically as 

A + B^ Jets +X^ + Rn{Qn)- (4.1) 

Here Xq stands for radiation into the regions between Vt and the jet axes, and 
Rq_ for radiation into VL. 

The subtlety associated with the computation of energy flow concerns the 
origin of logarithms, and is illustrated by Fig. 14.11 Gluon 1 in Fig. 14.11 is an 
example of a primary gluon, emitted directly from the hard partons near a jet 
axes. Phase space integrals for primary emissions contribute single logarithms 
per loop: (l/Qf^X ln"-^(Q/Qn) = (l/eQ)a"ln""^(l/e), n > I, and these 
logarithms exponentiate in a straightforward fashion [66]. At fixed for 
Eq. ()4.1|) . however, there is another source of potentially large logarithmic 
corrections in Qq. These are illustrated by gluon 2 in the figure, an example 
of secondary radiation in f2, originating a parton emitted by one of the leading 
jets that define the event into intermediate region VL. As observed by Dasgupta 
and Salam [68], emissions into Vt from such secondary partons can also result 
in logarithmic corrections, of the form (l/(5n)«" ln"~^(Qn/<5f7), n>2^ where 

is the maximum energy emitted into VL. These logarithms arise from strong 
ordering in the energies of the primary and secondary radiation because real 
and virtual enhancements associated with secondary emissions do not cancel 
each other fully at fixed Q^. 

If the cross section is fully inclusive outside of fi, so that no restriction 
is placed on the radiation into fi, can approach Q, and the secondary 
logarithms can become as important as the primary logarithms. Such a cross 
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section, in which only radiation into a fixed portion of phase space (Q) is 
specified, was termed "non-global" by Dasgupta and Salam, and the associated 
logarithms are also called non-global [68,70,71]. 

In effect, a non-global definition of energy fiow is not restrictive enough 
to limit final states to a specific set of jets, and non-global logarithms are 
produced by jets of intermediate energy, emitted in directions between region 
Q and the leading jets. Thus, interjet energy fiow does not always originate 
directly from the leading jets, in the absence of a systematic criterion for 
suppressing intermediate radiation. Correspondingly, non-global logarithms 
reflect color flow at all scales, and do not exponentiate in a simple manner. 
Our aim in this paper is to formulate a set of observables for interjet radiation 
in which non-global logarithms are replaced by calculable corrections, and 
which reflect the flow of color at short distances. By restricting the sizes of 
event shapes, we will limit radiation in region Q, while retaining the chosen 
jet structure. 

An important observation that we will employ below is that non-global 
logarithms are not produced by secondary emissions that are very close to a 
jet direction, because a jet of parallel-moving particles emits soft radiation 
coherently. By flxing the value of an event shape near the limit of narrow 
jets, we avoid flnal states with large energies in Cl away from the jet axes. At 
the same time, we will identify limits in which non-global logarithms reemerge 
as leading corrections, and where the methods introduced to study nongobal 
effects in Refs. [68,70,71] provide important insights. 

To formalize these observations, we study below correlated observables for 
e~^e~ annihilation into two jets. (In Eq. (j4.Hl A and B denote positron and 
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Figure 4.1: Sources of global and non-global logarithms in dijet events. Config- 
uration 1, a primary emission, is the source of global logarithms. Configuration 
2 can give non-global logarithms. 

electron.) In e^e~ annihilation dijet events, the underlying color flow pattern 
is simple, which enables us to concentrate on the energy flow within the event. 
We will introduce a class of event shapes, /(a) suitable for measuring energy 
flow into only part of phase space, with a an adjustable parameter. To avoid 
large non-global logarithmic corrections we weight events by exp[— z//], with 
V the Laplace transform conjugate variable. 

For the restricted set of events with narrow jets, energy flow is proportional 
to the lowest-order cross section for gluon radiation into the selected region. 
The resummed cross section, however, remains sensitive to color flow at short 
distances through anomalous dimensions associated with coherent interjet soft 
emission. In a sense, our results show that an appropriate selection of jet events 
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automatically suppresses nonglobal logarithms, and confirms the observation 
of coherence in interjet radiation [53,63]. 

In the next section, we introduce the event shapes that we will correlate 
with energy flow, and describe their relation to the thrust and jet broadening. 
Section 3 contains the details of the factorization procedure that characterizes 
the cross section in the two-jet limit. This is followed in Sec. 4 by a derivation 
of the resummation of logarithms of the event shape and energy flow, following 
the method introduced by Collins and Soper [72]. We then go on in Sec. 5 
to exhibit analytic results at leading logarithmic accuracy in Qq/Q and next- 
to-leading logarithm in the event shape. Section 6 contains representative 
numerical results. We conclude with a summary and a brief outlook on further 
applications. 

4.2 Shape/Flow Correlations 

4.2.1 Weights and energy flow in dijet events 

In the notation of Eq. ()4.H) . we will study an event shape distribution for the 
process 

e+ + e- ^ J,{pjJ + .Upj^) + Xn (/) + RniQn) , (4.2) 

at cm. energy Q ^ Qq ^ ^qco- Two jets with momenta pj^, c = 1, 2 emit 
soft radiation (only) at wide angles. Again, f2 is a region between the jets 
to be specified below, where the total energy or the transverse energy Qq of 
the soft radiation is measured, and denotes the remaining phase space (see 
Fig. 14.11) . Radiation into Cl is constrained by event shape /. We refer to cross 
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sections at fixed values (or transforms) of / and Qq as shape/flow correlations. 

To impose the two-jet condition on the states of Eq. ()4.2j) we choose weights 
that suppress states with substantial radiation into Cl away from the jet axes. 
We now introduce a class of event shapes /, related to the thrust, that enforce 
the two-jet condition in a natural way. 

These event shapes interpolate between and extend the familiar thrust [55] 
and jet broadening [58,59], through an adjustable parameter a. For each state 
N that defines process ()4.2|) . we separate into two regions, 0^, c = 1,2, 
containing jet axes, nc{N). To be specific, we let Cli and ^2 be two hemispheres 
that cover the entire space except for their intersections with region Q. Region 
Cli is centered on ni, and Q2 is the opposite hemisphere. We will specify the 
method that determines the jet axes fii and ^2 momentarily. To identify a 
meaningful jet, of course, the total energy within Cli should be a large fraction 
of the available energy, of the order ofQ/2 in dijet events. In e^e^ annihilation, 
if there is a well-collimated jet in f^i with nearly half the total energy, there 
will automatically be one in O2. 

We are now ready to define the contribution from particles in region Qc to 
the a-dependent event shape, 

/n,(iV,a) = -^ kl^ujf^ {l~n,.h,)'~\ (4.3) 

where a is any real number less than two, and where ^/s = Q is the cm. 
energy. The sum is over those particles of state N with direction that flow 
into Oc, and their transverse momenta ki^± are measured relative to fic- The 
jet axis ni for jet 1 is identified as that axis that minimizes the specific thrust- 
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related quantity fn^{N,a = 0). When flc in Eq. ()4.3p is extended to all of 
phase space, the case a = is then essentially 1 — T, with T the thrust, while 
a = 1 is related to the jet broadening. 

Any choice a < 2 in (j4.3|) specifies an infrared safe event shape variable, 
because the contribution of any particle i to the event shape behaves as 9'^~"' 
in the collinear limit, 6i = cos~^{ni ■ fic) 0. Negative values of a are clearly 
permissible, and the limit a —oo corresponds to the total cross section. 
At the other limit, the factorization and resummation techniques that we 
discuss below will apply only to a < 1. For a > 1, contributions to the event 
shape ()4.3j) from energetic particles near the jet axis are generically larger than 
contributions from soft, wide-angle radiation, or equal for a = 1. When this is 
the case, the analysis that we present below must be modified, at least beyond 
the level of leading logarithm [59] . 

In summary, once ni is fixed, we have divided the phase space into three 
regions: 

• Region Q, in which we measure, for example, the energy fiow, 

• Region Cli, the entire hemisphere centered on -hi, that is, around jet 1, 
except its intersection with Q, 

• Region ^2, the complementary hemisphere, except its intersection with 

n. 

In these terms, we define the complete event shape variable f{N, a) by 

fiN,a) = f^^{N,a) + f^^{N,a), (4.4) 
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with /q^, c = 1, 2 given by ()4.3|) in terms of the axes ni of jet 1 and n2 of jet 
2. We will study the correlations of this set of event shapes with the energy 
flow into Q, denoted as 

f{N) = ^ 5^ c^. . (4.5) 

The differential cross section for such dijet events at fixed values of / and 
/ is now 



V 

x6{e - f{N)) 6{e - f{N, a)) 6\n, - h{N)) ,(4.6) 



de de dm 2s 

N 



where we sum over all final states N that contribute to the weighted event, 
and where M{N) denotes the corresponding amplitude for e^e~ N. The 
total momentum is pi, with p'j = s = Q"^. As mentioned in the introduction, 
for much of our analysis, we will work with the Laplace transform of ()4.6|) . 



da{e,h',s,a) j- -ue ^'^i^^^^ ^y^) 



de dhi Jq de de dhi 

= |M(iV)pe-^^"(^''^) {2n)U\pj-p^) 

x6{e-f{N))S\n,-n{N)). (4.7) 

Singularities of the form (l/e) ln"(l/e) in the cross section ()4.6j) give rise to 
logarithms In""*"^ u in the transform ()4.7p . 

Since we are investigating energy flow in two-jet cross sections, we fix the 
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constants e and e to be both much less than unity: 

(4.8) 

We refer to this as the elastic limit for the two jets. In the elastic limit, the 
dependence of the directions of the jet axes on soft radiation is weak. We will 
return to this dependence below. Independent of soft radiation, we can always 
choose our coordinate system such that the transverse momentum of jet 1 is 
zero, 

PJi,± = 0, (4.9) 

with pj-^ in the X3 direction. In the limit e, e — 0, and in the overall cm., pj^ 
and then approach light-like vectors in the plus and minus directions: 




(4.10) 



As usual, it is convenient to work in light-cone coordinates, = {p'^ ,p~ ,p±), 
which we normalize as = {l/^/2){p^ ±p^). For small e and e, the cross 
section ()4.6|) has corrections in ln(l/e) and ln(l/e), which we will organize in 
the following. 

4.2.2 Weight functions and jet shapes 

In Eq. ()4.3|) . a is a parameter that allows us to study various event shapes 
within the same formalism; it helps to control the approach to the two-jet 
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limit. As noted above, a < 2 for infrared safety, although the factorization that 
we will discuss below applies beyond leading logarithm only to 1 > a > — oo. 
A similar weight function with a non-integer power has been discussed in a 
related context for 2 > a > 1 in [73] . To see how the parameter a affects the 
shape of the jets, let us reexpress the weight function for jet 1 as 

fn,{N,a) = ^ V a;,sin«^,(l-cos^,)'~", (4-11) 

where 9i is the angle of the momentum of final state particle i with respect to 
jet axis rii. As a — > 2 the weight vanishes only very slowly for 9^ — > 0, and at 
fixed /f^^, the jet becomes very narrow. On the other hand, as a ^ — oo, the 
event shape vanishes more and more rapidly in the forward direction, and the 
cross section at fixed /q^ becomes more and more inclusive in the radiation 
into Qi. 

In this paper, as in Ref. [66], we seek to control corrections in the single- 
logarithmic variable as{Q)l'n{l/e), with e — Qn/Q. Such a resummation is 
most relevant when 

a.W)ln(i)>l^.<exp(-±^). (4.12) 

Let us compare these logarithms to non-global effects in shape/flow correla- 
tions. At z/ = and for a — — oo, the cross section becomes inclusive outside 
fl. As we show below, the non-global logarithms discussed in Refs. [66, 68] 
appear in shape/flow correlations as logarithms of the form as{Q) ln(l/(£z/)), 
with i> the moment variable conjugate to the event shape. To treat these 
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logarithms as subleading for small e and (relatively) large we require that 
aM) lu (i;) < 1 - . > i exp (^) . (4,13) 

For large z/, there is a substantial range of e in which both ()4.12j) and ()4.13p 
can hold. When v is large, moments of the correlation are dominated precisely 
by events with strongly two-jet energy flows, which is the natural set of events 
in which to study the influence of color flow on interjet radiation. (The peak 
of the thrust cross section is at (1 — T) of order one-tenth at LEP energies, 
corresponding to v of order ten, so the requirement of large v is not overly 
restrictive.) In the next subsection, we show how the logarithms of {sv)^^ 
emerge in a low order example. This analysis also assumes that a is not large 
in absolute value. The event shape at fixed angle decreases exponentially with 
a, and we shall see that higher-order corrections can be proportional to a. We 
always treat \nv as much larger than \a\. 

4.2.3 Low order example 

In this section, we check the general ideas developed above with the concrete 
example of a two-loop cross section for the process ()4.2|) . This is the lowest 
order in which a non-global logarithm occurs, as observed in [68] . We normalize 
this cross section to the Born cross section for inclusive dijet production. A 
similar analysis for the same geometry has been carried out in [68] and [74]. 

The kinematic configuration we consider is shown in Fig. 14.21 Two fast 
partons, of velocities (3i and /32, are treated in eikonal approximation. In 
addition, gluons are emitted into the final state. A soft gluon with momentum 
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Figure 4.2: A kinematic configuration that gives rise to the non-global loga- 
rithms. A soft gluon with momentum k is radiated into the region Q, and an 
energetic gluon with momentum I is radiated into Q. Four- vectors /5i and (32, 
define the directions of jet 1 and jet 2, respectively. 

k is radiated into region Q and an energetic gluon with momentum / is emitted 
into the region Cl. We consider the cross section at fixed energy, uj^ = £\/s- 
As indicated above, non-global logarithms arise from strong ordering of the 
energies of the gluons, which we choose as ui ^ uJk- In this region, the gluon 
/ plays the role of a "primary" emission, while is a "secondary" emission. 

For our calculation, we take the angular region Q to he a "slice" or "ring" 
in polar angle of width 26, or equivalently, (pseudo) rapidity interval (— 77,77), 
with 

A, = 2, = lnfi±^). (4.14) 

\ 1 — sm d / 

The lowest-order diagrams for this process are those shown in Fig. 14.31 includ- 
ing distinguishable diagrams in which the momenta k and / are interchanged. 

The diagrams of Fig. 14.31 give rise to color structures Cp and CpCa, but 
terms proportional to Cp may be associated with a factorized contribution to 
the cross section, in which the gluon k is emitted coherently by the combina- 
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I I I 

a) b) c) 




d) e) f) 



Figure 4.3: The relevant two-loop cut diagrams corresponding to the emission 
of two real gluons in the final state contributing to the eikonal cross section. 
The dashed line represents the final state, with contributions to the amplitude 
to the left, and to the complex conjugate amplitude to the right. 
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tions of the gluon / and the eikonals. To generate the CfCa part, on the other 
hand, gluon k must "resolve" gluon / from the eikonal lines, giving a result 
that depends on the angles between I and the eikonal directions. 

The computation of the diagrams is outlined in Appendix IC.ll here we 
quote the results. We adopt the notation q = cos 6i, si = sin 6*;, with 9i the 
angle of momentum / measured relative to (3i, and similarly for k. We take, as 
indicated above, a Laplace transform with respect to the shape variable, and 
identify the logarithm in the conjugate variable u. We find that the logarithmic 
C^Cyi-dependence of Fig. l4.Hl mav be written as a dimensionless eikonal cross 
section in terms of one energy and two polar angular integrals as 

sin 



^ -T ^ \ I 

ae \ TT / e 



X 



in (5 J sin 5 

^ (4.15) 



1/1 1 \ 1 1 

+ 



Ck + ci 1 + Ck \ l + ci ' I - CkJ sll + Ci 



In this form, the absence of collinear singularities in the CfCa term at cos 6'; = 
+1 is manifest, independent of v. Collinear singularities in the / integral 
completely factorize from the k integral, and are proportional to C'j^. The 
logarithmic dependence on e for > 1 is readily found to be 

C,CJ^)'1\J1-\C(1^,), (4.16) 

V TT / e \eiy J 



de 

where C{Ar]) is a finite function of the angle 5, given explicitly in Appendix 

We can contrast this result to what happens when i/ = 0, that is, for an 
inclusive, non-global cross section. In this case, recalling that e = Qn/Q, we 
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find in place of Eq. ()4.16p the non-global logarithm 

As anticipated, the effect of the transform is to replace the non-global loga- 
rithm in Q/Qn, by a logarithm of l/leu). We are now ready to generalize this 
result, starting from the factorization properties of the cross section near the 
two- jet limit. 



4.3 Factorization of the Cross Section 

In this section we study the factorization of the correlations ()4.6p . The analysis 
is based on a general approach that begins with the all-orders treatment of 
singularities in perturbative cross sections [75,76], and derives factorization 
from the analyticity and gauge properties of high energy Green functions and 
cross sections [8]. The functions that appear in factorized cross sections are 
expressible in terms of QCD matrix elements [77], and the matrix elements 
that we will encounter are familiar from related analyses for heavy quark and 
jet production [78]. We refer in several places below to standard arguments 
discussed in more detail in [76,8]. The aim of this section, and the reason 
why a careful analysis is necessary, is to identify the specific dimensionless 
combinations of kinematic variables on which the factorized matrix elements 
may depend. We will use these dependences in the following section, when we 
discuss the resummation properties of our correlations. 
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4.3.1 Leading regions near the two-jet limit 

In order to resum logarithms of e and e (or equivalently the Laplace con- 
jugate of e) we have first to identify their origin in momentum space when 
£,£ 0. Following the procedure and terminology of [75], we identify "lead- 
ing regions" in the momentum integrals of cut diagrams, which can give rise to 
logarithmic enhancements of the cross section associated with lines approach- 
ing the mass shell. Within these regions, the lines of a cut diagram fall into 
the following subdiagrams: 

• A hard-scattering, or "short-distance" subdiagram H, where all compo- 
nents of line momenta are far off-shell, by order Q. 

• Jet subdiagrams, Ji and J2, where energies are fixed and momenta are 
collinear to the outgoing primary partons and the jet directions that 
emerge from the hard scattering. (For £ = £ = 0, the sum of all energies 
in each jet is one-half the total energy.) To characterize the momenta of 
the lines within the jets, we introduce a scaling variable, A <C 1. Within 
jet 1, momenta I scale as {t^ ~ Q,^" ~ \Q,li_ ~ A^/^Q). 

• A soft subdiagram, S connecting the jet functions Ji and J2, in which the 
components of momenta k are small compared to Q in all components, 
scaling as {k^ ~ AQ, k±_ ~ \Q). 

An arbitrary final state N is the union of substates associated with these 
subdiagrams: 

N^Ns®Nj,® Nj, . (4.18) 
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As a result, the event shape / can also be written as a sum of contributions 
from the soft and jet subdiagrams: 



f{N, a) 



/^(7V„ a) + fliNj,, a) + fl{Nj,, a) . 



(4.19) 



The superscript reminds us that the contributions of final-state particles 
associated with the soft and jet functions depend implicitly on the full final 
state, through the determination of the jet discussed in Sec. 2. In 

contrast, the energy flow weight, f{N), depends only on particles emitted at 
wide angles, and is hence insensitive to coUinear radiation: 



When we sum over all diagrams that have a fixed final state, the contribu- 
tions from these leading regions may be factorized into a set of functions, each 
of which corresponds to one of the generic hard, soft and jet subdiagrams. The 
arguments for this factorization at leading power have been discussed exten- 
sively [72,8,79]. The cross section becomes a convolution in £, with the sums 
over states linked by the delta function which fixes fii, and by momentum 



f{N) = f{N,) . 



(4.20) 
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conservation, 



da{e, e, s, a) 



daQ 



de de dh\ 



dfii 



x6{es-F{Ns,a)) 



c=l 



X (2nY6\pi-p{Nj,)-p{Nj,)-p{N,)) 



X 6'^{hi - h{N)) 6{6 - - - 6s) 



da< 







6{e) 6{e) + . 



(4.21) 



dn 



1 



Here dao/dhi is the Born cross section for the production of a single par- 
ticle (quark or antiquark) in direction ni, while the short-distance function 
H{s, hi) = 1 + 0{as), which describes corrections to the hard scattering, is an 
expansion in with finite coefficients. The functions J'c{NjJ, S{Ns) describe 
the internal dynamics of the jets and wide-angle soft radiation, respectively. 
We will specify these functions below. We have suppressed their dependence 
on a factorization scale. Radiation at wide angles from the jets will be well- 
described by our soft functions S{Ns), while we will construct the jet functions 
J'c{NjJ to be independent of e, as in Eq. (|4.2ip . 

So far, we have specified our sums over states in Eq. ()4.2H) only when 
all lines in A^^ are soft, and all lines in Nj^ have momenta that are collinear, 
or nearly collinear to pj^. As e and e vanish, these are the only final-state 
momenta that are kinematically possible. Were we to restrict ourselves to 
these configurations only, however, it would not be straightforward to make 
the individual sums over A^^ and Nj^ infrared safe. Thus, it is necessary 
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to include soft partons in Ns that are emitted near the jet directions, and 
soft partons in the Nj^, at wide angles. We will show below how to define 
the functions Jc{Nj^), S{Ns) so that they generate factoring, infrared safe 
functions that avoid double counting. We know on the basis of the arguments 
of Refs. [72,8,79] that corrections to the factorization of soft from jet functions 
are suppressed by powers of the weight functions e and/or e. 

4.3.2 The factorization in convolution form 

Although formally factorized, the jet and soft functions in Eq. 1)4.211) are still 
linked in a potentially complicated way through their dependence on the jet 
axes. Our strategy is to simplify this complex dependence to a simple convo- 
lution in contributions to e, accurate to leading power in e and e. 

First, we note that the cross section of Eq. ()4.21)1 is singular for vanishing 
e and e, but is a smooth function of s and hi. We may therefore make any 
approximation that changes s and/or hi by an amount that vanishes as a 
power of £ and e in the leading regions. 

Correspondingly, the amplitudes for jet c are singular in ej^, but depend 
smoothly on the jet energy and direction, while the soft function is singular 
in both e and Ss, but depends smoothly on the jet directions. As a result, at 
fixed values of e and e we may approximate the jet directions and energies by 
their values at £ = e = in the soft and jet functions. 

Finally, we may make any approximation that affects the value of e and/or 
Ej^ by amounts that vanish faster than linearly for £ — * 0. It is at this stage 
that we will require that a < 1. 
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With these observations in mind, we enumerate the replacements and ap- 
proximations by which we reduce Eq. ()4.2H1 . while retaining leading-power 
accuracy. 

1. To simplify the definitions of the jets in Eq. (j4.2H) . we make the replace- 
ments f^^{Nj^,a) fc{Nj^,a) with 

fl{Nj^,a)^^ E fc-±^'~'^(l-^.-^c)'"V (4.22) 

all hi€Nj^ 

The jet weight function fc{Nj^,a) now depends only on particles asso- 
ciated with Nj^. The contribution to fc{Nj^,a) from particles within 
region Clc, is exactly the same here as in the weight ()4.3|) . but we now in- 
clude particles in all other directions. In this way, the independent sums 
over final states of the jet amplitudes will be naturally infrared safe. The 
value of fc{Nj^, a) differs from the value of {Nj^, a), however, due to 
radiation outside Clc, as indicated by the new subscript. This radiation 
is hence at wide angles to the jet axis. In the elastic limit ()4.8|) . it is 
also constrained to be soft. Double counting in contributions to the total 
event shape, f{N, a), will be avoided by an appropriate definition of the 
soft function below. The sums over states are still not yet fully indepen- 
dent, however, because the jet directions fic still depend on the full final 
state N. 

2. Next, we turn our attention to the condition that fixes the jet direction 
ni. Up to corrections in the orientation of "hi that vanish as powers of e 



108 



and e, we may neglect the dependence of hi on Ng and Nj^: 

6{hi - h{N)) 6ihi - h{NjJ) . (4.23) 

In Appendix 1021 we show that this replacement also leaves the value of e 
unchanged, up to corrections that vanish as e^"". Thus, for a < 1, (j4.23j) 
is acceptable to leading power. For a < 1, we can therefore identify the 
direction of jet 1 with hi. These approximations simplify Eq. ()4.21|) by 
eliminating the implicit dependence of the jet and soft weights on the 
full final state. We may now treat hi as an independent vector. 

3. In the leading regions, particles that make up each final-state jet are 
associated with states Nj^, while A^^ consists of soft particles only. In 
the momentum conservation delta function, we can neglect the four- 
momenta of lines in Ng, whose energies all vanish as e,e ^ 0: 

S\pi-piNj,) -piNj,) -piNg)) ^ S\pj-pj^-pj^). (4.24) 

4. Because the cross section is a smooth function of the jet energies and 
directions, we may also neglect the masses of the jets within the momen- 
tum conservation delta function, as in Eq. ()4.10j) . In this approximation. 
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we derive in the cm., 



S\pi-Pj,-Pjj ^ SiV~s-co{NjJ-co{Nj,))5i\pj,\-\pj,\) ^ 



\PJi\ 



X + 712) 



s \ 2 ' '7 V 2 
x6\ni + n2). (4.25) 

Our jets are now back-to-back: 

n2 -fii . (4.26) 

Implementing these replacements and approximations for a < 1, we rewrite 
the cross section Eq. 1)4.211) as 



da{e,e,s,a) dao , . ^ / ^- - ^ 
dede-dn, = d^^^^'^^'^^ I ds.S{e,e.,a,^^) 
2 



X 



Yl / dsj^ Jc{£j,, a, fi)S{e- sj^ - sj^ - e^) , (4.27) 

c=l 



with (as above) H = 1 + 0{as)- Referring to the notation of Eqs. ()4.2ip and 
()4.22p . the functions S and are: 

5(£,£„a,/i) = 5^ 5(iV„;x)5(£-/(iV,))5(£,-/(iV„a)) (4.28) 

Ns 

UeJ^,a,^^) = ^(2vr)^ j;(iV^^, /i) 5(e^^ - /,(iV^^, a)) 5 - c.(iVj^ 



z2/ 



x6'{ni±n{NjJ), 

(4.29) 
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with the plus sign in the angular delta function for jet 2, and the minus for 
jet 1. The weight functions for the jets are given by Eq. ()4.22|1 and induce 
dependence on the parameter a. We have introduced the factorization scale 
/i, which we set equal to the renormalization scale. 

We note that we must construct the soft functions S{Ns,^) to cancel the 
contributions of final-state particles from each of the Jc{Nj^, fi) to the weight e, 
as well as the contributions of the jet functions to e from soft radiation outside 
their respective regions Cl^. Similarly, the jet amplitudes must be constructed 
to include coUinear enhancements only in their respective jet directions. Ex- 
plicit constructions that satisfy these requirements will be specified in the 
following subsections. 

To disentangle the convolution in ()4.27j) . we take Laplace moments with 
respect to e: 



Here and below unbarred quantities are the transforms in e, and barred quan- 
tities denote untransformed functions. 



^0 

and similarly for the jet functions. 

In the following subsections, we give explicit constructions for the functions 
participating in the factorization formula ()4.27|1 . which satisfy the requirement 




(4.30) 




(4.31) 
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of infrared safety, and avoid double counting. An illustration of the cross 
section factorized into these functions is shown in Fig. 14.41 As discussed above, 
non-global logarithms will emerge when ev becomes small enough. 




Figure 4.4: Factorized cross section ()4.27|1 after the application of Ward iden- 
tities. The vertical line denotes the final state cut. 

4.3.3 The short-distance function 

The power counting described in [75] shows that in Feynman gauge the sub- 
diagrams of Fig. 14.41 that contribute to H in Eq. ()4.27|1 at leading power in e 
and e are connected to each of the two jet subdiagrams by a single on-shell 
quark line, along with a possible set of on-shell, collinear gluon lines that carry 
scalar polarizations. The hard subdiagram is not connected directly to the soft 
sub diagram in any leading region. 

The couplings of the scalar-polarized gluons that connect the jets with 
short- distance subdiagrams may be simplified with the help of Ward identities 
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(see, e.g. [8]). At each order of perturbation theory, the coupUng of scalar- 
polarized gluons from either jet to the short-distance function is equivalent 
to their coupling to a path-ordered exponential of the gauge field, oriented in 
any direction that is not coUinear to the jet. Corrections are infrared safe, 
and can be absorbed into the short-distance function. Let h{pj^,ni, A) rep- 
resent the set of all short-distance contributions to diagrams that couple any 
number of scalar-polarized gluons to the jets, in the amplitude for the produc- 
tion of any final state. The argument A stands for the fields that create the 
scalar-polarized gluons linking the short-distance function to the jets. On a 
diagram-by-diagram basis, h depends on the momentum of each of the scalar- 
polarized gluons. After the sum over all diagrams, however, we can make the 
replacement: 

h{pj^, ni, ^(^'^)) ^ $® (0, -oo; 0) hipj^, h,,Q <f (0, -oo; 0) , (4.32) 

where /i2 is a short-distance function that depends only on the total momenta 
Pj^ and p j2 . It also depends on vectors that characterize the path-ordered 
exponentials $(0, — oo; 0): 

$fj(0, -oo; 0) = Pe-'sI-o.d^^c-Am(^ic) ^ (4 33) 

where the superscript (f) indicates that the vector potential takes values in 
representation f , in our case the representation of a quark or antiquark. These 
operators will be associated with gauge- invariant definitions of the jet functions 
below. To avoid spurious coUinear singularities, we choose the vectors 
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c = 1,2, off the light cone. In the full cross section ()4.3U|) the ^c-dependence 
cancels, of course. 

The dimensionless short- distance function H = \h2f in Eq. ()4.27p depends 
on ^/s and pj^ ■ ^c, but not on any variable that vanishes with e and e: 

H{pj^,^,,hi,fi) = Hl^,^^,h,,a,{fi)] , (4.34) 

where 

L^UVm- (4.35) 

Here we have observed that each diagram is independent of the overall scale 
of the eikonal vector 

4.3.4 The jet functions 

The jet functions and the soft functions in Eq. ()4.27|1 can be defined in terms 
of specific matrix elements, which absorb the relevant contributions to leading 
regions in the cross section, and which are infrared safe. Their perturbative 
expansions specify the functions S and J'c of Eq. (j4.29|) . We begin with our 
definition of the jet functions. 

The jet functions, which absorb enhancements collinear to the two outgoing 
particles produced in the primary hard scattering, can be defined in terms 
of matrix elements in a manner reminiscent of parton distribution or decay 
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functions [77]. To be specific, we consider the quark jet function: 



7 5^Tr[7'^(o|<|.g)t(o,-oo;0)g(0) N, 

X (iV^^|g(0)<l>g)(0,-oo;0) 

>i5\n,-h{Nj:)), 



(4.36) 



where Mc is the number of colors, and where he denotes the direction of 
the momentum of jet c, Eq. ()4.29|) . with h2 = —hi. q is the quark field, 
$^^^(0, — cxd; 0) a path-ordered exponential in the notation of ()4.33p . and the 
trace is taken over color and Dirac indices. We have chosen the normalization 
so that the jet functions J'^ in ()4.36|) are dimensionless and begin at lowest 
order with 



J'^(°)(ej„a,/i)=/5,^5(e,J, 



(4.37) 



with Pj! the lightlike velocities corresponding to the jet momenta in Eq. ()4.1()|1 : 



(4.38) 



The scalar jet functions of Eq. ()4.29p are now obtained by projecting out the 
component of J'^^ in the jet direction: 



(4.39) 
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where /3i = /32, /32 = /?i are the hghthke vectors in the directions opposite to 
Pi and /?2, respectively. By construction, the Jc are hnear in /?c- 

To resum the jet functions in the variables ej^, it is convenient to reexpress 
the weight functions (j4.22|) in combinations of light-cone momentum compo- 
nents that are invariant under boosts in the X3 direction, 

fANj.,a) = ^ E {2pj^ktY~\ (4.40) 

Here we have used the relation ■\/i/2 = uj^, valid for both jets in the cm. At 
the same time, we make the identification, 

-J - ^(iVjj) - HNjJ) = \6' ipj^ - piNjJ) , (4.41) 

which again holds in the cm. frame. The spatial components of each pj^ are 
thus fixed. Given that we are at small ej^, the jet functions may be thought 
of as functions of the light-like jet momenta pj^ of Eq. ()4.1()j) and of ej^. 
Because the vector jet function is constructed to be dimensionless, J'^ in Eq. 
()4.36|) is proportional to /5c rather than pj^. Otherwise, it is free of explicit 
/3c-dependence. 

The jet functions can now be written in terms of boost-invariant arguments. 



116 



homogeneous of degree zero in ^c- 

2e/?c-ec .(2) / pj.-^c ^ VI f \ 



1-a 



l~a 



(4.42) 



where J^^^ and J^^^ are independent functions, and where we have suppressed 
possible dependence on ^c,±- For jet c, the weight ej^ is fixed by 5{£j^ — 
fc{Nj^,a)), where on the right-hand side of the expression for the weight 
()4.40|) . the sum over each particle's momentum involves the overall factor 
{2p^^/y/sy^"'. After integration over final states at fixed ej^, the jet can 
thus depend on the vector pj^. At the same time, it is easy to see from 
the definition of the weight that pj^ can only appear in the combination 
(l/£j^v^)i/(i-'') {2pyy/^). This vector can combine with to form an in- 
variant, and all ^c-dependence comes about in this way. 

Expression ()4.42|1 can be further simplified by noting that 

2Pc-^cPc-^c=ec+ec,± ■ (4.43) 

Choosing ± = 0, we find a single combination. 



J, {ej^, a, /x) = Je I ej^ ^ (Cc)'"'' , a, a.(/x) ) , (4.44) 
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where, in the notation of Eq. (j4.42p . Jc = Jc + Jc , and we have defined 



In these terms, the Laplace moments of the jet function inherit dependence 
on the moment variable u through 



where the unbarred and barred quantities denote transformed and untrans- 
formed functions, respectively. We have constructed the jet functions to be 
independent of e, since the radiation into Q is at wide angles from the jet axes 
and can therefore be completely factored from the collinear radiation. This 
radiation at wide angles is contained in the soft function, which will be defined 
below in a manner that avoids double counting in the cross section. 

4.3.5 The soft function 

Given the definitions for the jet functions in the previous subsection, and the 
factorization ()4.27|) . we may in principle calculate the soft function S order by 
order in perturbation theory. We can derive a more explicit definition of the 
soft function, however, by relating it to an eikonal analog of Eq. ()4.27|) . 

As reviewed in Refs. [66,8], soft radiation at wide angles from the jets 
decouples from the collinear lines within the jet. As a result, to compute 
amplitudes for wide-angle radiation, the jets may be replaced by nonabelian 




2pj. ■ ec 



(4.45) 




(4.46) 
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phases, or Wilson lines. We therefore construct a dimensionless quantity, a^''''^\ 
in which gluons are radiated by path-ordered exponentials $, which mimic the 
color flow of outgoing quarks. 



^f{oo, 0; x) = Pe~^sIo-'dxp^-Ai^H^P.+-)^ (4.47) 



with Pc a light-like velocity in either of the jet directions. For the two-jet cross 
section at measured e and e^ik, we deflne 



iVcik 

0^ 



X (iV,,,|$|;^^)(oo,0;0)<l>g)(oo,0;0) 

= 6{e)6{e^,^) + 0{as). (4.48) 

The sum is over all final states N^^, in the eikonal cross section. The renormal- 
ization scale in this cross section, which will also serve as a factorization scale, 
is denoted jj,. Here the event shape function e^-^^^ is defined by /(A^eik, a) as in 
Eqs. ()4.3|) and ()4.4|) . distinguishing between the hemispheres around the jets. 
As usual, Afc is the number of colors, and a trace over color is understood. 

The eikonal cross section ()4.48p models the soft radiation away from the 
jets, including the radiation into Q, accurately. It also contains enhancements 
for configurations collinear to the jets, which, however, are already taken into 
account by the partonic jet functions in ()4.27|) . Indeed, ()4.48|) does not repro- 
duce the partonic cross section accurately for collinear radiation. It is also easy 
to verify at lowest order that even at fixed e^^j^ the eikonal cross section ()4.48p 
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is ultraviolet divergent in dimensional regularization, unless we also impose a 
cutoff on the energy of real gluon emission collinear to (3i or f]2- 

The construction of the soft function S from (t('="") is nevertheless possi- 
ble because the eikonal cross section (|4.48p factorizes in the same manner as 
the cross section itself, into eikonal jet functions and a soft function. The 
essential point [65] is that the soft function in the factorized eikonal cross 
section is the same as in the original cross section ()4.27|1 . The eikonal jets or- 
ganize all collinear enhancements in ()4.48p . including the spurious ultraviolet 
divergences. These eikonal jet functions are defined analogously to their par- 
tonic counterparts, Eq. ()4.H(ij) . but now with ordered exponentials replacing 
the quark fields. 



J(-)(e„a,/i) = ^ J] (o|$J)t(o,-oo;0)$j:ft(^^0;0) 

^ (oik) 



^(cik) 



<l>g)(oo,0;0)$f)(0,-oo;0) 

= 6{e,) + 0{a,), (4.49) 

where fc is a quark or antiquark, and where the trace over color is understood. 
The weight functions are given as above, by Eq. ()4.22|) . with the sum over 
particles in all directions. 

In terms of the eikonal jets, the eikonal cross section ()4.48jl factorizes as 



0-^"'') {e, e^ik, a, /i) 



c=l ^ 

X 6 -Cs-ei- 62) , (4.50) 
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where we pick the factorization scale equal to the renormahzation scale /i. As 
for the full cross section, the convolution in ()4.5()|1 is simplified by a Laplace 
transformation ()4.46|) with respect to e^n,, which allows us to solve for the soft 
function as 

S (e, u, a, /i) = ^^"^ = 5(6) + . (4.51) 

c=l 

In this ratio, collinear logarithms in u and the unphysical ultraviolet diver- 
gences and their associated cutoff dependence cancel between the eikonal cross 
section and the eikonal jets, leaving a soft function that is entirely free of 
collinear enhancements. The soft function retains //-dependence through soft 
emission, which is also restricted by the weight function e. In addition, be- 
cause soft radiation within the eikonal jets can be factored from its collinear 
radiation, just as in the partonic jets, all logarithms in u associated with wide- 
angle radiation are identical between the partonic and eikonal jets, and factor 
from logarithmic corrections associated with collinear radiation in both cases. 
As a result, the inverse eikonal jet functions cancel contributions from the 
wide-angle soft radiation of the partonic jets in the transformed cross section 
(ICTll . 

Given the definition of the energy flow weight function /, Eq. ()4.5p . the 
soft function is not boost invariant. In addition, because it is free of collinear 
logs, it can have at most a single logarithm per loop. Its dependence on e 
is therefore only through ratios of the dimensional quantities Ey/s with the 
renormahzation (factorization) scale. 
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As in the case of the partonic jets, Eq. ()4.46|) . we need to identify the vari- 
able through which u appears in the soft function. We note that dependence 
on the velocity vectors (3c and the factorization vectors S,c must be scale invari- 
ant in each, since they arise only from eikonal lines and vertices. The eikonal 
jet functions cannot depend explicitly on the scale-less, lightlike eikonal veloc- 
ities /3c, and independent of the ^c- Dependence on the factorization 
vectors enters only through the weight functions, ()4.4()|1 for the eikonal jets, 
in a manner analogous to the case of the partonic jets. This results in a de- 
pendence on (Cc)^"'^, as above, with defined in Eq. ()4.45|1 . In summary, we 
may characterize the arguments of the soft function in transform space as 



4.4 Resummation 

We may summarize the results of the previous section by rewriting the trans- 
form of the factorized cross section (j4.3(jp in terms of the hard, jet and soft 
functions identified above, which depend on the kinematic variables and the 
moment v according to Eqs. ()4.34|) . ()4.46|) and ()4.52|) respectively. 




(4.52) 
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The natural scale for the strong coupling in the short- distance function H is 
-y/s/2. Setting fi = \/s/2, however, introduces large logarithms of e in the 
soft function and large logarithms of u in both the soft and jet functions. The 
purpose of this section is to control these logarithms by the identification and 
solution of renormalization group and evolution equations. 

The information necessary to perform the resummations is already present 
in the factorization ()4.53|1 . The cross section itself is independent of the fac- 
torization scale 



dfi dedrii 



and of the choice of the eikonal directions, ^c, used in the factorization, 

d da {e, u, s, a) 



. (4.55) 



The arguments of this section closely follow the analysis of Ref. [80]. We will 
see that the dependence of jet and soft functions on the parameter a that 
characterizes the event shapes (3) is reflected in the resummed correlations, 
so that the relationship between correlations with different values of a is both 
calculable and nontrivial. 



4.4.1 Energy flow 

As a first step, we use the renormalization group equation ()4.54|) to organize 
dependence on the energy flow variable e. Applying Eq. ()4.54j] to the factorized 
correlation ()4.53|1 . we derive the following consistency conditions, which are 
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themselves renormalization group equations: 



/i— In S ,eiy, — {Q °, a, = -7s (4.56) 

In J, (pjAl^ }ll (^Ji-, a, as{ii) \ = -7j. , (4.57) 

ail \ jjL {JLU 

/i In I — , ^-^^ ^' ,ni,a,(/i) I = 7s («s(/i)) + Y] 7Jc (".(z^)) • 

(4.58) 



The anomalous dimensions 7^, d = s, Jc can depend only on variables held 
in common between at least two of the functions. Because each function is 
infrared safe, while ultraviolet divergences are present only in virtual diagrams, 
the anomalous dimensions cannot depend on the parameters u, e or a. This 
leaves as arguments of the 7^ only the coupling as(/i), which we exhibit, and 
(c, which we suppress for now. 

Solving Eqs. (p3K|l and (p37jl we find 



S {^^,eu, — {Q^ " ,a,as{fi)] = S ( ^^,eu, — {Q^ " ,a,as{fio) 

- J f >{"4A)) 

X e "0 ^ (4.59) 

T I PJc ■ , . xl_a / N \ J I PJc ■ . . xl-a /- N 

Jc , — (Cc) ,a,«s(/i) = -'c — z ,- — (Cc) ,a,a4/io) 

- / t7J.(a.(A)) 

X e ''u , (4.60) 



for the soft and jet functions. As suggested above, we will eventually pick 
H ~ y/s to avoid large logs in H. Using these expressions in Eq. ()4.53jl we 
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can avoid logarithms of e or z/ in the soft function, by evolving from /io = e^/s 
to the factorization scale /i ~ y/s. No choice of /io, however, controls all 
logarithms of u in the jet functions. Leaving ftQ free, we find for the cross 
section ()4.5H|1 the intermediate result 



da {e, u, s, a) dao ( y/s pj^ ■ ic . , . 

H — , , ni, asifi) 



de dhi dhi \ /i fi 

X S (l, eu, iCc)^'", a, asie^/s)) expl- J 



^IsiasiX))} (4.61) 



X Jc I ^"^1 ^ (Cc)^ " , a, a,(/io) I exp J - /" ^7j^ ("^(A)) 



We have avoided introducing logarithms of e into the jet functions, which 
originally only depend on by evolving the soft and the jet functions inde- 
pendently. The choice of /xq = S\/s or \fslv for the soft function is to some 
extent a matter of convenience, since the two choices differ by logarithms of ev. 
In general, if we choose /iq = y/s/i'-, logarithms of ev will appear multiplied by 
coefficients that reflect the size of region Vt. An example is Eq. ()4.15p above. 
When VL has a small angular size, /io = \/s/v is generally the more natural 
choice, since then logarithms in ev will enter with small weights. In contrast, 
when VL grows to cover most angular directions, as in the study of rapidity 
gaps [82], it is more natural to choose /io = e^/s. 
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4.4.2 Event shape transform 

The remaining unorganized "large" logarithms in ()4.fiH) . are in the jet func- 
tions, which we will resum by using the consistency equation ()4.55|) . The 
requirement that the cross section be independent of pj^ ■ implies that the 
jet, soft and hard functions obey equations analogous to (j4.56|) - (j4.58p . again 
in terms of the variables that they hold in common [80]. The same results may 
be derived following the method of Collins and Soper [72], by defining the jets 
in an axial gauge, and then studying their variations under boosts. 

For our purposes, only the equation satisfied by the jet functions [72,80] is 
necessary. 



The functions Kc and Gc compensate the ^c-dependence of the soft and hard 
functions, respectively, which determines the kinematic variables upon which 
they may depend. In particular, notice the combination of u- and ^c-dependence 
required by the arguments of the jet function, Eq. ()4.4fj|l . 

Since the definition of our jet functions ()4.36p is gauge invariant, we can 
derive the kernels and Gc by an explicit computation of d Jc/ din {p,j^ ■ 
in any gauge. The multiplicative renormalizability of the jet function, Eq. 
fl4.57|) . with an anomalous dimension that is independent of pj^ ■ ensures 
that the right-hand side of Eq. ()4.fi2j) is a renormalization-group invariant. 
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Thus, Kc + Gc are renormalized additively, and satisfy [72] 



ail ynu 

li^cJ^^^asiii)] = 7K.(«s(/i)). (4.63) 



Since Gc and hence jKc^ i^^iy be computed from virtual diagrams, they do not 
depend on a, and •jkc is the universal Sudakov anomalous dimension [72,83]. 

With the help of these evolution equations, the terms Kc and Gc in Eq. 
()4.62|) can be reexpressed as [84] 

- / ^7... («. (A')) 

— A;(ci,a,a, (A')) , 
(4.64) 

where in the second equality we have shifted the argument of the running 









, a 
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coupling in Kc^ and have introduced the notation 



5^ (ci,C2,a,a5 (/i)) = 
2A'Jci,a,a3 (/i)) = 



-Kc ( — , a, as (/i) ) - Gc ( — , (/i) 



(4.65) 



The primes on the functions A'^ and B'^ are to distinguish these anomalous 
dimensions from their somewhat more famihar versions given below. 
The solution to Eq. ()4.(i2|l with = /io is 



X exp < 



dX 
T 



v^/(2Co) 



Jr ( 777^, ^ (Co)^ " , a, as(/io) 

C2 A 



/(iA' 
— A'^ (ci,a,a, (A')) 



^ 1/(2 A)!-"! 



(4.66) 



where we evolve from \fs/{2 Co) to pj^ ■ = y/s/ (2 Cc) (see Eq. ()4.45|l ) with 



z/\i/(2-a) 



(4.67) 



After combining Eqs. (jCTH) and (jOHIl . the choice jlo = A/i/(2Co) = 
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allows us to control all large logarithms in the jet functions simultaneously: 



X exp 



X exp 



v/5/(2Co) 

- / 

V5/(2Co) 



dX 

T 



C2 A 

— ^^(ci,a,Q;s (A')) 



^1/(2 A)i-» 



(4. 



As observed above, we treat a as a fixed parameter, with \a\ small compared 
to In (l/e) and Inu. 
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4.4.3 The resummed correlation 



Using Eq. ()4.fi8j) in ()4.f)l|l . and setting fi = y/s/2, we find a fully resummed 
form for the correlation, 



da {e, u, s, a) dao ^ ( 2pj^-^, 



de dfii 



dhi 



(v/i/2 
- J y7s(as(A)) 



n 



c=l 



2 Co 



v^/2 



exp 



d\ 
T 



X exp < 



dX 
T 



v^/(2Co) 



v^/{2Co) 
C2 A 



/dX' 
— (ci,a,a, (A')) 



^i^(2A)l-<' 



> . 



(4.69) 



Alternatively, we can combine all jet-related exponents in Eq. ()4.69p in 
the correlation. As we will verify below in Section 14.5.21 the cross section is 
independent of the choice of ^c- As a result, we can choose 



(4.70) 
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This choice allows us to combine and B'^ in Eq. ()4.69p . 



da (e, z/, s, a) 
de dhi 



5* (l, eu, 1, a, ^^(eyi)) exp 



dX 



7s MX)) 1,1, a 



c=l 



a.. 



2 Co 



X exp 



y [7Jc + K (ci, C2, a, a, (caA)) 



v^/(2Co) 
C2 A 



/dX' 
— A',(ci,a,a,(A')) 



gl-a/2 
1^(2 A)!-" 



(4.7 



with Co given by Eq. ()4.(j7|l . 

In Eqs. fl4.69p and ()4.7ip . the energy flow e appears at the level of one 
logarithm per loop, in S, in the first exponent. Leading logarithms of e are 
therefore resummed by knowledge of 7!^'' , the one-loop soft anomalous dimen- 
sion, where we employ the standard notation. 



00 



n=0 



(4.72) 



for any expansion in as. At the same time, u appears in up to two logarithms 
per loop, characteristic of conventional Sudakov resummation. To control 
i^-dependence at the same level as e-dependence, it is natural to work to next- 
to- leading logarithm in u, by which we mean the level a J In'^ u in the exponent. 
As usual, this requires one loop in B'^ and 7j^, and two loops in the Sudakov 
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anomalous dimension A'^, Eq. ()4.65p . These functions are straightforward to 
calculate from their definitions given in the previous sections. Only the soft 
function S in Eqs. fl4.69p and ()4.7H) contains information on the geometry of 
Q. The exponents are partially pro cess- dependent, but geometry- independent. 
In Section 1^31 we will derive explicit expressions for these quantities, suitable 
for resummation to leading logarithm in e and next-to-leading logarithm in u. 

4AA The inclusive event shape 

It is also of interest to consider the cross section for e"'"e~-annihilation into 
two jets without fixing the energy of radiation into Q, but with the final state 
radiation into all of phase space weighted according to Eq. ()4.4j) , schematically 

e+ + ^ Ji{pj,j\) + J2{pj,Jn,) , (4.73) 

where f^i and ^2 cover the entire sphere. This cross section can be factorized 
and resummed in a completely analogous manner. The final state is a convo- 
lution in the contributions of the jet and soft functions to e as in Eq. ()4.27|) . 
but with no separate restriction on energy flow into Q. All particles contribute 
to the event shape. We obtain an expression very analogous to Eq. ()4.69|) for 
this inclusive event shape in transform space, which can be written in terms 
of the same jet functions as before, and a new function S*™"^^ for soft radiation 
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as: 



H I ^,ni,a^(Vs/2) 



dni 



dhi 



- I y7.(«.(A)) 



n jcUa 



c=l 



2 Co 



yCs/2 



exp 



A 



X exp < 



dX 
T 



v^/(2Co) 



v^/(2Co) 

C2 A 



/(iA' 
— A'^ (ci,a,as (A')) 



> . 



^ 1/(2 A)l-'i 



(4.74) 



Here the soft function S'™'^' = 1 + (9(as). The double-logarithmic dependence 
of the shape transform is identical to our resummed correlation, Eq. ()4.69|) . 
We will show below, in Sec. 14. 5. 3^ that Eq. ()4.74p coincides at NLL with the 
known result for the thrust [57] when we choose a = 0. 



4.5 Results at NLL 

4.5.1 Lowest order functions and anomalous dimensions 

In this section, we describe the low-order calculations and results that provide 
explicit expressions for the resummed shape/flow correlations and inclusive 
event shape distributions at next-to-leading logarithm in u and leading log- 
arithm in e (we refer to this level collectively as NLL below). We go on to 
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verify that for the case a = we rederive the known result for the resummed 
thrust at NLL, and we exhibit the expressions for the correlation that we will 
evaluate in Sec. 14.61 



The soft function 

The one-loop soft anomalous dimension is readily calculated in Feynman gauge 
from the combination of virtual diagrams in Eq. (IOH|l . and J^"'"^), Eq. 

()4.49|1 . in Eq. ()4.5H1 . The calculation and the result are equivalent to those of 
Ref. [65], where the soft function was formulated in axial gauge, 



7, 



2 / 

^In - In 

c=l ^ 



(4.75) 



The first, ^c-dependent logarithmic term is associated with the eikonal jets, 
while the second is a finite remainder from the combination of a^^^^^ and J^^^^^ 
in ()4.51|) . Whenever C,c,± =0, the logarithmic terms cancel identically, leaving 
only the final term, which comes from the C,c eikonal self-energy diagrams in 
the eikonal jet functions. 

The soft function is normalized to S^^\e) = 5{e) as can be seen from ()4.5ip . 
For non-zero e, da/de is given at lowest order by 

is ^0,n) = C- [ dPS, ^ /yf^ , (4.76) 



where PS2 denotes the two-dimensional angular phase space to be integrated 
over region fl, and k = k/uk. We emphasize again that the soft function 
contains the only geometry-dependence of the cross section. Also, S^^^ for 
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£ 7^ is independent of v and a. 

As an example, consider a cone with opening angle 25, centered at angle 
a from jet 1. In this case, the lowest-order soft function is given by 

5(1) [e ^ 0, a, b) = Cf- In ( \~ • (4.77) 

e \cos^ a — cos^ o J 

Similarly, we may choose as a ring extending angle 5i to the right and 82 to 
the left of the plane perpendicular to the jet directions in the center-of-mass. 
In this case, we obtain 

S") (e ^ 0. S,) = C.i In ( il±£i!iMii±£HiM) = C.? A, . (4.78) 

e V (1 ~ sm Oi) (1 — sm ^2) / e 

with A?7 the rapidity spanned by the ring. For a ring centered around the 
center-of-mass (5i = 82 = 8) the angular integral reduces to the form that we 
encountered in the example of Sec. I4.2.3[ and that we will use in our numerical 
examples of Sec. 14. 6| with Ar] given by Eq. ()4.14|1 . 

The jet functions 

Recall from Eq. ()4.39j) that the lowest-order jet function is given by J(°) = 1. 
The anomalous dimensions of the jet functions are found to be 

lf^=-lcF, (4.79) 

the same for each of the jets. The jet anomalous dimensions are process- 
independent, but of course flavor-dependent. The same anomalous dimensions 
for final-state quark jets appear in three- and higher-jet cross sections. 
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The ii'-G-decomposition 

The anomalous dimension for the ii'-G- decomposition is, as noted above, the 
Sudakov anomalous dimension. 



7g 



KCf, 



(4.80) 
(4.81) 



also independent of the jet-direction. The well-known coefficient K (not to be 
confused with the functions K^) is given by [85] 



\18 6 J ^ 9 



(4.82) 



with the normalization Tp = 1/2 and Nf the number of quark flavors. 

Kc and Gc, the functions that describe the evolution of the jet functions 
in Eq. ()4.(j2|l . are given at one loop by 



. , , Ql-a/2 / „ ^ a-1 



1) / PJc ■ Co 



-Cf In [e 
/ 

-Cf In e 



2(l-a) 



(4.83) 



Evolving them to the values of fi with which they appear in the functions A'^ 
and B'^, Eq. (|4.65p . they become 



Cl 



-Cf In {e^^^'^'-^^cl) 



C2 



-C^ln( 

^2 



(4.84) 
(4.85) 
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Recall that Gc is computed from virtual diagrams only, and thus does not 
depend on the weight function. It therefore agrees with the result found in [72]. 
The soft-gluon contribution, Kc-, which involves real gluon diagrams, does 
depend on the cross section being resummed. 

With the definitions ()4.65p of and B'^ we obtain 



(ci,a) 
B'}'^ (ci,C2,a) 



2 



2Cf In e 



3 7E 



2 

C2 



I) 



(4.86) 
(4.87) 

(4.88) 



Here Pq is the one-loop coefficient of the QCD beta-function, /?o = | (HA/c — ATpNf) 
The hard scattering, and the Born cross section 

At NLL only the lowest-order hard scattering function contributes, which is 
normalized to 

ij(°)(a,(v^/2)) = l. (4.89) 

At this order the hard function is independent of the eikonal vectors C,c, al- 
though it acquires ^c-dependence at higher order through the factorization 
described in Sec. 14.3.31 For completeness, we also give the electromagnetic 
Born cross section at fixed polar and azimuthal angle: 



dao 
dhi 



(4.90) 
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where 6 is the cm. polar angle of ni, eQf is the charge of quark flavor f, and 
Q^em = e^/(47r) is the fine structure constant. 

4.5.2 Checking the ^c-dependence 

It is instructive to verify how dependence on the eikonal vectors cancels in 
the exponents of the resummed cross section ()4.69|) at the accuracy at which 
we work, single logarithms of e, and single and double logarithms of u. In 
these exponents, ,^c-dependence enters only through the combinations {Pc ■ ic) 
and {pj^ ■ ic)- 

Let us introduce the following notation for the exponents in Eq. ()4.(j9|l . to 
which we will return below: 

/ y%(«.(A))-XJ / y7j.(«.(A)), (4.91^ 



E, = - 



9 P./c-?c 

y^/(2Co) 



C2 A 

/d\' 
— A'Jci,a,a, (A')) 



1 1/(2 A)l-" 



(4.92) 



At NLL, explicit dependence is found only in 7^, Eq. (j4.75j) . for Ei, and in 
the upper limit of the A integral of E2. We then find that 

v/i/2 

^ {El + E2) = 2Cf / ^ -2Cp — + NNLL . 



(4. 
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Here the second term stems entirely from A'^^\ Eq. ()4.86p : other contribu- 
tions of E2 are subleading. The ^c-dependence in the exponents begins only 
at the level that we do not resum, at ln(l/£:z/), which is compensated by 
corrections in S'(ez/, a^). The remaining contributions are of NNLL order, that 
is, proportional to ag{y/s) In*^^^ (^i^ Pc ■ icj, as may be verified by expanding 
the running couplings. Thus, as required by the factorization procedure, the 
relevant ^c-dependence cancels between the resummed soft and jet functions, 
which give rise to the first and second integrals, respectively, in Eq. ()4.93p . 

4.5.3 The inclusive event shape at NLL 

We can simplify the differential event shape, Eq. ()4.74j) . by absorbing the soft 
anomalous dimension 7^ into the remaining terms. We will find a form that can 
be compared directly to the classic NLL resummation for the thrust (a = 0). 
This is done by rewriting the integral over the soft anomalous dimension as 





^l-a/2/[j,(2A)l-a] 



dX' 



f3i9{X')) — ^Aa,iX')) . 

(4.94) 
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In the first equality we split the A integral so that the limits of the first term 
match those of the B'^ integral of Eq. ()4.74|1 . In the second equality we have 
changed variables in the second term according to 



so that the limits of the second integral also match. In the third equality of 
Eq. ()4.94|) . we have reexpressed the running coupling at the old scale A in 
terms of the new scale. This is a generalization of the procedure of Ref. [86], 
applied originally to the threshold-resummed Drell-Yan cross section [87]. 

Using Eq. ((Olj), and identifying pj^ ■ 4 with y/s/2 (Eq. dUTDI)) in the 
inclusive event shape distribution, Eq. ()4.74|) . we can rewrite this distribution 
at NLL as 




(4.95) 



da''''-^ {u, s, a) 



dao 



dfii 



dhi 




(4.96) 
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where we have rearranged the contribution of 7s as: 

1 d 
Ac(ci,a,as (/i)) = A'^ici,a,asifi))--il-a)f3{g{iJ,))—-fsiasilJ')), 

Be (ci, C2, a, as (/i)) = 7j, + ^ ^) ("^(/^)) + K (ci, Ca, a, (/i)) • 

(4.97) 

Next, we replace the lower limit of the A'-integral by an explicit 6'-function. 
Then we exchange orders of integration, and change variables in the term con- 
taining A from the dimensionful variable A to the dimensionless combination 



« = ^ . (4.98) 



We find 
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dX' 
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dfii dhi 1 J A 

Vi/[2(!//2)l/{2-a 
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Be (ci,C2,a,as (A)) 



n-P|-2/ ^ - ^ (c,V-^ -ij A. (c„a,a. (AO) 

(4.99) 



Here, the ^-function vanishes for small A', and the remaining effects of re- 
placing the lower boundary of the A' integral by are next-to-next-to-leading 
logarithmic. 

A further change of variables allows us to write the NLL resummed event 
shapes in a form familiar from the NLL resummed thrust. In the first line of 
Eq. ()4.99|) . we replace A^ — > us /A. In the second line we relabel A' a/^. 
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and exchange orders of integration. Finally, choosing 



Cl 
C2 



-IE 



(4.100) 



we find at NLL 



(z/, s, a) 







dni 



dhi 

2 



X 



+ -B, {as{us/A)) (e 



-u(i//2)2/(2-<!)e-7B 



and reproduce the well-known coefficients 



5(1) 



Cf, 

-CpK, 
2 



independent of a. In Eq. ()4.101|) . we have made use of the relation 



(4.101) 



(4.102) 
(4.103) 
(4.104) 



(4.105) 



which is valid at NLL in the logarithmic integrals. With these choices, when 
a = we reproduce the NLL resummed thrust cross section [57]. 
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The choices of the q in Eq. ()4.100j) cancel all purely soft NLL compo- 
nents (7s and Kc). The remaining double logarithms stem from simultane- 
ously soft and collinear radiation, and single logarithms arise from coUinear 
configurations only. At NLL, the cross section is determined by the anomalous 
dimension A^, which is the coefficient of the singular 1/[1 — term in the 
nonsinglet evolution kernel [88], and the quark anomalous dimension. All radi- 
ation in dijet events thus appears to be emitted coherently by the two jets [57]. 
This, however, is not necessarily true beyond next-to-leading logarithmic ac- 
curacy for dijets, and is certainly not the case for multijet events [65]. Similar 
considerations apply to the resummed correlation, Eq. ()4.69|1 . 

4.5.4 Closed expressions 

Given the explicit results above, the integrals in the exponents of the resummed 
correlation, Eq. ()4.69|) . may be easily performed in closed form. We give the 
analytic results for the exponents of Eq. ()4.69|1 . as defined in Eqs. ()4.9H) and 
(11321) • As in Eq. (jTZOl), we identify pj^ ■ ic with y/s/2. 




6Cp 



(4.106) 
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1 -K^F I \ 



X 



(4.107) 
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with 



Ki{a) 
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(4.108) 



(4.109) 



We have used the two-loop running couphng, when appropriate, to derive Eqs. 
fl4.106p - ()4.109|) . The results are expressed in terms of the one- loop running 
coupling 



27r 



/?o In 



Aqcd 



(4.110) 



and the first two coefficients in the expansion of the QCD beta-function, /3o 
and 

34 /20 \ 

iix = -C\-[-CA^^CA TpNf. (4.111) 



Combining the expressions for the exponents, Eqs. ()4.10(j|l and ()4.107|1 . for the 
Born cross section, Eq. ()4.90|) . and for the soft function, Eq. ()4.76|) . in Eq. 
()4.69|) . the complete differential cross section, at LL in e and at NLL in z/, is 
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given by 
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vr '^T^ I3i ■ k 132 ■ k 
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4Cp , s 1 4Cp , > 



a, I ) / U.J ^ 



2 Co 



X I 1 . (4.112) 



as 



2 Co 



These are the expressions that we will evaluate in the next section. We note 
that this is not the only possible closed form for the resummed correlation 
at this level of accuracy. When a full next-to-leading order calculation for 
this set of event shapes is given, the matching procedure of [57] may be more 
convenient. 



4.6 Numerical Results 

Here we show some representative examples of numerical results for the cor- 
relation, Eq. ()4.112|1 . We pick the constants q as in Eq. ()4.1()()|1 . unless stated 
otherwise. The effect of different choices is nonleading, and is numerically 
small, as we will see below. In the following we choose the region f2 to be a 
ring between the jets, centered in their center-of-mass, with a width of At] = 2, 
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or equivalently, opening angle 5 ~ 50 degrees (see Eq. ()4.14|) ). The analogous 
cross section for a cone centered at 90 degrees from the jets (Eq. ()4.77j) ) has a 
similar behavior. In the following, the center-of-mass energy Q = ^/s is chosen 
to be 100 GeV. 




Figure 4.5: Differential cross section dil/dhi ^ normalized by the Born cross 
section, at Q = 100 GeV, as a function of e and a at fixed z/: a) u = 10, b) 
I' = 50. f2 is a ring (slice) centered around the jets, with a width of At] = 2. 



Fig. 14.51 shows the dependence of the differential cross section ()4.69|) . mul- 
tiplied by £ and normalized by the Born cross section, ^^^^y^^^, on the mea- 
sured energy e and on the parameter a, at fixed u. In Fig. 14.51 a), we plot 
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^ dcTa/^dhi ^ ~ ^^■> -^^S- 14.51 b) for V = 50. As z/ increases, the radiation 
into the complementary region is more restricted, as illustrated by the com- 
parison of Figs. 14.51 a) and b). Similarly, as a approaches 1, the cross section 
falls, because the jets are restricted to be very narrow. On the other hand, as 
a assumes more and more negative values at fixed e, the correlations ()4.69p 
approach a constant value. For a large and negative, however, non-global de- 
pendence on Ine and \a\ will emerge from higher order corrections in the soft 
function, which we do not include in Eq. ()4.112j) . 




Figure 4.6: Differential cross section ^ dXldtii ^o^^^^-^i^^d by the Born cross 
section, at Q = 100 GeV, as a function of a at fixed z/ = 20 and e = 0.05. Q is 
chosen as in Fig. 14 .51 Solid line: ci = e"'''^, C2 = 2, as in Eq. ()4.100p . dashed 
line: Ci = C2 = 1, dotted line: Ci = C2 = 2. 



In Fig. 14.61 we investigate the sensitivity of the resummed correlation, Eq. 
fl4.112j) . to our choice of the constants q. The effect of these constants is 
of next-to-next-to-leading logarithmic order in the event shape. We plot the 
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differential cross section dio/d^i ^ at Q = 100 GeV, for fixed e = 0.05 and 
= 20, as a function of a. Tlie effects of changes in the Cj are of the order of 
a few percent for moderate values of a. 

Finally, we illustrate the sensitivity of these results to the flavor of the 
primary partons. For this purpose we study the corresponding ratio of the 
shape/flow correlation to the cross section for gluon jets produced by a hy- 
pothetical color singlet source. Fig. 14.71 displays the ratio of the differential 
cross section da'^{e,a)/{dedni), Eq. ()4.112j) . normalized by the lowest-order 
cross section, to the analogous quantity with gluons as primary partons in 
the outgoing jets, again at Q = 100 GeV. This ratio is multiplied by Ca/Cf 
in the figure to compensate for the difference in the normalizations of the 
lowest-order soft functions. Gluon jets have wider angular extent, and hence 
are suppressed relative to quark jets with increasing v or a, as can be seen by 
comparing Figs. 14.71 a) and b). Fig. 14.71 a) shows the ratio at u = 10, and Fig. 
14. 71 b) at z/ = 50. These results suggest sensitivity to the more complex color 
and flavor flow characteristic of hadronic scattering [65,66]. 

4.7 Summary and Outlook 

We have introduced a general class of inclusive event shapes in e"'"e~ dijet 
events which reduce to the thrust and the jet broadening distributions as 
special cases. We have derived analytic expressions in transform space, and 
have shown the equivalence of our formalism at NLL with the well-known 
result for the thrust [57]. Separate studies of this class of event shapes in the 
untransformed space, at higher orders, and for nonperturbative effects [60] are 
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. 02"^ ■ ^ 

Figure 4.7: Ratios of differential cross sections for quark to gluon jets 
Ca. ^ = ^ function of e and a 

at fixed u: a) u = 10, h) u = 50. as in Fig. 14.51 Ci and C2 as in Eq. ()4.100p . 

certainly of interest. We reserve these studies for future work. 

We have introduced a set of correlations of interjet energy flow for the 
general class of event shapes, and have shown that for these quantities it is 
possible to control the influence of secondary radiation and nonglobal loga- 
rithms. These correlations are sensitive mainly to radiation emitted directly 
from the primary hard scattering, through transforms in the weight functions 
that suppress secondary, or non-global, radiation. We have presented analytic 



149 



and numerical studies of these shape/flow correlations at leading logarithmic 
order in the flow variable and at next-to-leading-logarithmic order in the event 
shape. The application of our formalism to multijet events and to scattering 
with initial state hadrons is certainly possible, and may shed light on the 
relationship between color and energy flow in hard scattering processes with 
non-trivial color exchange. 



150 



Appendix A 



A.l Power counting with contracted vertices 

In this appendix we will include the possibility of contracted vertices in the 
reduced diagram in Fig. 12.1b . These are associated with internal lines (col- 
lapsed to a point) which are off-shell by y/s. Our analysis closely follows [39] 
and [43]. 

If we go back to the argument that led us to Eq. ()2.15|) for the superficial 
degree of IR divergence for the soft part, we see that the same reasoning as in 
the case of elementary vertices applies to the case of contracted vertices since 
the result ()2.15|) has been obtained by means of dimensional counting. 

The analysis of contracted vertices connecting jet lines only is, however, 
more subtle. We have to demonstrate that the suppression factors correspond- 
ing to the contracted vertices are at least as great as the ones for the elementary 
vertices. The expression (I2.22j) tells us that we can restrict ourselves to the 
two and three point vertices. For these cases, we analyze the full two and 
three-point subdiagrams, by studying the tensor structures that are found 
after integration over their internal loop momenta. 
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Before we discuss all the possible structures, we state some results which 
will be essential for the upcoming analysis. The first one is the simple Dirac 
matrix identity 

^ ^ ^ = 2(a-6) ^-a^ ^. (A.l) 

The other two follow from Eqs. ()2.7|) and ()2.11|) for the gluon propagator 
in Coulomb gauge, and hold for any jet momenta scaling as Ia ^ I'a ^ 
A~, A^/^) coUinear to the momentum pA defined in Eq. fl2.2|) 

lTNapiU,v) = 0(A'/'v^), (A.2) 

for all components of (3. We now proceed to discuss the particular cases. 

Ghost self-energy: The most general covariant structure is, using p-p = p"^, 

1 

where /x is a scale introduced by a UV/IR regularization of Feynman diagrams 
and p is the momentum of an internal jet line. Strictly speaking, the covariants 
should be formed from the vectors p and rj, but since p has nonzero light-cone 
components, we can use Eq. ()2.8j) . to express t] in terms of p. The maximum 
degree of divergence for the ghost self-energy occurs when the internal lines 
become either parallel to the external momentum p or soft. The most general 
pinch singular surface consists of a subdiagram of collinear lines moving in 

a direction of the external ghost. This subdiagram can interact with itself 

^In the rest of this subsection we are concerned the momentum factors only, and we omit 
dependence on the color structure. 
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through the exchange of soft quanta. Power counting arguments similar to 
the ones given in Sec. 12.3.21 show, however, that there is no IR divergence for 
these pinch singular points. This shows that the dimensionless function / in 
Eq. ()A.3|) is IR finite. Hence the combination [tree level ghost propagator] - 
[ghost self-energy] - [tree level ghost propagator], [l/{p ■ p)]Il{p,p) [l/{p-p)], 
is suppressed at least as much as a single tree level ghost propagator, l/{p-p). 
Therefore the contracted two point ghost vertex within a jet subdiagram con- 
tributes at least the same suppression as a single tree level ghost propagator. 

Gluon self- energy: With external momentum p, its most general tensor 
decomposition has the form 

n^i.(p, P) = g^,u h + PpiPu f2 + Pf^Pu h + iPt^Pu + PfiPu) U ■ (A.4) 

As verified by explicit one-loop calculations in Refs. [49] and [50] the gluon 
self-energy in Coulomb gauge is not transverse. In Eq. ()A.4|) . the fi = 
fiij)^ / jji^ , p^ / jji^ , as{ij)) are dimensionless functions. Contracting 11^,^ with 
tree level gluon propagators, and using Eq. ()2.10|1 . the last two terms in Eq. 
flA.4|) drop out and the first and the second terms give at least one factor of p^ 
in the numerator, which cancels one of the denominator factors. Since 

the maximum degree of IR divergence for the gluon self-energy occurs when all 
the internal lines become either collinear to the external momentum p or soft, 
we can use the results of the power counting of Sec. 12.3.21 to demonstrate that 
the dimensionless functions are at worst logarithmically divergent. There- 
fore the combination: gluon jet line - 2 point gluon contracted vertex - gluon 
jet line, behaves the same way as a gluon jet line for the purpose of the jet 
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power counting. 

Fermion self-energy: In the massless fermion limit, the most general matrix 
structure of the fermion self-energy is 

=^gi+ pg2, (A.5) 

with dimensionless functions gi = gi{p'^/fi^, P^//^^? ^ = 1,2. When we 

sandwich the fermion self-energy between the tree level fermion denominators, 
the first term in Eq. ()A.5|1 behaves the same way as the tree level fermion 
propagator, modulo logarithmic enhancements due to the function gi. The 
second term, however, is absent from the fermion self-energy as was shown 
in Ref. [43] using the method of induction and Ward identities. The idea 
was to study a variation of the fermion self-energy by making an infinitesimal 
Lorentz boost on the external momentum. This implies a relationship between 
the (r + 1) and the r-loop self energy. Assuming that the term proportional 
to ff' is absent from the r-loop expansion Sen shows that it is also absent from 
the (r + l)-loop expansion. So the first term in Eq. ()A.5|) is the only possible 
structure of the fermion self-energy when its external momentum is jet like 
and approaches mass shell. 

Now let us investigate the 3 point functions. 
Fermion- gluon- fermion vertex function: F^, can depend on vectors that scale 
as I A, I'a Eq. ()A.2jl . provided all momenta external to the contracted vertex 
are coUinear to momentum given in Eq. ()2.2|) . It has one Lorentz index, /i, 
and neglecting the fermion masses, it contains an odd number of gamma ma- 
trices. This implies that the most general tensor and gamma matrix expansion 
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of involves 

2. Ip^yAVAl {U ■ Ia) and all permutations of 7^, ^, ^, 

3. yAi^A/il. hil/{iA-iA), VaI^aHIa-Ia), hill A- 

The differences between the listed set of structures and other possible combi- 
nations are C(A^/^ yi), as can be shown using Eqs. ()A.H) - ()A.2jl . The listed 
gamma matrix structures are multiplied by dimensionless functions, which can 
depend on the combinations l\, l\, l'^, l'^, besides the renormalization scale 
and the running coupling. Using the arguments similar to the ones leading to 
Eq. ()2.23|) . we easily verify that the above mentioned dimensionless functions 
are at most logarithmically divergent. Next we analyze the possible Dirac 
structures. 

1. The first case has the same structure as the elementary vertex, and 
therefore causes the same suppression as the elementary vertex. 

2. The fermion-gluon-fermion composite 3-point vertex is sandwiched be- 
tween the factors and /a, originating from the numerators of the 
fermion propagators external to the composite vertex. Therefore the 
terms from case 2 where J/a is on the first or third position in the string 
of the gamma matrices provide a suppression a/^- On the other hand in 
the case, when I/a is in the middle of this string of three gamma matrices, 
we encounter the combination ^AliiVA after taking into account the nu- 
merators of the external fermions. Using Eq. ()A.1|) . we can immediately 
recognize that this combination provides a suppression A^/^. 
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3. Based on the preceding arguments it is obvious that also the structures 
included in item 3 supply at least the same suppression factor as the 
elementary vertex. 

Therefore, we conclude that the composite 3-point fermion-gluon-fermion ver- 
tex behaves as the elementary vertex for the purposes of the jet power counting. 

Three gluon vertex: V^i^p, with external momenta collinear to momen- 
tum Pa. This vertex can depend on momenta Ia-, Ia defined above and the 
metric tensor g^p. Taking into account the dimension of the 3 gluon Green 
function, its only possible tensor structure involves combinations of the form 
[g^,v l^A + perm. + 0{\^l^ ^i)] and [/^ l^Jll + 0{\^/^ v^)], with all possible 
replacements of Ia — * Ia- These tensor structures are multiplied by dimension- 
less functions. The former is the same as in the case of an elementary vertex 
and it therefore supplies the same suppression factor as the elementary vertex. 
The latter also provides the same suppression as the elementary vertex, since 
the two momenta, say after being contracted with the propagators of 

the external gluons, give suppression factors, as in Eq. ()A.2j) . which cancel the 
1//^ enhancement. The leftover momentum provides the same suppression 
factor as the elementary vertex. Using the collinear power counting of Sec. 
12.3.21 one can immediately see that the IR divergence of the dimensionless 
functions multiplying these tensor structures is not worse than logarithmic. 
Hence, there is a suppression factor A^/^ associated with every contracted 3 
gluon vertex. 

Ghost - gluon - antighost three point vertex: When all lines external to the 
contracted vertex are of the order Ia-, the most general tensor structure for this 
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contracted vertex is 



l^A hi + l^A h2 + Oi\ 



(A.6) 



with dimensionless functions hi = hi jjp', Ja/^"^, as{^)), i = 1,2, which are 
at most logarithmically IR divergent. Using Eq. ()A.2|1 . we see that when the 
momenta in Eq. ()A.6|) are contracted with the tree level gluon propagator, we 
get a suppression of the order of the transverse jet momentum, and that this 
contracted vertex gives the same suppression as the elementary three point 
vertex, at least. 

A. 2 Varying the Gauge-Fixing Vector 

In this appendix we study the effect of an infinitesimal boost, performed on 
the gauge fixing vector r], on an expectation of a time ordered product of fields, 
denoted by O, taken between physical states. The gauge- fixing and the ghost 
terms in the QCD lagrangian are 



respectively. In Eq. ()A.7|) . 6 A is a Grassmann parameter defining the BRS 
transformation, ba{x) is an antighost field and 




•ghost 



ha{x) (5brs ga{x) /6A, 



(A.7) 



{d- {ri-d) 7]) ■ Aa{x). 



(A.8) 
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Let us consider an infinitesimal boost with velocity 6(3 on a gauge fixing vector 
7] performed in the plus-minus plane 



T] 1]' = 1] + T] Sp, (A. 9) 

where the vectors 1] and fi are defined in Eqs. (j2.9j) and (j2.35j) . respectively. 
Since only the gauge fixing and the ghost terms in the QCD lagrangian depend 
on 7], we can write to accuracy 0{6P^) 

6 <0> = < 0{r]') > - < 0{r]) > = < fj"^ 6(3 > 
= j d(^x < 0{r]) ga{x) 5ga{x) > 
- i j d^x< Oiji) baix) 5 ((5brs gaix)/5A) > . 

(A.IO) 

Using the BRS invariance of the QCD lagrangian and the BRS transformation 
rule for an antighost field 



SBRsba{x)/6A = ^Qaix), {A.ll] 



we arrive at 



6<0>=-tJd^x< ((5brsO/(5A) baix)6gaix) > . (A.12) 
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Taking a variation of ga{x) in Eq. ()A.12|) . we obtain 

5<0>=-i j d^x< ((5brsO/5A)6„(x) > . (A.13) 

Substituting for O a product of n gluon fields, we can use Eq. ()A.13|) . together 
with the rule for the BRS transformation of a gluon field 

5brs Al{x)/5K = d,c\x) + gJ-'^A'^{xY{x), (A.14) 

with c°'{x) representing the ghost field, to derive the gauge variation for a con- 
nected Green function. However, our jet functions are one-particle irreducible 
in external soft lines and we therefore cannot apply Eq. ()A.13|) directly, and 
must find an analog for this subset of diagrams. The modification of Eq. ()A.13jl 
due to the restriction to IPI diagrams is, however, not difficult to identify. 

Let us consider the graphical analog of the derivation of Eq. ()A.13j) just 
outlined. The variation in 77 may be implemented as a change in the gluon 
propagator and, in Coulomb gauge, the ghost-gluon interaction, which is also 
ry-dependent. This is the viewpoint that was taken in axial gauge in Ref. 
[44] . At lowest order in the variation, the modified gluon propagator produces 
scalar-polarized gluon lines, which decouple through repeated applications of 
tree-level Ward identities to the sum over all diagrams. The relevant tree-level 
identities are given in [46]. We need not describe these identities in detail here. 
We need only note that they are to be applied to any diagram in which a scalar 
polarized gluon appears at an internal vertex. Every such application produces 
a sum of diagrams, each of which fall into one of two sets: 1) diagrams in which 
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an internal gluon line is transformed to a yet another ghost line ending in a 
scalar polarization, and 2) diagrams in which one gluon line is contracted to 
a point. The new vertex formed in the former case is the ghost term, and in 
the latter case it is the ghost-gluon vertex of the BRS variation ()A.14|) . Eq. 
()A.13|) must result from the cancellation of all diagrams, set 2), in which an 
internal gluon line is contracted. Contracted external lines provide the ghost- 
gluon terms, and the ghost lines of set 1) eventually provide the ghost terms 
of the BRS variations ()A.14|) of external fields in Eq. ()A.13|1 . 

The simplicity of the tree level Ward identities puts strong limitations on 
the sets of diagrams that can combine to form different diagrammatic contri- 
butions to Eq. ()A.13|) . For diagrams of set 1), the topology of the original 
diagram is unchanged, and a IPI diagram remains IPI. For diagrams of set 
2), generally IPI diagrams remain IPI, except in the special case of a dia- 
gram that is two-particle reducible, with these two lines separated by a single 
propagator. In this case, the contraction of the internal line that separates 
the other two will bring those two lines together at a single vertex, producing 
a diagram precisely of the topology shown in Fig. 12.41 On the one hand, by 
Eq. (jA.13j) all such diagrams must cancel in the full perturbative sum. On 
the other hand, the same topology results from a diagram that is one-particle 
reducible with respect to a single line, which is then contracted as a result of 
the tree-level Ward identity. The latter diagram, however, is not included in 
the set of IPI diagrams with which we work. The application of the Ward 
identity to IPI diagrams only, therefore, results in terms that would cancel 
this special set of one-particle reducible diagrams, in which the only line that 
spoils irreducibility is contracted to a point. These are the diagrams shown 
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a) b) c) 

Figure A.l: Two-loop diagram illustrating the idea of the tulip-garden for- 
mahsm. Ti,T2,T3 are the possible tulips. 

in Fig. 12.41 in which the ghost-gluon vertex of Eq. ()A.14|) is inserted between 
one-particle irreducible subdiagrams in all possible ways. The ghost line end- 
ing at this composite vertex is continuously connected to the variation of a 
gluon propagator, according to Eq. ()A.13|1 . The full composite vertex of the 
Ward identity in Eq. ()A.13|) appears only at true external lines of the IPI jet. 
This vertex is given by the momentum factor in Eq. (j2.37j) and is represented 
by the double line crossing a gluon line in Fig. lA. 21 below. Diagrams that are 
reducible in one or more internal lines can be treated in a similar manner. The 
"left-over" terms in the Ward identities for each set of diagrams of definite re- 
ducibility properties (IPI, 2PI, etc.), must cancel in the full sum, reproducing 
the identity for Green functions, Eq. ()A.13|) . 

A. 3 Tulip- Garden Formalism 

In this appendix we illustrate how a given Feynman diagram contributing to 
the process (12.11) in the leading power can be systematically written in the 
form ()2.27p . For concreteness let us consider a two loop diagram where the 
quarks interact via the exchange of a one rung gluon ladder as in Fig. lA.ll 
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The important contributions of this diagram come from the regions when all 
of the exchanged gluons are soft, Fig. lA.lb or when the gluons attached to the 
A quark line are soft, while the rest of the gluons carries momenta parallel to 
the — direction (they belong to jet B), Fig. lA.lb . or when the two gluon lines 
attached to the B quark line are soft and the other gluons are coUinear to the + 
direction (they belong to jet A), Fig. lA.lb . The possible central soft exchange 
parts are called tulips. In our case the possible tulips are denoted as Ti, T2, T3 
in Fig. lA.ll The garden is defined as a nested set of tulips {Ti, . . . , T„} such 
that Ti C Ti+i for z = 1, . . . ,n - 1. In Fig. EH {Ti}, {T2}, {T3}, {Ti,T3}, 
{T2,T3} are the possible gardens. 

For a given tulip we make the soft approximation, consisting of attaching a 
soft gluon to jet A via the — component of its polarization only and to jet B via 
the the + component of its polarization. The result of this soft approximation 
for a given Feynman diagram F corresponding to a tulip T is denoted S{T)F. 
It has obviously the form of Eq. ()2.27j) . Following the prescription given in 
Refs. [44] and [16] we write the contribution to a given diagram F in the form 

F = Y.i-'^r^'SiT,) . . . 5(T„) F + Fr, (A.15) 

G 

where the sum over inequivalent gardens, as defined bellow, G in Eq. ()A.15|) 
is understood. The meaning of this expression is the following. For a given 
garden consisting of a set of tulips {Ti, . . . , T„}, we start with the largest tulip 
Tn and make the soft approximation for the gluon lines coming out of it. Then 
for Tn-i we proceed the same way as for T„. If some of the lines coming out 
of T^-i are identical to the ones coming out of T„ we leave them untouched. 
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For instance, if we consider a garden {T2,T3} from Fig. lA.ll we first perform 
the soft approximation on tulip T3 and then proceed to tuhp T2. However the 
hnes coming out of T2 and T3 which attach to the B quark hne are identical so 
when performing S{T2)S(T3)F we leave these gluon lines out of the game and 
make soft approximations only on the gluon lines attaching to the ladder's 
rung. Two gardens are equivalent if the soft approximation is identical for 
both of them. is defined by Eq. ()A.15j) . The contribution to comes 
from the integration region where \k\ > ^/s for all gluons coming out of the 
central soft part. As a result, the contribution to is suppressed by positive 
powers of -\/— Therefore we can ignore the contribution from Fr within 
the accuracy at which we are working. 
We can now rewrite Eq. ()A.15jl . as 

F = J2[ Yl i-ir^'SiT,)...SiT^^,)]s{T)F + Fn. (A.16) 

T \G,T„=T / 

This expression is indeed in the form of Eq. ()2.27j) since the term S{T)F is of 
that form and the subtractions ^(— l)"'+^S'(Ti) . . . S'(T„_i) modify only the 
soft function S in Eq. ()2.27|) . but do not alter the form of the equation. We 
can therefore conclude that the contribution to a given Feynman diagram in 
leading power can be expressed in the first factorized form given by Eq. ()2.27|) . 

A. 4 Feynman Rules 

In Fig. IA.21 we list the Feynman rules for the lines and the vertices encoun- 
tered in the text. The double lines are eikonal lines, while the dashed lines 
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i 

VA,B -k^ie 

—i 
k ■ k + ie 

— IV 

-i Saik) = -i [{ri ■ k)r]a + (r] ■ k)ria] 
Sap - kakp) 

~iSsfabc Onv 
^Qsfabc 9fiv 
dsfabc k^ 

Figure A. 2: Feynman rules for the eikonal lines, ghost lines and special ver- 
tices. 

represent ghosts. The four vectors f], fj are defined in Eqs. ()2.9|) and ()2.35p . 
respectively. The conventions for the gluon-ghost and gluon-eikonal vertices 
(third and second from the bottom of Fig. IA.2|) are the following. We start 
with a color index of a gluon external to the diagram defining the evolution 
kernel, see for instance Fig. VA.'Ah . then proceed to the gluon internal to the di- 
agram and finally to the ghost /eikonal line in order to assign the color indices 
of fabc- For the three point antighost - gluon - ghost vertex at the bottom of 
Fig. IA.21 we start with an antighost (arrow flowing out of the vertex) then 
proceed to the ghost and finally we reach the gluon line. 
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Pb 



PA 




q — k 



Figure A. 3: Two loop diagram demonstrating the origin of the Glauber 
(Coulomb) region. 

A. 5 Origin of Glauber Region 

In this appendix we exhibit the origin of the Glauber (Coulomb) region using 
the two-loop diagram shown in Fig. IA.31 Consider a situation when the upper 
gluon loop is a part of J a- Momentum k of the exchanged gluon flows through 
jet lines with momenta I2 = I — k and Z3 = pa — i — q + k. The components 
of k can be pinched by double poles coming from the denominators of the 
gluon propagators fc^ + ie and (g — k)"^ + ie. In addition to these pinches, the 
component k~ can be pinched by the singularities of the jet lines I2 and /s, at 
values 



The two poles given by Eq. ()A.17|1 are located in opposite half planes since 
in the region considered (2^3 > 0- This indicates that we must consider the 
possibility that the different components of the soft momentum k can scale 
differently. For instance, we can have fc"*" ~ /c_|_ ~ o^^/s and k~ ~ X^/s where 



k 



k 



r - 



r + 




(A. 17) 
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A ^ cr ^ 1. Indeed, the power counting performed in Sec. 12.3.21 shows 
that the singularities originating from these regions can produce a logarithmic 
enhancement. We also note that it is only minus components that are pinched 
in this way by the lines in J a-, and plus components by the lines in Jb- 
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Appendix B 



B.l Symmetry property of the jet functions 

Here we prove Eq. ()H.45|1 . In order to do so, we first shift tlie momentum 
of J^"'' in tlie transverse direction: Pa ^ Pa + Q- Tliis does not effect tlie 
logaritlimic beliavior of tlie jet function, but it makes easier tlie comparison 
of diagrams from and J^"^ 

The general diagram contributing to the jet function has the form 
shown in Fig. IB. 11 There are N gluons connecting the quark line to the rest 
of the diagram, which contains L loops. For every diagram contributing to 
we identify the corresponding diagram from J^""*. The guide is that they 
should be a mirror image of each other, but they should have the same color 
structure. For instance, in the case of a one loop contribution to J\ shown 
in Fig. IB. 2b the related diagram from J^^^ is the one in Fig. IB. 2b . 

Let us now make a transformation in {pA + q,"!]): 

kf — ^ —kf, ki± ^ ki± for all z = 1, . . . , n . (B-l) 
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Figure B.l: General class of diagrams contributing to J\ 



(n) 



In the denominators of the propagators this transformation can be compen- 
sated by a similar transformation of loop momenta if —if and lj± lj± 
for j = 1, . . . , {L + N — 1). The three-point vertices in the L-loop subdiagram 
generate a factor (—1)''^ under this set of transformations. Using the relation- 
ship between the number of loops, vertices and external lines on the L-loop 
subdiagram: 

vs = N + n + 2{L-V4-l), (B.2) 



we see that 



N+n 



(B.3) 



We also have to identify the relationship between the strings of gamma ma- 
trices on the fermion line 



(B.4) 



i=N-\ 
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Figure B.2: Identification of a particular one loop diagram between the jet 



functions J^"* and J^''. 



?(2) 



in J^"^ and 



J^a'^^ =v{pA + q,XA)lt,^ [ Yl {- ) '"(PaAi), (B.5) 



i=N-l 



in J^"\ Using the charge conjugation matrix C and its properties, Ref. [51]: 



i-l)'^'/'C^pVpip,X), 



(B.6) 



we obtain 

j("-°) = (-l)^(-ir-i(-l)^-+^^^(p^,Ai) (^n7.. ^,^l,^u{pA + q,XA). 

(B.7) 

In the Regge limit, we can neglect the momentum q in the numerators. Since 
all the momenta pi are coUinear to pa, Ja''^^ is symmetric under the exchange of 
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its external Lorentz indices. Therefore, up to the power-suppressed corrections, 
we have: 

J(n,0) ^ (_i)A^-l (_l)A^+Ai jinfl) ^g_g^ 

Combining Eqs. ()B.3|1 and ()B.8|) . we immediately recover Eq. ()3.45p . 

B.2 Evolution kernel 

In this appendix we provide the calculational details of the results given in 
Sec. EZl 

The contribution to Fig. 13.9b is: 

Fig.Hnb = 

X 
X 

Note that the overall sign in Eq. ()B.9|) reflects the minus sign explicit in Fig. 

We analyze the diagram in Fig. 13.9b in a similar way. It takes the form: 
Fig.ESfc = 

X 
X 
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Wl t facb ffbd fcfe 



d^k d^/ 



(27r)^ (27r) 



D 



/2 



N-P{k)S,{k) 



(q - If 
1 1 



{l-kf 



{k+ = 0) 



k"^ k - k k^ — itiq — kf 



0. 



(B.9) 



t facb ffbd fi 



cfe 



d^k d"/ 



{2n)D (27r) 
N-f'il).- -. N- + {q-l) 1 



J 



(2)de 



D -"A 



{pA,q-J~^ = 0,1 ,i±) 



N~P{k)Sp{k) 



A;2 



k ■ k k~ — ie 



:b.io) 



Note, again, that the overall sign in Eq. (lBlO|l reflects the minus sign explicit 
in Fig. IH.9b . The imaginary part of the eikonal propagator in Eq. (jB.lOj) gives 
vanishing contribution since the resulting integrand is an odd function under 
—k^. Hence only the principal value part contributes. In the Glauber 
region, ^ ~ /_|_, the integral vanishes due to the antisymmetry of the 
integrand under the transformation — > — k^ — * —k^. In the region 

(1') 

l± we can factor the gluon with momentum / from the jet J)^ . 
The contribution to Fig. 13.10b from a G gluon with momentum k corre- 
sponding to the first term in Eq. ()3.5()j) and to Fig. 13.8b is: 

/d^k d^l 

X ^"(pA,g;fc+ = 0,fc-,fc^,/+ = 0,r,/^), (B.ll) 
with the soft function Sq in Eq. ()B.11|) defined as 



iV-"(fc) S^{k) N-^q-k) 



Since, by construction k k^ ^ k± for an external G gluon with momentum 
fc, we can factor the gluon with momentum q — k from J)^ in Fig. I3.1()b . Then 
the resulting integrand in Eq. ()B.11|) : Scik^l^q) jk" , with Sq defined in Eq. 
flB.12|) : is an antisymmetric function under the transformation k^ —>■ —k^, 
—>■ — in the Glauber region /~ ^ ~ Z^. This indicates that the gluon 



with momentum / can be factored from . Therefore the contribution to 



r(3) 

A ■ 

Fig. I3.1()b from a G gluon agrees with a gluon reggeization at NLL. 
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The contribution to Fig. l3.1Ub is: 



(27r)^ (27r)^ ■^a'^^Pa^T, = 



P k-J {l-kf 



X Vux-iq-l,l-k,k-q) 



N^-{q-l) 



= 0) 



X 



N-''{k)S^{k) N- + {q-k) 
k"^ k ■ k {q — ky 



(B.13) 



In the Glauber region the integral vanishes due to the antisymmetry of the 
integrand. In the region ~ ~ the gluon with momentum / factors out 
and the contribution takes a factorized form. 

The diagram in Fig. I3.10l i for a G gluon with momentum k is: 



(Fig.Eim)G = -Wltfaed X j 



d^k d^l 



{27t)D (27r) 

(q - If '^"-^ V 
N~^{k) Sgjk) N- + {q-k) 
k"^ k -k (g — k)"^ 



D 



jf'\pA,q:l^ = ^,l-,lA) 



:b.i4) 



where 



bedc 
lip-iy 



+ 



fpbe fpdc i,9p.— 9pu 9pv9p—) ~l~ fpbd fpce i9pu9p— 9pp9v—) 

fpbc fped {.9pp9v- - 9p~9pu) , (B.15) 



is the Lorentz and the color part of the four point gluon vertex. In the Glauber 
region the integral over vanishes due to the antisymmetry of the integrand 
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Ill Eq. (pl4|l . In the region / ~ /+ ~ we can factor the gluon with 

(2) 

momentum / from the jet function J\ . Therefore the contribution from Fig. 
13. ion is in accordance with gluon reggeization. 

The contribution to Fig. l3.1Ub . from the second term in Eq. (|3.5U|) used in 
Eq. (pOnjl is: 



(Fig.Em)soft = git fade febc j -^-^S^;\k,l,q) 

X ^V,g;^+ = 0,A;-,A;^,/+ = 0,r,/^), (B.16) 



with the soft function S'i"^ defined by 



Vf,py{l, k,k I) ^f^^_ j^y^^ (f^^) 

k"^ k ■ k (g — ky 



[B.17) 



Eqs. fpi6|) and (PTfl) immediately follow form Eqs. (IHTTTl) and dEUll)- We 
apply the following identity 

s{i-,k~) = s{r = o,k- = o)e{M -\i-\)e{M -\k-\) 

+ [Sir, k~ = 0)- s{i- = o,k- = 0) e{M - \r\)] e{M - \k-'\) 

+ [S{1- = 0, A;-) - S{1- =0,k- = 0) e{M - |A;-|)] e{M - |/-|) 

+ [{s{i-, k-) - s{i-,k- = 0) e{M - \k-\)-} 

{s{i- = 0, k-) - s{i- = o,k- = 0) e{M - e{M - |r I)], 

(B.18) 
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to the soft function Ss"'\ Eq. ()B.17|) . to analyze the contribution to Eq. ()B.16|) . 
After using the first term of Eq. ()B.18|) for the soft function Si""^ in Eq. ()B.16|) . 
we arrive at 



(Fig.|Sm);;;t = git fade febcj 



X 



X 



(27r)^ (27r)« 

(^j dfc+ d/+ s^;\k^, /+, k- = 0,1- = 0, k±, q)^ 

rf''\pA,q;k^,h). (B.19) 



Performing the k~^ and Z+ integrals, we obtain 



d^-2A;x d^-2/^ (|A;x| 






\q± - k±- l±\) 


(27r)^-2 (27r)^-2 




2 


g_L - /c± - l±\ 


2 



X 



i\k±\ 




\q± - L 


1 + 


\q± - k± - l±\) i\k^\ + 


\l±\+2\, 


i± - 1±\ + \q± -k±- 1±\) 






\q± - 


1 


{\l±\ + \q±-l±\)H\k±\ 

-p(3) bed 


+ \l±\ + \q±-k^-h\) 
iPA,q;k^,h). (B.20) 


i\k±\ 


+ 


\q± - h\ 


+ 


\q±-k^-h\) ^ 



When analyzing the second (third) term in Eq. ()B.18jl . after used in Eq. 
(1111, we can factor the gluon with momentum I (k) from the jet function 
The resulting k~^ {k^, integral is over an antisymmetric function 
under the transformation —^^) k^ —k^ {k^ —k^, — > — /^) 

and therefore it vanishes. In the case of the last term in Eq. ()B.18|) . we can 
factor both gluons with momenta k and I from the jet J\ . So this part is in 
agreement with gluon reggeization. 

The contribution to Fig. 13.10b from the second term in Eq. ()3.50p used in 
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Eq. (jCTjl is: 



(Fig.lIIDb)soft = -gltfaebfced I -^-^Sf\k,l,q) 



{2ti)d (27r) 



Jf '"\PA, q; k+ = 0, k-,k^, 1+ = 0, r , h), 



(B.21) 



with the soft function Ss'^\k, I, q) defined by 



(fc_^)2 -l,l-q,q-k) -J- ^ ) (A;^ 



X 



/2 /.[ ■ 



(g - A;)2 

- /, / - g, g - fc) ) = 0) 



(g - fc)^ 



(B.22) 



We apply the identity ()B.18|) to the soft function Si'^\ Eq. ()B.22|) . to analyze 
the contribution to Eq. ()B.21j) . After using the first term of Eq. ()B.18j) in Eq. 
(lB.21jl . we arrive at a contribution 



(Fig.iiiiib);;i 



-ds t faeb f, 



d^-2A;, d^-2Z 



(27r)^ (27r)^ 
j dA;+ d/+ Sf^ ik+, t,k~ = 0, l 



0,k±,l±,q) 



X r 



(3) bed 



(B.23) 
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Performing the and l'^ integrals, we obtain 



{27r)D (27r) 



D 



X 





+ \l± 


+ \q± - 


k± 


+ \q±- 


k±-l±\) 


\k±\\h 






h 




+ \q± - kx 


+ \q±-k±-l±\) 





'-\h\i\qi 




2 + 


q±-k±^- /_L 


2) _ 2\q^ _ 


'\q±-k^-h\) 


{\k^ 




+ \q±-k^-h\] 


H\k± 


+ \q± - 


k±\r 



y<rf''\pA.q;k^,L). (B.24) 



When analyzing the second (third) term in Eq. ()B.18jl . when applied to Eq. 
()B.2H1 . we can factor the gluon with momentum I (k) from the jet function 
jjf"*. The resulting {k^, integral is over a function antisymmetric 

under the transformation —i^, k~^ —>■ —k^ {k^ ~k^, — 

and therefore it vanishes. In the case of the last term in Eq. ()B.18|) . we can 
factor both gluons with momenta k and / from the jet J\ . So this part is in 
agreement with the gluon reggeization. 

The contribution to Fig. 13.10b from the second term in Eq. ()3.50p used in 
Eq. (nnn|l is: 

(27r)^ (27r)^ ^a^'^'Ipa, g; = 0,/",/^) 
X 1^— ^ -k, k - I) ^ KA-(g -lj-k,k-q) -J^^Jy^) (k^) 

-k, k - ^(7!^~,^V ..„(g -i,i-k,k-q) ik^ = 0) 

N-''{k)Sa{k) N- + {q-k) 
k"^ k ■ k {q — ky 

(B.25) 
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This follows immediately from Eq. ()B.13|) . In the Glauber region, ^ ~ 
l±, the integrand is an antisymmetric function of and and therefore the 
k^, integral vanishes. In the region ~ ~ /_|_, we can factor the gluon 

(2) 

with momentum / from the jet function Jj^ . Hence this contribution is in 
agreement with gluon reggeization. 

The contribution to Fig. 13.10b from the third term in Eq. ()3.50|) used in 
Eq. (Kmf is: 



(Fig. liiJilb)soft' = g'^t fade febc J ^2ti)D (27r)^ '^''''^ 

X '"'(pa, q; k+ = 0,k-, k^, /+ = 0, /x), 

(B.26) 

with the soft function S^'^^ defined by 



S^^ ^ ^:L^(V,,.(l,-k,k-l)-g+,V,^,il,-k,k-l))^'^'^^ 



2 



NP-{k) N- + {q-k) 

This follows immediately from Eqs. ()B.11|) and ()B.12|) . We apply the identity 
fIB.lSj) to the soft function S^^\ Eq. ()B.27|1 . to analyze the contribution to Eq. 
()B.26|) . After using the first term of Eq. ()B.18|) in Eq. ()B.26|) , we arrive at 



(Fig.Em)St' = 9s t fade febc j 



d^-2^ , d^-2/ 



(277)^ (27r)^ 

X (^j dk+dl+ S''f{k-^,l+,k~ = 0,1- = 0,k^,l^,q)^ 
X rf''\pA,q;k^,h). (B.28) 
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The integrand S^V in Eq. (in:28|) is an antisymmetric function of /+ and 
and therefore the integral over and in Eq. ()B.28|) vanishes. In order to 
show this antisymmetry, we have used that the component 

Ar+"(A;- = 0, k+, k±) Sa{k~ = 0, k+, k±) = 0, (B.29) 

in the Glauber region, fc~ ~ 0. We have also used the Ward identity 

-k, k-l){k- ly = 0, (B.30) 

when analyzing the contribution from the transverse polarization of the gluon 
with momentum k in Fig. IH.lOb . When studying the second (third) term in Eq. 
(IRT8|) . after applied to Eq. (jR26l), we can factor the gluon with momentum 
I [k) from the jet function J^^\ The resulting k~^ {k^, integral is over 
an antisymmetric function under the transformation — > — k~^ —k^ 
{k^ —k^, ~^^) S'lid therefore it vanishes. In the case of the last term 

in Eq. ()B.18j) . we can factor both gluons with momenta k and / from the jet 
Therefore the contribution from fFig. l3.l0k )„„ft'. Eq. ()B.26|) . is in an 
agreement with gluon reggeization at NLL. 

The contribution to Fig. IH.lOb from the third term in Eq. ()3.5()j) used in 
Eq. (pOnjl is: 

(Fig. |3JU|D),oft' = -git faeb feed j -^^^ 

X Jf ''\pA, q; k+ = 0,k~,k±,l+ = 0, (B.31) 
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with the soft function S'y ^ defined by 



•^^^ = ^(fc ^-ly^ (^^^^(^ -iJ-q,q-k)- g+pV^^^{k - l,l - q,q - k)) 

This follows from Eqs. (IK;2T|1 and (1^:2211 . We apply the identity (IHlSll to the 
soft function, Eq. ()B.32|) . to analyze the contribution to Eq. ()B.31|) . After 
using the first term of Eq. ()B.18|) in Eq. ()B.31|) , we arrive at 



(Fig.mnh)!! = -g'tfaebfced J 



X 



(27r)^ (27r)^ 

(^j dk+ dz+ s^\k+, /+, k- = 0, r = 0, k^, h, q)^ 

Tf''''{pA,q;k^,h). (B.33) 



The integrand S^^'' in Eq. (in33|) is an antisymmetric function of /"*" and k~^ 
and therefore the integral over k~^ and in Eq. ()B.33|) vanishes. In order to 
show this antisymmetry, we have used that the component N^^{q — /) = in 
the Glauber region, /~ ~ 0. We have also used the Ward identity 

{k - ini - qyV^p.{k -l,l-q,q-k){q-ky = 0, (B.34) 

when analyzing the contribution from the transverse polarization of the gluon 
with momentum g — / in Fig. 13.10b . When analyzing the second (third) term 
in Eq. ()B.18j) . applied to Eq. ()B.3H) . we can factor the gluon with momentum 
/ (k) out the jet function J^^\ The resulting k^ (A;^, /+) integral is over 
an antisymmetric function under the transformation k"^ —k^ 
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{k^ —k^, ~^^) S'lid therefore it vanishes. In the case of the last term 

in Eq. ()B.18j) . we can factor both gluons with momenta k and / from the jet 

. Therefore the contribution from (Fig. EHb) soft', Eq. jEM}, is in an 
agreement with gluon reggeization at NLL. 

The contribution to Fig. 13.10b from the third term in Eq. ()3.50|) used in 
Eq. (pOnjl is: 



d^k d^l 



(27r)^ (27r) 



D 



X Ja iPA,q;r = 0,1 

X {V,^^{1, -k, k-l)- -k, k - 0) ^ ^ 

X (KA7(g - 1,1 - k,k - q) - g^^Vu\-{q - 1,1 - k,k - q)) 
N^'iq - I) NP^{k) Sa{k) iV^ + (g - k) 
{q-iy k-k {q-ky ' 



B.35) 



In the Glauber region, ^ ~ the integrand is an antisymmetric 
function of k^ and l'^ and therefore the k^, integral vanishes. In the region 
~ /+ ~ /_L, we can factor the gluon with momentum / out the jet function 

(2) 

J\ . Hence this contribution is in agreement with gluon reggeization. 

The diagram in Fig. I3.10l i from the third term in Eq. ()3.50p used in Eq. 
(ICTll is: 

(Fig.Einii)soft' = -igltU,,X ^>^,g;/+ = 0,r,Z^) 

S.jk) N^^q-k) 

k-k {q-kf ' ^^-^^^ 
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In the Glauber region, ^ ~ the integrand is an antisymmetric 
function of and and therefore the k^, integral vanishes. In the region 
~ ~ we can factor the gluon with momentum / out the jet function 

(2) 

J\ . Hence this contribution is in agreement with gluon reggeization. 
The contribution from the first term in Eq. ()3.5H) is: 

(27r)^ (27r)^ 'S'jet(fc, ?) 
X ''\PA, q; k+ = 0,k-, k^, 1+ = 0,1-, h), 

(B.37) 

with a soft function S^et defined 



Sictik,l,q) 



F k-k {q-l-kf (l-q).(l-q) P ' 

(B.3^ 



Next we use an identity ()B.18|1 for Sjct{k ,1 ), defined in Eq. ()B.38j) . After 
performing the /c"*" and Z"*" integrals in the first term of Eq. ()B.18jl , when used 
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in Eq. ()B.37|) . we obtain the following contribution: 



X 





+ \l± 


+ \q± - h 


+ 


2\q± -k^- 




\k±\ \q± - 
i\k± 


k±-l±\{\l± 
'-\q±-h 


\ + \q± - l±\) {\k± 
^ + \q±-k±- h 


+ \q± - k± 


+ \q±-k^-h\y 



{\k±\ + \q± - l±\ + \q± -k±- h\y \q^ -k^-lj 

+ 



\l± 


+ 2|g±- 




l± 






+ \q± -k_i- 


l±\)i\k± 




+ \q±-k^~h\)^ 



\q± - h 


+ 2\q±-k^-h 




{\q±-L 


+ \q±-k^-h\) {\k^ 


+ \q±-i±\ 


+ \q±-k^-h\r 



i\k± 


2 _ 


- \q± - l± 


2 _ 


■ \q± -k±- /_Lp 


') 


\k±\i\l± 


'-\q±- 


l± 


^) \q± -k^-h 


3 



(B.39) 



In the second (third) term of Eq. ()B.18|) . we can factor the gluon with momen- 
tum I (k) from the jet function J^^\ The resulting k^ {k^, /+) integrals 
are over antisymmetric functions and therefore they vanish. In the last term 
of Eq. ()B.18|) . we can factor both gluons with momenta k and / from the jet 
function Hence this contribution is in a factorized form, in agreement 

with the gluon reggeization. 

The contribution to the second term in Eq. ()3.5ip comes from a diagram 
in Fig. imik : 

{27r)D {2-k)d '^A^'^\PA,q] = 0,l-,l±) 

X Vp(? ^' k,k + l q) 

Nf^[q-l-k) jl-q), iV- + (0 
{q-l-ky {l-q).(l-q) P ' ^ ' ^ 

In the Glauber region l±, the integrand is an antisymmetric function 
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under the transformation — > —k^ and — so the integral over 

k^ and vanishes. In the region /~ ~ /"•" ~ lj_ we can factor the gluon 

(2) 

with momentum I from the jet function J\ using the K-G decomposition. 
Therefore the contribution from Fig. l3.1Ub takes a factorized form. 

The same reasoning apphes to the diagram in Fig. I3.10t '. which is due to 
the third term in Eq. 1)3.511) : 

/d^k d^l 
(27r)^ (27r)^ •>^a^'^'(pa, g; k^ = 0,k~,k±) 



P k-k {q-kf ^"'^^ ' ' ' P 

N'^^jq-l-k) jl-q)^ 
{q-l-kf {l-q).{l-qy ^ ' ^ 



It also factorizes in agreement with the gluon reggeization hypothesis. 
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Appendix C 



C.l Eikonal Example 

In this appendix, we give details of the calculation of the logarithmic behavior 
in the diagrams of Fig. 14.31 We choose the reference frame such that the 
momenta of the final state particles are given by: 

A = (1,0,0,1), 

(52 = (1,0,0,-1), 

/ = u;i(l,sz,0, q), 

k = cjfc(l,SfcCos0,SfcSin0, Cfc). (C.l) 

Here we define = sin6';fc and q^^ = cos6i^k. 6i is the angle between the 
vectors / and f3i, 9^ is the angle between the vectors k and j3i and is the 
azimuthal angle of the gluon with momentum k relative to the plane defined 
by (32 and I. The available phase space in polar angle for the radiated 
gluons is Ok e in/2 - 5, n/2 + 5) and 9i E (0, n/2 - 5) U {n/2 + 6, n). 

Using the diagrammatic rules for eikonal lines and vertices, as listed for ex- 
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ample in [8] , we can write down the expressions corresponding to each diagram 
separately. For example, diagram 14.31 a) gives 



a) + (k^l) = [Uc^r{tMtc)][-^gtPlP2Pl) V^^^{k + I, -k, -I) 



pi- (k + l) 2k-l (3^-1 (32-k 
+ {k^l). (C.2) 

Va/3^{k+l, —k, —I) = [{2k+l)jga/3 + {l — k)ag/3'y — {2l + k)/3ga'y] is the momentum- 
dependent part of the three gluon vertex. Using the color identity fabc^T^(tatbtc) '- 
iCpAfcCA/2, and the approximation Pj ■ I ^ (3j ■ k for j = 1, 2, which is valid 
due to the strong ordering of the final state gluon energies, we arrive at 

a,Uk»l)^\c,McC.^-l^[-^-±^^^^-^^. (C,3) 

We proceed in a similar manner for the rest of the diagrams. The results are: 

1^ *r ^ 1 1 



c) = -CpMcCAg. 



4 " " k-l (3i-l(32-k- 

1^ ^ 4A-/?2 1 1 



e) = CFUc{CF-CAl2)g, 



APl-P2? 1 



+ ^ C.^c(C.-C./2)s^|l^^^-^. (C.4) 



The color factors in the last two equations of ()C.4|) are obtained from the 
identity Tiitahtatb) = CpAfciCp — Ca/2). Combining the terms proportional 
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to the color factor CfMcCa-, and including the complex conjugate diagrams, 
we find for the squared amplitude 



\Mf = 2gtCFMcCAPi-l32 

X ( 1 I 1 /^i-/^^ ^ 

\k-l(3i-k(32-l k-l(3i-l(32-k ■ I p2 ■ I Pi ■ k p2 ■ k J ' 

(C.5) 

Having determined the amplitude, we need to integrate |Mp over the phase 
space corresponding to the geometry given in Fig. 14.21 Specifically, we have 
to evaluate: 

(C.6) 

where the weight function /(/, a) is given, as in Eqs. ()4.4p and (HHH), by 

1^(1- CiY-'' s? : 01 e (0, n/2 - 6) 
f{l,a) = \ ' (C.7) 

with a < 1. 

Using the equalities: /3i ■ = 2, /5i ■ Z = uji{l — q), [32 ■ I = uji(l + q), 
Pi - k = Ukil - Ck), P2 - k = uJk{l + Ck) and k - l = UkU^iil - CkQ - SfcS;cos0) 
in Eq. ()C5|) . performing the integration over 0, and changing the integration 
variable q — > — q in the angular region 6i E {n/2 + 6, vr), we easily arrive at 
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the following three-dimensional integral: 



/ — CfCa 
1 



sin (5 fl 

dcfc / dQ / — e 

sin 5 J sin 5 J £\/i 

1 1 \ 1 1 

+ 



Cfc + Q 1 + Cfe V 1 + Q ' i- CkJ sl 1 + Ci 



(C.8) 



We are interested in the (l/£)ln(l/£) behavior of /. This is obtained after 
performing the Ui integral with the replacement e~^^'-^^~'^'-^^ "■sf/^/s _^ q^^^ _ 
vuji{l — ciY'"" si / Remainders do not contain terms proportional to Ins. 
In this approximation, the ci integration can be carried out, and we obtain the 

integral representation for the term containing {1/e) ln(l/£): 



si — cos^ 5 



-In 



1 + sin 5 



In 



1 + sin 5 
1 — sin 5 



(C.9) 



The potential non-global logarithm of e is replaced by hi{ei/). The angular 
integral over can be expressed in terms of dilogarithmic functions. The 
final expression for the term proportional to ln(ez/)/e takes the form: 



\ TT / £ \ev J 



tt"^ In ( ^ (1 + sin 5) \ / 1 + sin S 
"6"^^V 4 ) ^ U - sin5 



+ Li2 



1 — sill S 



- Lio 



1 + sill S 



- Lio 



2 sin 5 
1 — sin (5 



- Lio 



1 — sin d 



1 + smS 



(C.IO) 



Equivalently, we can express our results in terms of the rapidity width of the 
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region f2, Eq. ()4.14|) . and we obtain 



IT / e yeu 



y + A77 - In (2 sinh(A77)) 



+ Li2 — r- — -— - Li2 



2 cosh(Ar//2) J \2 cosh(Ar//2)^ 

Li2 (-2 sinh(A?7/2) e^"/') - Liale"^")] . (C.ll) 



The coefficient 



\ 2 T 



C(Ar/)^- - ^ ^ , , . (C.12) 



as a function of A77 is shown in Fig. 1(111 Naturally, C is a monotonically 
increasing function of At]. For Arj 0, 

C~ 0(Ar7 InAr^), (C.13) 

and the cross section vanishes, as expected. On the other hand, as the size of 
region fl increases, C rapidly saturates and reaches its limiting value [68] 



71^ 

lim C=—. (C.14) 



C.2 Recoil 

In this appendix, we return to the justification of the technical step represented 
by Eq. ()4.23|) . According to this approximation, we may compute the jet 
functions by identifying axes that depend only upon particles in the final 
states Nj^ associated with those functions, rather than the full final state N. 
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Figure C.l: C{Ari), as defined in ()C.12|) . as a function of rapidity width Ar] 
of the region fl. The dashed fine is its hmiting value, C{Ar} —>■ oo) = 

Intuitively, this is a reasonable estimate, given that the jet axis should be 
determined by a set of energetic, nearly collinear particles. When we make 
this replacement, however, the contributions to the event shape from energetic 
particles near the jet axis may change. This change is neglected in going from 
the original factorization, Eq. ()4.2H) . to the factorization in convolution form, 
Eq. ()4.27j] . which is the starting point for the resummation techniques that we 
employ in this paper. The weight functions /^(A^j, a) in Eq. 1)4.2111 are defined 
relative to the unit vector ni corresponding to a = 0, the thrust-like event 
shape. The factorization of Eq. ()4.2H) applies to any a < 2, but as indicated 
by the superscript, individual contributions to /^(A^j,a) on the right-hand 
side continue to depend on the full final state A^, through the identification of 
the jet axis. 

To derive the factorization of Eq. ()4.27|) in a simple convolution form, we 
must be able to treat the thrust cLXlS, SI'S cL fixed vector for each of the 
states Ng, Nj^. This is possible if we can neglect the effects of recoil from soft. 
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wide-angle radiation on the direction of the axis. Specifically, we must be able 
to make the replacement 



fl{Nj^,a)^MNj^,a), 



(C.15) 



where fc{Nj^,a) is the event shape variable for jet c, in which the axis fic is 
specified by state Nj^ only. Of course, this replacement changes the value of 
the weight, e, {Nj^, a) 7^ fc{Nj^, a). As we now show, the error induced by 
this replacement is suppressed by a power of e so long as a < 1. In general, the 
error is nonnegligible for a > 1. The importance of recoil for jet broadening, 
at a = 1, was pointed out in [59]. We now discuss how the neglect of such 
radiation affects the jet axis (always determined from a = 0) and hence the 
value of the event shape for arbitrary a < 2. 

The jet axis is found by minimizing f{a = 0) in each state. The largest 
influence on the axis "he for jet c is, of course, the set of fast, collinear particles 
within the state Nj^ associated with the jet function in Eq. 1)4.211) . Soft, wide- 
angle radiation, however, does affect the precise direction of the axis. This is 
what we mean by 'recoil'. 

Let us denote by Us the energy of the soft wide-angle radiation that is 
neglected in the factorization ()4.27p . Neglecting this soft radiation in the 
determination of the jet axis will result in an axis ni(A^j^), which differs from 
the axis ni{N) determined from the complete final state (A^) by an angle As4>: 



At the same time, the soft, wide-angle radiation also contributes to the total 



Q 



(C.I6) 
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event shape f{N,a) ~ {l/Q)k1{k-y-'' at the level of 




(C.17) 



because for such wide-angle radiation, we may take ~ ks^± ~ t^s- In 
summary, the neglect of wide-angle soft radiation rotates the jet axis by an 
angle that is of the order of the contribution of the same soft radiation to the 
event shape. 

In the factorization ()4.27j) . the contribution of each final-state particle is 
taken into account, just as in Eq. ()4.21|) . The question we must answer is how 
the rotation of the jet axis affects these contributions, and hence the value of 
the event shape. 

For a wide-angle particle, the rotation of the jet axis by an angle of order 
As4> in Eq- ()(^.16jl leads to a negligible change in its contributions to the 
event shape, because its angle to the axis is a number of order unity, and 
the jet axis is rotated only by an angle of order Eg- Contributions from soft 
radiation are therefore stable under the approximation ()4.23|1 . The only source 
of large corrections is then associated with energetic jet radiation, because 
these particles are nearly coUinear to the jet axis. 

It is easy to see from the form of the shape function in terms of angles, Eq. 
()4.11|) . that for any value of parameter a, a particle of energy Ui at a small 
angle 6i to the jet axis ni{N) contributes to the event shape at the level 




2~a 
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The rotation of the jet axis by the angle A^^ due to neglect of soft radiation 
may be as large as, or larger than, 9i. Assuming the latter, we find a shift in 
the Si of order 

(C.19) 

The change in Ei is thus suppressed by at least a factor e^^"" compared to 
which is the contribution of the wide-angle soft radiation to the event shape. 
The contributions of nearly-collinear, energetic radiation to the event shape 
thus change significantly under the replacement ()4.23|) . but so long as a < 1, 
these contributions are power-suppressed in the value of the event shape, both 
before and after the approximation that leads to a rotation of the axis. For 
this reason, when a < 1 (and only when a < 1), the value of the event shape 
is stable whether or not we include soft radiation in the determination of the 
jet axes, up to corrections that are suppressed by a power of the event shape. 
In this case, the transition from Eq. ()4.2ip to Eq. ()4.27p is justified. 
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